I am having a problem coding a Matlab code that solves a second ODE equation which I give below:

x^3*(1-2*x*M)d^2J(x)/dx^2+2*(2*x^2+i*nu*x-7*x^3*M)*dJ(x)/dx

-2*(2*x+8*M*x^2+i*nu)*J(x)=0.

where

M = 1 (Mass of a black hole),

nu = 0.74734+0.17792*i,

J is a function of x,

i represents a complex number.

This equation is very complex to solve. Fistly I tried the Matlab shooting

method and did not work because this equation is singular at the boundries i.e

0 and 1. Now I am trying to solve it using Matlab Finite Difference Method. I must point out that

this problem has no analytic solutions and we hope that it can be solved numerically.

Can any one out there help me with this problem.

Thanks

Amos

x^3*(1-2*x*M)d^2J(x)/dx^2+2*(2*x^2+i*nu*x-7*x^3*M)*dJ(x)/dx

-2*(2*x+8*M*x^2+i*nu)*J(x)=0.

where

M = 1 (Mass of a black hole),

nu = 0.74734+0.17792*i,

J is a function of x,

i represents a complex number.

This equation is very complex to solve. Fistly I tried the Matlab shooting

method and did not work because this equation is singular at the boundries i.e

0 and 1. Now I am trying to solve it using Matlab Finite Difference Method. I must point out that

this problem has no analytic solutions and we hope that it can be solved numerically.

Can any one out there help me with this problem.

Thanks

Amos

Last edited: