Question

方程求解：如何用matlab r2018中for循环语句与解方程结合解决二维数据？

Ka1=1.3E-7; Ka2=7.1E-15;Ksp1=6.3E-36; Ksp2=1.07E-21;a=1;for pH=0.5:0.1:4; H=10.^-pH; b=1; for M=0.5:0.05:1.5; syms c(a,b);c(a,b)=solve(c.*pH+M,c);c(a,b)=double(c);index(a,b)=find(c>0);R(a,b)=c(index);b=b+1;end a=a+1;end

Answer 1

次求得的解应该是相同的制大小的一维数组，因为不知道具体的解的维度，大概写一下，应百该像下面度这样吧：

for i=1:20

x(:,i)=fsolve(@energe,x0,[],Qxi(i),Qyi(i),cx(i),cy(i));

end

注意要把i＝i+1去掉，for循环会自动问加的。。。

如果不对的话请把具体维数告诉我，才能写出程答序～～