Equations différentielles

Avant de commencer, on pourra consulter la vidéo suivante présentant la méthode de variation de la constante.

1. On commence par résoudre l'équation homogène $y'+y=0$ dont la solution générale est $y(x)=\lambda e^{-x},$ $\lambda\in\mathbb R.$ On cherche une solution particulière sous la forme $y(x)=\lambda(x)e^{-x}$, de sorte que $y'(x)=\lambda'(x)e^{-x}-\lambda(x)e^{-x}$. On introduit ceci dans l'équation différentielle et on trouve $$\lambda'(x)e^{-x}-\lambda(x)e^{-x}+\lambda(x)e^{-x}=\frac{1}{1+e^x}.$$ Après simplification, ceci donne : $$\lambda'(x)e^{-x}=\frac{1}{1+e^x}\implies \lambda'(x)=\frac{e^x}{1+e^x}.$$ Une solution particulière est donc donné par $y(x)=\ln(1+e^x)e^{-x}$. Finalement, la solution générale de l'équation avec second membre est donnée par $$x\mapsto \ln(1+e^x)e^{-x}+\lambda e^{-x},\ \lambda\in\mathbb R.$$
2. On commence par résoudre l'équation homogène $(1+x)y'+y=0$, dont la solution générale est donnée par $y(x)=\frac{\lambda}{1+x}$, $\lambda\in\mathbb R$. On cherche une solution particulière par la méthode de variation de la constante, en posant $y(x)=\frac{\lambda(x)}{1+x}$, de sorte que $$y'(x)=\frac{\lambda'(x)}{1+x}-\frac{\lambda(x)}{(1+x)^2}.$$ En introduisant ceci dans l'équation différentielle, on trouve $$(1+x)\left(\frac{\lambda'(x)}{1+x}-\frac{\lambda(x)}{(1+x)^2}\right)+\frac{\lambda(x)}{1+x}=1+\ln(1+x)$$ ce qui donne après simplifications $$\lambda'(x)=1+\ln(1+x).$$ Une primitive est donnée par $\lambda(x)=(1+x)\ln(1+x)$, et la solution générale de l'équation avec second membre est donc donnée par $$x\mapsto\frac{\lambda}{1+x}+\ln(1+x),\ \lambda\in\mathbb R.$$
3. On commence par résoudre l'équation sans second membre $y'-\frac yx=0$. On remarque que $x\mapsto x$ est une solution. L'ensemble des solutions de l'équation sans second membre est donc les fonctions $x\mapsto Cx$, $C\in\mathbb R$.
   On cherche ensuite une solution particulière sous la forme $y(x)=\lambda(x) x$. Reportant dans l'équation différentielle, on trouve l'équation $\lambda'(x)=x$, ce qui donne $\lambda(x)=\frac{x^2}2+C$. L'ensemble des solutions de l'équation différentielle est donc donné par les fonctions $$x\mapsto Cx+\frac{x^3}{2},\ C\in\mathbb R.$$
4. On commence par résoudre l'équation homogène $y'-2xy=0$. Pour chercher une solution, on peut remarquer que si $y$ est une solution qui ne s'annule pas, $$y'-2xy=0\iff \frac{y'}{y}=2x\iff \ln|y|=x^2+C.$$ Ainsi, ceci nous conduit à observer que $x\mapsto e^{x^2}$ est solution de l'équation homogère et donc que la solution générale de l'équation homogène est $x\mapsto \lambda e^{x^2},$ $\lambda\in\mathbb R.$ On cherche ensuite une solution particulière de l'équation en utilisant la méthode de variation de la constante. On pose donc $y(x)=\lambda (x) e^{x^2}$ et introduisant $y$ dans l'équation avec second membre, on trouve $$\lambda'(x)e^{x^2}=(-2x+1)e^x\iff \lambda'(x)=(-2x+1)e^{-x^2+x}.$$ Une primitive est donnée par $\lambda(x)=e^{-x^2+x}$ et donc une solution particulière de l'équation avec second membre est donnée par $$x\mapsto e^{-x^2+x}e^{x^2}=e^{x}.$$ Finalement, les solutions de l'équation sont les fonctions $x\mapsto \lambda e^{x^2}+e^x,$ $\lambda\in\mathbb R.$
5. On commence par résoudre l'équation homogène $y'-\frac2t y=0$. On trouve que les solutions sont les fonctions de la forme $y(t)=\lambda t^2,$ $\lambda\in\mathbb R.$ On cherche une solution particulière par la méthode de variation de la constante en posant $y(t)=\lambda(t)t^2$. L'équation devient : $$t^2=y'(t)-\frac2t y(t)=\lambda'(t) t^2.$$ Dès lors, $\lambda'(t)=1$ soit $\lambda(t)=t+C$. Finalement, les solutions sur $]0,+\infty[$ de l'équation de départ sont les fonctions $$t\mapsto t^3+Ct^2,\ C\in\mathbb R.$$

On va procéder par analyse/synthèse. Pour cela, on écrit qu'une solution $f$ s'écrit $$f(x)= \frac{C+x}{1+x^2}\iff C=(1+x^2)f(x)-x.$$ Le but étant d'éliminer $C$, on dérive la relation précédente, et on trouve $$0=(1+x^2)f'(x)+2xf(x)-1.$$ C'est l'équation différentielle recherchée! En effet, réciproquement, si l'on considère l'équation différentielle $$(1+x^2)y'(x)+2xy(x)-1=0,$$ on prouve facilement que ses solutions sont les fonctions $ x\mapsto \frac{C+x}{1+x^2}.$

1. 1. Il est clair que $\lim_{x\to 0^+}f(x)=0$, indépendamment de la valeur de $C$, et que $\lim_{x\to 0^-}f(x)=\pm\infty$ si $D\neq 0$, et $\lim_{x\to 0^-}f(x)=0$ si $D=0$. Ainsi, on a un prolongement continu en $0$ si et seulement si $D=0$. Dans ce cas, on a $f(0)=0$.
   2. On suppose donc que $D=0$. La fonction $f$ étant identiquement nulle à gauche de 0, elle est dérivable à gauche en 0 et sa dérivée est nulle. Pour $x>0$, on a $$\frac{f(x)-f(0)}x=\frac{C}x\exp(-1/x).$$ Posons $u=\frac 1x$. Lorsque $x$ tend vers $0^+$, $u$ tend vers $+\infty$ et $$\frac{1}x\exp(-1/x)=u\exp(-u).$$ Par comparaison des fonctions polynomes et exponentielle, on en déduit que $\frac{f(x)-f(0)}x$ tend vers 0 lorsque $x$ tend vers $0^+$, et donc $f'$ est dérivable à droite en 0, de dérivée nulle. Ainsi, on a bien que $f$ est dérivable en 0, avec $f'(0)=0$. La continuité à gauche de $f'$ en 0 ne pose alors pas de problèmes. En ce qui concerne la dérivée à droite, on remarque que, pour $x>0$, on a $$f'(x)=\frac C{x^2}\exp(-1/x)=Cu^2\exp(-u)$$ toujours avec le même changement de variables. Comme précédemment, on en tire que $\lim_{x\to 0^+}f'(x)=0$, et donc que $f'$ est continue en 0.
2. Sur l'intervalle $]0,+\infty[$, la fonction $x^2$ ne s'annule pas et l'équation est équivalente à $$y'=\frac1{x^2}y.$$ Les solutions de cette équation sont les fonctions $y(x)=C\exp(-1/x)$, où $C\in \mathbb R$. La résolution sur l'intervalle $]-\infty,0[$ donne exactement le même ensemble de solutions.
3. Soit $y$ une solution sur $\mathbb R$. Sa restriction à $]0,+\infty[$ est solution sur $]0,+\infty[$, et donc il existe une constante $C\in\mathbb R$ telle que, pour tout $x>0$, $y(x)=C\exp(-1/x)$. La restriction de $y$ à $]-\infty,0[$ est aussi solution sur $]-\infty,0[$, et donc il existe une constante $D\in\mathbb R$ telle que, pour tout $x<0$, $y(x)=D\exp(-1/x)$. Remarquons ici que $C$ et $D$ n'ont aucune raison d'être égaux. En effet, les résolutions sur $]-\infty,0[$ et $]0,+\infty[$ se font totalement indépendamment.
   D'ailleurs, le résultat des premières questions entraîne que, pour que $y$ soit continue en $0$, il est nécessaire que $D=0$. Dans ce cas, la fonction $y$ est de classe $C^1$, et elle vérifie bien l'équation différentielle : c'est clair pour $x\neq 0$, et c'est aussi vrai en 0 par continuité de $y$ et $y'$ en 0.

On va chercher une équation différentielle linéaire du second ordre admettant exactement cette famille de solutions. Pour cela, il suffit de trouver une équation différentielle dont l'équation caractéristique admet $2$ et $-1$ comme racine, c'est-à-dire que l'on souhaite qu'elle se factorise en $(r+1)(r-2)=r^2-r-2.$ Une équation différentielle qui convient est donc $$y''-y'-2y=0.$$

D'abord, l'équation $z''=0$ donne facilement $z(t)=at+b$, avec $a$ et $b$ des constantes. Ensuite, on a $$u'=x''+iy''=\omega y'-i\omega x'=-i\omega u.$$ On en déduit que $u(t)=(c+id) e^{-i\omega t}$. En prenant les parties réelles et imaginaires, on trouve que \begin{eqnarray*} x'(t)&=&c\cos(\omega t)+d\sin(\omega t)\\ y'(t)&=&d\cos(\omega t)-c\sin(\omega t). \end{eqnarray*} Il suffit d'intégrer une nouvelle fois pour trouver les valeurs de $x$ et de $y$ : \begin{eqnarray*} x(t)&=&c'\sin(\omega t)-d'\cos(\omega t)+x_0\\ y(t)&=&d'\sin(\omega t)+c'\cos(\omega t)+y_0 \end{eqnarray*} où on a posé $c'=c/\omega$ et $d'=d/\omega$. Il s'agit d'un mouvement hélicoïdal.

L'accroissement de la population est mesuré par $P'(t)$, cette fonction est donc proportionnelle à $P(t)$. Autrement dit, $P$ vérifie une équation différentielle du type $P'(t)=kP(t)$. Sa solution est de la forme $P(t)=Ce^{kt}$. De l'autre information (la population double tous les 50 ans), on déduit que $P(t+50)=2P(t)$ pour tout $t\geq 0$. Ainsi, on a $$Ce^{kt}e^{50k}=2Ce^{kt}\implies 50k=\ln 2\implies k=\ln2/50.$$ On cherche $x$ tel que $P(t+x)=3P(t)$ pour tout $t\geq 0$. Ceci donne $$Ce^{kt}e^{kx}=3Ce^{kt}\implies kx=\ln 3\implies x=\frac{50\ln 3}{\ln 2}\simeq 79,24.$$ La population triple environ tous les 79 ans et un quart.

On note $x(t)$ la quantité restante en fonction du temps, de sorte que $x(0)=20$. Du fait que la vitesse de dissolution est proportionnelle à la quantité restante, on tire qu'il existe $\alpha>0$ de sorte que $x$ est solution de l'équation différentielle $$-x'(t)=\alpha x(t).$$ La résolution de cette équation différentielle donne $x(t)=Ce^{-\alpha t}$, avec $C\in\mathbb R$. De $x(0)=20$ on tire $C=20$. Puis, de $x(5)=10$, on tire $$20e^{-5\alpha}=10\iff \alpha =\frac15\ln 2.$$ La quantité restante vaut donc $x(t)=20\exp(-(\ln 2)t/5)$. On cherche enfin $t_0$ de sorte que $x(t_0)=1$. Il vient $$20\exp(-(\ln 2)t_0/5)=1\iff -(\ln 2)t_0/5=-\ln 20\iff t_0=\frac{5\ln 20}{\ln 2}.$$ Le temps d'attente restant vaut donc $t_0-5$ minutes. Une application numérique donne environ 16 minutes et 36 secondes.

Soit $M=(t,f(t))$ un point de $C_f$. Le point $P$ a pour coordonnées $(t,0)$. La tangente à $C_f$ en $M$ a pour équation $$y-f(t)=f'(t)(x-t).$$ L'intersection de la tangente avec l'axe $(O,\vec i)$ a pour coordonnées $(a,0)$ où $a$ vérifie $$-f(t)=f'(t)(a-t)\iff a=t-\frac{f(t)}{f'(t)}.$$ Le vecteur $\overrightarrow{TP}$ a donc pour coordonnées $\left(\frac{f(t)}{f'(t)},0\right)$. On cherche les fonctions $f$ pour lesquelles ce vecteur est constant, c'est-à-dire pour lesquelles il existe $k\in\mathbb R$ de sorte que $\frac{f(t)}{f'(t)}=k$. $k$ ne peut pas être nul, sinon la fonction $f$ serait identiquement nulle et sa dérivée aussi. Posons $\lambda=1/k$. Alors on cherche les fonctions $f$ solutions de l'équation différentielle $f'-\lambda f=0$. Ces fonctions sont de la forme $f(t)=Ce^{\lambda t}$ avec $C\neq 0$ pour la même raison que ci-dessus. On a donc prouvé que les fonctions solutions sont à rechercher parmi les fonctions $f(t)=Ce^{\lambda t}$ avec $C,\lambda\in\mathbb R^*$.
Réciproquement, on vérifie facilement que ces fonctions sont solutions (il ne faut pas oublier de traiter la réciproque, ou alors il faut faire très attention de raisonner plus haut par équivalence).

Faisant $s=t=0$, on remarque que $f(0)^2=f(0)$, et donc $f(0)=0$ ou $f(0)=1$. Si $f(0)=0$, faisant $s=0$, on trouve que pour tout $t\in\mathbb R$, on a $$f(t)=f(0)f(t)=0.$$ On peut donc supposer $f(0)=1$. On repart alors de la relation initiale, on fixe la variable $s$ à une valeur quelconque de $\mathbb R$, et on dérive par rapport à $t$. On obtient $$f'(s+t)=f(s)f'(t).$$ On évalue en $t=0$, et on trouve $$f'(s)=f(s)f'(0).$$ Ainsi, $f$ est solution d'une équation différentielle $y'=\alpha y$. On en déduit que $f$ est de la forme $f(x)=Ce^{\alpha x}.$ Pour que $f(0)=1$, il est nécessaire que $C=1$. Réciproquement, on vérifie aisément que la fonction nulle et les fonctions $x\mapsto e^{\alpha x}$, $\alpha\in\mathbb R$, sont bien solutions de l'équation fonctionnelle.

On pose $g=f+f'$. Alors $f$ est solution de l'équation différentielle $f+f'=g$. On résoud cette équation. L'équation homogène est $f'+f=0$ dont la solution générale est donnée par $\lambda e^{-x}$. On résout l'équation avec second membre par la méthode de variation de la constante : en posant $f(x)=\lambda(x)e^{-x}$, on trouve $$\lambda'(x)e^{-x}=g(x),$$ et une solution particulière est donnée par $$f_0(x)=e^{-x}\int_0^{x}g(t)e^tdt.$$ Finalement, toute fonction $f$ vérifiant $f+f'=g$ s'écrit $$f(x)=\lambda e^{-x}+e^{-x}\int_0^{x}g(t)e^tdt.$$ Pour montrer que $f$ tend vers 0 en $+\infty$, il suffit de prouver que $e^{-x}\int_0^{x}g(t)e^tdt$ tend vers 0 lorsque $x$ tend vers $+\infty$. Pour cela, on va utiliser que $g$ tend vers 0 en $+\infty$, et on va couper l'intégrale en 2. Voici l'idée. Soit $A>0$ arbitraire pour le moment. Alors, pour $x\geq A$, on a $$ \left|e^{-x}\int_0^x g(t)e^{t}dt\right|\leq e^{-x}\int_0^A |g(t)e^{t}|dt+e^{-x}\int_A^x |g(t)|e^{t}dt.$$

• Si $A$ est fixé, alors on peut rendre le premier terme petit en choisissant $x$ suffisamment grand.
• Sur $[A;x]$, on peut rendre $|g(t)|$ petit à condition d'avoir choisi $A$ assez grand, ce qui nous permettra de rendre le deuxième terme petit.

Reste à effectuer cette démarche dans le bon ordre. Soit $\varepsilon>0$. Il existe $A>0$ tel que, pour $t>A$, on a $|g(t)|\leq \varepsilon$. Soit $M=\int_0^A |g(t)e^{t}|dt$ et soit $B\geq A$ tel que, pour $x\geq B$, on a $e^{-x}M\leq\varepsilon$. Alors, pour tout $x\geq B$, il vient \begin{eqnarray*} \left|e^{-x}\int_0^x g(t)e^{t}dt\right|&\leq& e^{-x}\int_0^A |g(t)e^{t}|dt+e^{-x}\int_A^x |g(t)|e^{t}dt\\ &\leq&e^{-x}M+e^{-x}\int_A^x \varepsilon e^{t}dt\\ &\leq&2\varepsilon. \end{eqnarray*} On a bien prouvé que $\lim_{+\infty}f=0$.

On va procéder par analyse synthèse. Commençons par considérer $f$ une fonction solution de l'équation. Par récurrence, $f$ de classe $C^\infty$. Il est donc légitime de dériver cette équation, et on trouve $$f''(x)=-f'(\lambda-x)=-f\big(\lambda-(\lambda-x)\big)=-f(x).$$ Ainsi, $f$ est solution de l'équation différentielle $f''+f=0$. Il existe donc des constantes $C$ et $\theta$ telles que $$f(x)=C\cos(x+\theta).$$ On va chercher maintenant des conditions sur $C$ et $\theta$ pour que $f$ soit effectivement solution de l'équation. De $f'(x)=f(\lambda-x)$, on trouve, pour $C\neq 0$, $$-\sin(x+\theta)=\cos(\lambda-x+\theta)$$ soit $$\cos(\pi/2+\theta+x)=\cos(\lambda-x+\theta).$$ Ceci entraîne, pour un $x\in\mathbb R$ fixé, $$\frac\pi2+\theta+x=\lambda-x+\theta\ [2\pi]\textrm{ ou } \frac\pi2+\theta+x=-\lambda+x-\theta\ [2\pi].$$ Maintenant l'égalité $$\frac\pi2+\theta+x=\lambda-x+\theta\ [2\pi]$$ entraîne $$x=\frac\lambda 2-\frac \pi4\ [\pi].$$ Elle ne peut donc pas être satisfaite pour toutes les valeurs de $x$, et donc on sait qu'il existe des réels $x$ pour lesquels l'égalité $$\frac\pi2+\theta+x=-\lambda+x-\theta\ [2\pi]$$ est vraie. Ceci implique qu'il existe $k\in\mathbb Z$ tel que $$\theta=-\frac{\lambda}{2}-\frac{\pi}{4}+k\pi.$$ Ainsi, si $f$ est solution, il existe $C\in\mathbb R$ et $k\in\mathbb Z$ tels que $$f(x)=C\cos\left(x-\frac{\lambda}{2}-\frac{\pi}{4}+k\pi\right).$$ Tenant compte de la périodicité de la fonction cosinus et du fait que $\cos(u+\pi)=-\cos(u)$, on peut simplifier ceci en disant que si $f$ est solution, alors il existe $C\in\mathbb R$ tel que $$f(x)=C\cos\left(x-\frac{\lambda}{2}-\frac{\pi}{4}\right).$$ Réciproquement (c'est la synthèse!), on vérifie facilement que les fonctions définies par la formule précédente satisfont l'équation.

Soit $f$ une solution de cette équation. Alors, on a pour tout réel $x$, $f'(x)=e^x-f(-x)$. Ainsi, la fonction $f'$ est de classe $C^1$ et donc $f$ est de classe $C^2$. Dérivant cette équation, on trouve $$f''(x)-f'(-x)=e^x.$$ Or, $f'(-x)=e^{-x}-f(x),$ et donc l'équation devient $$f''(x)+f(x)=e^x+e^{-x}.$$ On résoud alors cette équation différentielle linéaire du second ordre. Les solutions de l'équation sans second membre sont les fonctions $x\mapsto \lambda\cos(x)+\mu\sin(x)$, $\lambda,\mu\in\mathbb R$. Une solution particulière est donnée par la fonction $\cosh$. Ainsi, si $f$ est solution de l'équation initiale, on sait qu'il existe $\lambda,\mu\in\mathbb R$ tels que, pour tout $x\in\mathbb R$, on a $$f(x)=\lambda\cos x+\mu\sin x+\cosh(x).$$ Ici, il faut prendre garde au fait que l'on n'a pas procédé par équivalence. En effet, on a dérivé l'équation, et rien ne dit que toute solution de l'équation différentielle du second ordre est une solution de l'équation initiale. On doit vérifier pour quelles valeurs de $\lambda$ et $\mu$ c'est effectivement le cas. On a : $$f'(x)=-\lambda\sin x+\mu\cos x+\sinh x$$ et on doit donc avoir, pour tout réel $x$ $$(-\lambda-\mu)\sin x+(\mu+\lambda)\cos x+e^x=e^x.$$ Si on fait $x=0$, on trouve que $\mu=-\lambda$. Et réciproquement, toute fonction $x\mapsto \lambda(\cos x-\sin x)+\cosh x$, $\lambda\in\mathbb R$, est solution de l'équation.

Soit $y$ une solution de l'équation. Sa tangente au point d'abscisse $x_0$ a pour équation $$y-y(x_0)=\big(a(x_0)y(x_0)+b(x_0)\big)(x-x_0).$$ Par tout point $(x_0,\lambda)$ il passe une courbe intégrale (ou encore il y a une unique solution avec $y(x_0)=\lambda$), et il s'agit de démontrer que toutes les droites de la famille $(D_\lambda)_{\lambda\in\mathbb R}$, où $(D_\lambda)$ est la droite d'équation $$y-\lambda=\big(a(x_0)\lambda+b(x_0)\big)(x-x_0)$$ sont parallèles ou concourantes. Si $a(x_0)=0$, il est clair qu'elles sont parallèles, et parallèles à $y=b(x_0)x$. Sinon, cherchons le point d'intersection, pour $\lambda\neq \mu$, de $D_\lambda$ et $D_\mu$. On doit résoudre le système $$\left\{ \begin{array}{rcl} y-\lambda&=&\big(a(x_0)\lambda+b(x_0)\big)(x-x_0)\\ y-\mu&=&\big(a(x_0)\mu+b(x_0)\big)(x-x_0)\\ \end{array}\right.$$ On trouve que le point d'intersection a pour coordonnées $\left(x_0-\frac{1}{a(x_0)},\frac{-b(x_0)}{a(x_0)}\right)$. Il ne dépend pas de $\lambda$ ni de $\mu$. Toutes les droites $D_\lambda$ passent par ce point!

La solution générale de l'équation s'écrit $$y(x)=\left(\alpha+\int_0^x b(t)e^{-A(t)}dt\right)e^{A(x)}$$ où $A(x)=\int_0^x a(t)dt$. Notons que $$A(x+1)=A(x)+\int_x^{x+1}a(t)dt=A(x)+\int_0^1 a(t)dt,$$ et posons $\lambda=\int_0^1 a(t)dt$. On a ainsi : \begin{eqnarray*} y(x+1)&=&\left(\alpha+\int_0^1b(t)e^{-A(t)}dt+\int_1^{x+1}b(t)e^{-A(t)}dt\right)e^{A(x+1)}\\ &=&\left(\alpha+\int_0^1b(t)e^{-A(t)}dt+\int_0^{x}b(t)e^{-\lambda-A(t)}dt\right)e^{A(x)+\lambda}\\ &=&y(x)+\left(\alpha(e^\lambda-1)+\mu e^{\lambda}\right)e^{A(x)} \end{eqnarray*} où on a posé $\mu=\int_0^1 b(t)e^{-A(t)}dt$. Autrement dit, $f$ est $1-$périodique si et seulement si $$\alpha(e^\lambda-1)+\mu e^{\lambda}=0.$$ Si $\lambda\neq 0$, l'équation admet une unique solution $1-$périodique, donnée par $$\alpha=\frac{\mu e^{\lambda}}{1-e^{\lambda}}.$$ Si $\lambda=0$ et $\mu=0$, alors toute solution est 1-périodique.
Si $\lambda=0$ et $\mu\neq 0$, alors il n'y a aucune solution 1-périodique.

Il y a deux clés pour résoudre cet exercice :

• toute fonction impaire vaut $0$ en $0$;
• l'équation différentielle $(E)$ admet une unique solution $y_0$ vérifiant $y_0(0)=0$.

Ceci montre déjà l'unicité : s'il y a une fonction $y$ impaire solution de $(E)$, elle vérifie $y(0)=0$ et doit donc être égale à $y_0$. Réciproquement, on doit prouver que $y_0$ est impaire. On va poser $z(t)=-y_0(-t)$. $z$ est solution de $(E)$. En effet, $$z'(t)+a(t)z(t)=y_0'(-t)-a(t)y_0(-t)=y_0'(-t)+a(-t)y_0(-t)=b(-t)=b(t),$$ car $y_0$ est solution de $(E)$, $a$ est impaire et $b$ est paire. $z$ est donc solution de $(E)$, et satisfait de plus $z(0)=0$. Ainsi, par unicité au problème de Cauchy, $z$ est égale à $y_0$, et donc $y_0$ est impaire.
On pouvait aussi prouver que $y_0$ est impaire, en cherchant à résoudre l'équation différentielle par la méthode usuelle (solution de l'équation homogène à l'aide de l'exponentielle et d'une primitive de $a$, puis méthode de variation de la constante).

L'équation caractéristique $r^2+ar+b=0$ a pour discriminant $\Delta=a^2-4b$. On distingue alors les trois cas :

• Si $\Delta>0$, les solutions sont de la forme $x\mapsto C_1\exp(\lambda_1 x)+C_2\exp(\lambda_2 x)$, avec $\lambda_1,\lambda_2$ des réels distincts. L'un des deux, disons $\lambda_1$, est non-nul. La solution $x\mapsto \exp(\lambda_1 x)$ n'est jamais bornée.
• Si $\Delta=0$, alors les solutions sont de la forme $x\mapsto (C_1 x+C_2)\exp(\lambda x)$ avec $\lambda$ un réel. Dans tous les cas (même si $\lambda=0$), la solution $x\mapsto x\exp(\lambda x)$ n'est pas bornée.
• Si $\Delta<0$, alors posons $\delta>0$ tel que $\delta^2=-\Delta$. Les racines de l'équation caractéristique sont $\frac{-a\pm i\delta}2$ et les solutions de l'équation sont les fonctions de la forme $x\mapsto C_1\sin(\delta x/2)\exp(-ax/2)+C_2\cos(\delta x/2)\exp(-ax/2)$. Or, les fonctions $x\mapsto \sin(\delta x/2)\exp(-ax/2)$ et $x\mapsto \cos(\delta x/2)\exp(-ax/2)$ sont bornées si et seulement si $a=0$.

Ainsi, en résumant, toutes les solutions de $y''+ay'+by=0$ sont bornées si et seulement si $\Delta<0$ et $a=0$, c'est-à-dire si et seulement si $a=0$ et $b>0$. En physique, cela correspond à l'équation différentielle obtenue pour un oscillateur harmonique.