Write a Python code to estimate the integral of the following function from a=input, b=input with n=input using Simpson's 1/3 method. Don't use import fx=0.2+25x-200x^2+675x^3-900x^4+400x numpy.
Write the program to estimate the integration using Eq.(5.1) and Eq.(5.2)
(5.1)
∫[a,b] f(x)dx ≈ (1/3n) [f(x0) + 4∑f(xi) + 2∑f(xj) + f(xn)] where i=1,3,5 and j=2,4,6, n=3n
(5.2)