Matlab fsolve和for连用,关于利用for循环和fsolve函数求解一系列五次方程根的疑问（每次得到的结果都是初值）...

format long;

r1=1e-3;    % 圆心到器件外界面层1上表面的距离

hc=1e-4;    %陶瓷板厚度

pai_p=pi/6;    % 热电支p腿中心角

delta=1e-3;  % 热电支腿厚度

T1=400;   % 热端温度

T2=300;   % 冷端温度

alpha_p=2e-4;   % P型热电材料seebeck系数

alpha_n=-2.5e-4;   % n型热电材料seebeck系数

%alpha_n=-2.13e-4;   % n型热电材料seebeck系数

alpha=alpha_p-alpha_n;   % 热电器件seebeck系数

kp=1.6;  % p型热电材料热导率

kn=0.85;  % n型热电材料热导率

%kn=0.7;  % n型热电材料热导率

sigma_p=1.1e+5;  % p型热电材料电导率   Bi2Te3

sigma_n=0.38e+5;  % n型热电材料电导率   Bi0.5Sb1.5Te3

%sigma_n=0.277e+5;

r=0.4;   %r=kp/kc

kc=kp/r;

n=1e-4;    %n=sigma_p/sigma_c

sigma_c=sigma_p/n;

er=0.1:0.01:5;   %er=theta_n/theta_p

%er=[0.5,1,2];   %er=theta_n/theta_p

sr=1.05:0.05:5;     %sr&