Odeint is a modern C++ library for numerically solving Ordinary Differential Equations.

www.odeint.com/

ya Odeint C++ library only i used...with integrate method. But graphs are not matching with my matlab values

What method did you use in Matlab? How big is the difference?

ode15s in matlab

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.