En fait, pour la solution générale d'une équation différentielle du type :
y'(x) + p(x)y(x) = q(x), tu fais ça.
Tu considères une fonction u.
u(x)y'(x) + u(x)p(x)y(x) = q(x)u(x).
En posant (u(x)y(x))' = u(x)y'(x) + u(x)p(x)y(x), tu as :
u(x)y'(x) + u(x)p(x)y(x) = u'(x)y(x) + y'(x)u(x)
Tu élimines, il te reste :
u'(x) = u(x)p(x)
=> u'(x)/u(x) = p(x)
=> ln|u(x)| = int(a,x,p(s)ds)
=> u(x) = exp(int(a,x,p(s)ds))
Et avant, tu as :
(u(x)y(x))' = q(x)u(x)
=> u(x)y(x) = int(b,x,q(s)*exp(-int(a,x,p(s)ds))ds)
=> y(x) =
exp(-int(a,x,p(s)ds)*int(b,x,q(s)*exp(-int(a,x,p(s
)ds))ds)