求助各位matlab大侠们~
请教大侠们,我下面的这个函数哪里错了?调用的时候,如r(0,4)总是说自变量太多。我是新手,matlab语言不懂啊~ function y=intf(x)
y=exp(-x.^2);
function [s,n]=r(a,b,eps)
if nargin<3
eps=1e-4;
end
s=10;
s0=0;
k=2;
t(1,1)=(b-a)*(intf(a)+intf(b))/2;
while (abs(s-s0)>eps)
h=(b-a)/2^(k-1);
w=0;
if(h~=0)
for i=1:(2^(k-1)-1)
w=w+intf(a+i*h);
end
t(k,1)=h*(intf(a)/2+w+intf(b)/2);
for p=2:k
for i=1:(k-p+1)
t(i,p)=(4^(p-1)*t(i+1,p-1)-t(i,p-1))/(4^(p-1)-1);
end
end
s=t(1,k);
s0=(t(1,k-1));
k=k+1;
n=k;
else s=s0;
n=-k;
end
end
r.rar
(703 Bytes)