def f(x): return 400*x**5-900*x**4+675*x**3-200*x**2+25*x+0.2

Ireal=1.640533
a=0.0
b=0.8
n=4


sum1=0
sum2=0

h=(b-a)/n

for i in range(1,n,2):
    x=a+i*h        # x=0.2 and x=0.6
    sum1=sum1+f(x)
for i in range (2,n-1,2):
    x=a+i*h         # x=0.4  
    sum2=sum2+f(x) 
I=h*(f(a)+4*sum1+2*sum2+f(b))/3
print(I)

Et=abs((Ireal-I)/Ireal)*100
print("Et=%",Et)

def fd3(x): return 24000*x**2-21600*x+4050

Ea=-(((b-a)**5)/(180*n**4))*(fd3(b)-fd3(a))
print("Ea=",Ea)