0
$\begingroup$

I have a differential equation and when I try to solve it by using the "NDSolve" command, the Mathematica gives an error. The differential equation and error are as follows

NDSolve[{Y''[a] + 51/20 Y'[a]/a - 9/20 a^-3 Y[a] == 0, Y[1] == 1, 
Y'[1/1000] == 1}, Y, {a, 1/1000, 1}]

NDSolve: At a == 0.001`, step size is effectively zero; singularity or 
stiff system suspected.

I will be thankful if someone help.

Improve this question
$\endgroup$
1
  • $\begingroup$ Cannot reproduce the mentioned warning in v11.3 and v12. Which version are you in? $\endgroup$
    – xzczd ♦
    Commented Oct 4, 2019 at 6:21

2 Answers 2

Reset to default
2
$\begingroup$ 
DSolve[{y''[a] + 51/20 y'[a]/a - 9/20 a^-3 y[a] == 0, y[1] == 1, 
  y'[1/1000] == 1}, y, {a, 1/1000, 1}]

enter image description here

Improve this answer
$\endgroup$
2
$\begingroup$ 
NDSolve[{Y''[a] + 51/20 Y'[a]/a - 9/20 a^-3 Y[a] == 0, Y[1] == 1, 
Y'[1/1000] == 1}, Y, {a, 1, 1/1000}, WorkingPrecision -> 24]
Plot[Evaluate[Y[a] /. %], {a, 1, 1/1000}]

enter image description here

Improve this answer
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.