import numpy as np
import scipy.integrate
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
fi0d = 5*0.0001
fi1d = 0.6
fi0g = 5*0.001
figd = 0.03
fi1g=3.3
fi2g = 6.5
fi3g = 2*(1.2**2)
gd = 3.31
md = 1.05
sigd = 0.96
gg = 9.1
mg = 10.5
sigg = 9.6

e = [ -0.2 -0.2,  -0.1  -0.02
def ode(h):
        return iscipy.ntegrate.quad(lambda h :(fi1d*mg + figd*md)*(((1/2)*sigg*fi1d+(1/2)*sigd*figd)*(e**2)+(gg*fi1d+gd*figd)*e+fi0g*fi1d+fi0d*figd))

def fun(h):
    return  figd*(fi0d/fi1d)-(fi1g-figd**2/fi1d)*h+fi2g*h**2-fi3g*h**3
t = np.linspace(-1, 1, 100)
plt.plot(ode(t) , fun(t))
plt.show()