Files
FYS-STK4155/doc/Programs/LifeScience/diffeq.py
T
2018-04-06 12:01:43 -04:00

13 lines
232 B
Python

import numpy as np
t = np.linspace(0, 10, 21) # 20 intervals in [0, 10]
dt = t[1] - t[0]
N = np.zeros(t.size)
N[0] = 1
r = 0.5
for n in range(0, N.size-1, 1):
N[n+1] = N[n] + r*dt*N[n]
print 'N[%d]=%.1f' % (n+1, N[n+1])