Question: how to solve first order ode numerically

Hi, I want to solve first order ode by using a numerical technique called shooting method, can anyone help me to write an algorithm for the
 

restart

sigma := x^2+1

x^2+1

(1)

M := .1

.1

(2)

N := .5

.5

(3)

``

diff(P(x), x) = -(3*(M-1))/(2*sigma^2)-3*N/sigma^3

diff(P(x), x) = 1.350000000/(x^2+1)^2-1.5/(x^2+1)^3

(4)

``

NULL

ics := P(-infinity) = 0

P(-infinity) = 0

(5)

``


 

Download Shooting_Help.mw
 

restart

sigma := x^2+1

x^2+1

(1)

M := .1

.1

(2)

N := .5

.5

(3)

``

diff(P(x), x) = -(3*(M-1))/(2*sigma^2)-3*N/sigma^3

diff(P(x), x) = 1.350000000/(x^2+1)^2-1.5/(x^2+1)^3

(4)

``

NULL

ics := P(-infinity) = 0

P(-infinity) = 0

(5)

``


 

Download Shooting_Help.mw

 

shooting method. please find attachment

Please Wait...