title: Solving nonlinear differential first order equations using Python
body: <p>I would like to solve a nonlinear first order differential equation using Python. </p>
  
  <p>For instance, </p>
  
  <p>df/dt = f**4</p>
  
  <p>I wrote the following program, but I have an issue with matplotlib, so I don't know if the method I used with scipy is correct.</p>
  
  <pre><code>from scipy.integrate import odeint
  import numpy as np
  import matplotlib.pyplot as plt
  derivate=lambda f,t: f**4
  f0=10
  t=np.linspace(0,2,100)
  f_numeric=scipy.integrate.odeint(derivate,f0,t)
  print(f_numeric)
  plt.plot(t,f_numeric)
  plt.show()
  </code></pre>
  
  <p>Which results in the following error:</p>
  
  <pre><code>AttributeError: 'float' object has no attribute 'rint'
  </code></pre>