常微分方程解法（2）例题

用SciPyのscipy.integrate.odeint解联立微分方程式，以时间t（tの范围在0~几秒就可以，例：t=0~2秒）和解x(t), y(t), z(t)作图。
 
作为t=0的初期条件，x(t=0) = -10, y (t=0) = 0, z (t=0) = 35.0。还有，关于系数a,b,c，可以尝试a = 40, b = 5, c = 35 的情况和a = 40, b = 10, c = 35的情况。还有时间分割Δt要适当取小一点的值。另外，Δt如果取的过小，计算量会就变得很大。根据实际实行后的错误进行调整吧。

下面是解决的python程序：
import numpy as np
import scipy.integrate as sciin
import matplotlib.pyplot as plt

#把与自变量t相关的导函数的函数名字放在F里
def f(F, t,params):
x,y,z = F
f_values = [a*(y-x),(c-a)*x-x*z

常微分方程解法（2）例题最先出现在Python成神之路。