# -*- coding: utf-8 -*-

a=0
b=0.8

xd1=-0.5773
xd2=0.5773

xx1= ((b+a)+(b-a)*xd1)/2
xx2= ((b+a)+(b-a)*xd2)/2

dx= ((b-a)/2)


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


I=(f(xx1)*dx)+(f(xx2)*dx)
print('I:',I)

GaussInt=1.640533

Et=((GaussInt-I)/GaussInt)*100
print('Et:',Et)