Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function I=trapez(a,b,c)
- a=0;
- b=1;
- h=(b-a)/2
- f=@(x) 2./(1+x.^2);
- fp2=@(x) 16*x.^2/((x.^2+1).^3)-4/((x.^2+1).^2);
- r=feval(fp2,c)./(12).*(b-a).^3
- I=h*(f(a)+f(b))-r
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement