实验5.2

抛物线法求解积分。

clear
clc
syms x;
a = zeros(1,6);
b = zeros(1,6);
disp(vpa(int(exp(-x.^2),1,2)));
a(1) = 0.1;
for i =2:6
    a(i) = vpa(a(i-1)*0.1,15);
end
for i =1:6
    b(i) = vpa(quad('exp(-x.^2)',1,2,a(i)),10);
end
disp(vpa(a),10);
disp(vpa(b),15);
plot(log(a),b);

数值一直在0.1352左右。

右边太长不截了，大概明白就好了。