利用python画微分方程（组）的数值曲线

# This is a simple numerical example using matplotlib to simulate differential equatons, here we take y'= -sin(t), y(0)= 1. It's easy to compare with the standard plot method.
# Filename: sin.py

# Author: Qiang Xiao
# Time: 2015-06-19

import matplotlib.pyplot as plt    #matplotlib deal with plotting
import numpy as np    #numerical computation

t0= 0    #start point
te= 10    #end point
delta= 0.1    #step length

y= np.zeros(101)    #declare funtion value
dy= np.zeros(101)    #declare derivative value
y[0]= 1;    #initial value

for i in np.arange(100):
    dy[i]= -np.sin(i*delta);    #computing the first derivative
    y[i+1]= y[i]+ delta*dy[i];    #estimating function value by differential property

time= np.arange(t0,te+delta,delta);    

plt.