singularité à l'équation différentielle avec MATLAB

Question

Je ne peux pas résoudre cette équation différentielle par ode45 parce qu'elle a la sigularité.

xy"=3xcos(x)+sin(x) ; x(0)=0 , x'(0)=0 

pouvez-vous m'aider à écrire ode45 fonction?

Answer 1

Vous pouvez utiliser la fonction sinc(x), qui est définie comme sin(π*x)/(π*x), sauf à x=0 où sa valeur est 1. Ainsi, vous pouvez réécrire votre ODE comme:

y'' = 3*cos(x) + sinc(x/π) 

qui ode45 ne devrait pas avoir de problème résoudre.