na-lab-simpsons

code

clc;

close all;

f=input("Enter function");

a=input("Enter lower limit");

b=input("Enter upper limit");

n=input("Enter no of subintervals");

h=(b-a)/n;

sum=0;

for i=1:n-1

x=a+h*i;

if rem(i,2)==0

sum=sum+2*f(x);

else

sum=sum+4*f(x);

end

end

sum=sum+f(a)+f(b);

sum=sum*(h/3)