Simpson's three by eight rule

Simpson's three by eight rule

Code is Here

clear;clc
fxn=input("Enter the fuction : (eg:y=(x^2+5*x+6)) : ")
deff('y=f(x)',fxn);
printf("The given function is %s\n\n",fxn);
a=input('Inter the lower limit : ')
b=input('Inter the upper limit : ')
n=input('Inter the value of n : ')
h=(b-a)/n
for i=0:n
    X(1,i+1)=a+h*i
    Y(1,i+1)=f(a+h*i)
end
printf("x=")
disp(X)
printf("y=")
disp(Y)
a1=Y(1,1)+Y(1,n+1);
a2=0
for j=3:3:n-1
    a2=a2+Y(1,j+1);
end
a3=0
for k=0:n-1
    p=modulo(k,3)
if p==0
else
    a3=a3+Y(1,k+1);
end
end
i=(3*h/8)*(a1+2*a2+3*a3);
printf("The value of Integral is : ")
disp(i)

Simpson's three by eight rule