主函数:
clc
clear
close all
tic
a=1; b=2; c=3;
for i = 1 : 30
x0 = fsolve(@(x) myfun(x,a,b,c),0);
a = 1 + x0; b=2*(b+x0); c=3*(c+x0);
X(i) = x0;
end
plot(X)
toc
子函数:
function F = myfun(x,a,b,c)
F = a - b*exp(c*x);
运行结果:
[attach]257864[/attach]