用复化梯形公式、复化Simpson公式计算下列积分中的一个,要求误差不超过:
(1)对lnx在1到10上求定积分
(2)对1/exp(x)在0到1上求定积分
我编的程序是
>> a=1;
>> b=10;
>> quad('log(x)',a,b)
ans =
14.0259
>> clear
>> format long
>> a=1;
>> b=10;
>> m=0;
>> n=1000;
>> h=(b-a)/n;
>> f=inline('log(x)');
>> for i=1:(n-1)
x=0;
x=a+h*i;
m=m+f(x);
end
>> th=h/2*(f(a)+2*m+f(b));
>> th
th =
14.025844854958663
我的问题是,就在我这个方法的基础上,如果不给n的值,用程序来算出满足误差小于10^(-5)下的n值呢?