1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33
|
function [number,error]=assignment2(a,b,e,n,t,number,error)
if n==1 y=(sin(b)/b+sin(a)/a)*(b-a)/2; else y=t/2; q=(b-a)/n; for ii=1:n/2 y=y+q*sin(a+(2*ii-1)*(b-a)/n)/(a+(2*ii-1)*(b-a)/n); end number=[number n]; error=[error abs(y-t)/3]; end if abs(y-t)<e*3 else [number, error]=assignment2(a,b,e,n*2,y,number,error); end
|