Solving differential equations in c++ ?

Why C++? Can you test with Python?

You learn how to do Differential Equations in Python using the numpy and scipy packages on this website:

http://www.udacity.com/overview/Course/cs222/CourseRev/1

The course is Differential Equations in Action.

For example, you want to solve the following differential equation:

z''(t) = −z(t).

For a second order differential equation, you need to specify z(tini)

and z'(tini)! In the example, lets assume z(0) = 0.0005 and

z'(0) = 0.2.

Your code in python is:

#################################

import numpy as np

from scipy import integrate

import matplotlib.pyplot as plt

def deriv(y,t):

uprime = y[1]

wprime = - y[0]

yprime = np.array([uprime,wprime])

return yprime

start=0

end=10

numsteps=1000

time=np.linspace(start,end,numsteps)

y0=np.array([0.0005,0.2])

y=integrate.odeint(deriv,y0,time)

plt.plot(time,y[:])

plt.show()

##################################

Output is attached.