Page 1 sur 1

dsolve

Publié : ven. avr. 21, 2017 7:54 am
par alb
supposons(x>-pi/2 et x<pi/2)
dsolve(y''+y=1/cos(x))
renvoie
c_0*cos(x)+c_1*sin(x)+cos(x)*ln(2*cos(x))

Or une solution particuliere est
x*sin(x)+cos(x)*ln(cos(x))

Re: dsolve

Publié : ven. avr. 21, 2017 1:26 pm
par parisse
oui, ca vient d'une erreur dans l'integration de a:=-2/(2*i*exp(2*i*x)+2*i)*exp(2*i*x); b:=int(a); simplify(diff(b)-a);
il ne devrait pas y avoir de ln(abs()) quand on integre apres un changement de variables complexes.
En corrigeant ca, on a une solution qui verifie l'equation mais son expression est assez compliquee (presence de sign un peu partout).

Re: dsolve

Publié : ven. avr. 21, 2017 4:20 pm
par parisse
J'ai reussi a ameliorer la solution sous une hypothese plus restrictive, a tester sur Xcas pour Firefox pour le moment http://www-fourier.ujf-grenoble.fr/%7ep ... os(x))%3B&

Re: dsolve

Publié : ven. avr. 21, 2017 4:57 pm
par alb
l'algorithme de Xcas n'a rien à voir avec cette methode ?
https://fr.wikipedia.org/wiki/M%C3%A9th ... constantes

Re: dsolve

Publié : ven. avr. 21, 2017 5:26 pm
par parisse
Si, c'est exactement la methode utilisee. Mais les calculs sont faits avec des exponentielles complexes comme base de solutions.