daily pastebin goal
39%
SHARE
TWEET

Untitled

a guest Dec 16th, 2018 67 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. import numpy as np
  2. import scipy.integrate
  3. import matplotlib.pyplot as plt
  4. from mpl_toolkits.mplot3d import Axes3D
  5. fi0d = 5*0.0001
  6. fi1d = 0.6
  7. fi0g = 5*0.001
  8. figd = 0.03
  9. fi1g=3.3
  10. fi2g = 6.5
  11. fi3g = 2*(1.2**2)
  12. gd = 3.31
  13. md = 1.05
  14. sigd = 0.96
  15. gg = 9.1
  16. mg = 10.5
  17. sigg = 9.6
  18.  
  19. e = [ -0.2 -0.2,  -0.1  -0.02
  20. def ode(h):
  21.         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))
  22.  
  23. def fun(h):
  24.     return  figd*(fi0d/fi1d)-(fi1g-figd**2/fi1d)*h+fi2g*h**2-fi3g*h**3
  25. t = np.linspace(-1, 1, 100)
  26. plt.plot(ode(t) , fun(t))
  27. plt.show()
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top