%1.算法主函数如下:
function y=paiabc(Eps);
y =3; i=0;j=1;
while (abs(y-pi)>=Eps);
j=j+1;
y=y+(-1)^j*4/((i+2)*(i+3)*(i+4));
i=i+2;
end;
%调用函数的代码:
format long g;
a=16;
x=1./logspace(1,a,a),
line(1:length(x),ones(1,a)*pi);
axis([1 a 3.13 3.17]);
xlabel('点号');
ylabel(['计算值与','{\pi}','值']);
title('计算精度示意图');
hold on;
for k=1:a;
y=paiabc(x(k));
pause(0.5);
plot(k,y,'r.');
end;
legend('精确值','计算值');
y,
%结果如下:
%2.算法主函数如下:
function y=paiabcd(N);
j=0;
for i=1:N;
r=rand(1,2);
if r(1)^2+r(2)^2<=1;
j=j+1;
end;
end;
y=4*j/N;
%调用函数的代码:
format long g;
a=7;
x=logspace(1,a,a),
line(1:length(x),ones(1,a)*pi);
axis([1 a 2.5 3.9]);
xlabel('点号');
ylabel(['计算值与','{\pi}','值']);
title('计算精度示意图');
hold on;
for k=1:a;
y=paiabcd(x(k));
pause(0.5);
plot(k,y,'r.');
end;
legend('精确值','计算值');
y,
%结果如下:
迈特莱博
2024-08-07 广告