我用matlab学习梯形法求定积分,把要求的函数定义为f1.m
function y=f1(x)
y=exp(-x*x);
我在另一个M文件ts.m中计算f1在0-1上面的定积分:
a=0;
b=1;
n=1;
h=b-a;
eps=0.00001;
syms x
I1=h*(f1(a)+f1(b));
while (n<1000)
for k=0:n-1
I=I1/2+h/2*f1(a+k*h+h/2);
end
if abs(I-I1)<eps
break
end
n=n+n;
h=h/2;
I1=I;
end
I
为什么报错