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


a=0
b=0.8
n=2

Ireal=1.640533367

h=(b-a)/(2)
sum1=0
for iteration in range(1,n):
    sum1+=f(a+iteration*h) 
    
I=(b-a)*(f(a)+2*sum1+f(b))/(2*n)
print("I=",I)


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


def df(x): return 2000*x**4-3600*x**3+2025*x**2-400*x

ff=(df(b)-df(a))/(b-a)

print(ff)

Ea=-(ff*((b-a)**3))/(12*n**2)
print("Ea=",Ea)