Matlab利用循环求数值解（wzl）

下面展示 代码。

%循环结构
s=0;
for n=1:10
    for m=1:10
        s=s+m*n;
    end
end
disp(s)
%输出结果：3025

下面展示 代码。

% for循环求解
x0=0;
for k=1:100
	x=x0;  %赋给x原始值x0
	x0=cos(x);  %再用x得cos(x)值赋给x0
end  %x=x0=cos(x)，利用中间值求解
>>disp(x)
% 输出结果：0.7391

%while循环求解
x=0;
y=cos(x);
while abs(x-y)>5e-5  %只要大于就继续循环
	x=y;
	y=cos(x);       %利用中间值y
end                 %直到x与cos(x)非常接近
>>disp(x)
%输出结果：0.739


%利用画图求解
y1=x;
y2=cos(x);
x=0:0.1:pi/2;
>>plot(x,y1,'r--')
>>hold on
>>plot(x,y2,'b--')

输出结果