# Décomposition d'Adomian

La décomposition d'Adomian est une méthode semi-analytique de résolution d'équations différentielles développée par le mathématicien américain George Adomian (en) durant la seconde partie du XXe siècle. On rencontre fréquemment l'utilisation d'ADM pour Adomian Decomposition Method.

## Généralités

On considère le problème de Cauchy suivant :
${\dfrac {dy}{dt}}=f(t,y)\quad {\text{et}}\quad y(0)=y_{0}$
Cette équation vérifiée par y est générale dans la mesure où y peut être à valeurs vectorielles et que nous n'avons pas de condition sur f. Il faut bien noter que dans cette méthode, il est plus commode de considérer y comme un vecteur et f comme une fonction pour éviter des confusions :
${\begin{matrix}f:&\mathbb {R} \times \mathbb {R} ^{d}&\to &\mathbb {R} ^{d}\\&(t,y)&\mapsto &f(t,y)\end{matrix}}$
Considérant que f est analytique proche de y=y$_0$ et t=0, résoudre le problème initial revient à résoudre :
$y(t)=y_{0}+\int _{0}^{t}f\left(s,y(s)\right)\,\mathrm {d} s$

## Méthode d'Adomian

La méthode d'Adomian consiste à décomposer y comme une série :
$y=y_{0}+\sum _{n=1}^{\infty }y_{n}$
et à décomposer de la même manière la fonction f :
$f(t,y)=\sum _{n=0}^{\infty }A_{n}(t,y_{0},y_{1},...,y_{n}),$
où les fonctions A$_n$ sont les polynômes d'Adomian, calculés formellement comme-ci :
$A_{n}={\frac {1}{n!}}{\frac {\mathrm {d} ^{n}}{\mathrm {d} \varepsilon ^{n}}}f\left(t,\sum _{i=1}^{\infty }y_{i}\varepsilon ^{i}\right){\Bigg |}_{\varepsilon =0}.$
En injectant les deux premières décompositions dans l'équation intégrale, on peut en déduire une méthode itérative de calcul des y$_n$ :
$y_{n+1}=\int _{0}^{t}A_{n}\left(s,y_{0}(s)...,y_{n}(s)\right)\,\mathrm {d} s$
On part de la formulation intégrale
    $y(t)=y_{0}+\int _{0}^{t}f\left(s,y(s)\right)\,\mathrm {d} s\quad \Rightarrow \quad y(t)-y_{0}=\int _{0}^{t}f\left(s,y(s)\right)\,\mathrm {d} s$
Or la décomposition de y donne :
    $y(t)-y_{0}=\sum _{n=1}^{\infty }y_{n}(t)$
Soit par égalité :
    $\sum _{n=1}^{\infty }y_{n}(t)=\int _{0}^{t}f\left(s,y(s)\right)\,\mathrm {d} s$
Il suffit d'injecter la décomposition de la fonction f :
    $\sum _{n=1}^{\infty }y_{n}=\int _{0}^{t}\sum _{n=0}^{\infty }A_{n}(t,y_{0},y_{1},...,y_{n})\,\mathrm {d} s$
Dans un premier temps on ré-indice la série de gauche. Ensuite sous réserve que la convergence de la série des A$_n$ vers f est uniforme sur le segment [0,t] (ce qui est valide avec l'hypothèse d'analyticité de f), on intervertit somme et intégrale :
    $\sum _{n=0}^{\infty }y_{n+1}=\sum _{n=0}^{\infty }\int _{0}^{t}A_{n}(t,y_{0},y_{1},...,y_{n})\,\mathrm {d} s$
Puis on identifie les termes. Évidemment ce n'est pas l'unique solution mais elle nous satisfait entièrement. En effet, le calcul de y$_n$ en lui-même n'est pas très important, c'est la série qui découlera de ce calcul itératif qui nous intéressera (celle-ci construisant la solution y) :
$y_{n+1}=\int _{0}^{t}A_{n}(t,y_{0},y_{1},...,y_{n})\,\mathrm {d} s$
Pouvant calculer à la suite chaque y$_n$, on construit au fur et à mesure la solution finale y par sommation.

## Polynômes en dimension 1

Dans cette section on se restreint à d = 1. Les premiers polynômes d'Adomian (les dérivées correspondent à des dérivées partielles par rapport à y) sont alors :
    $A_{0}={\frac {1}{0!}}\ f\left(t,\sum _{i=1}^{\infty }y_{i}\varepsilon ^{i}\right){\Bigg |}_{\varepsilon =0}=f\left(t,y_{0}\right),$
    $A_{1}={\frac {1}{1!}}\ {\frac {\mathrm {d} }{\mathrm {d} \varepsilon }}f\left(t,\sum _{i=1}^{\infty }y_{i}\varepsilon ^{i}\right){\Bigg |}_{\varepsilon =0}=y_{1}{\frac {\partial }{\partial y}}f\left(t,y_{0}\right),$
    $A_{2}={\frac {1}{2!}}\ {\frac {\mathrm {d} ^{2}}{\mathrm {d} \varepsilon ^{2}}}f\left(t,\sum _{i=1}^{\infty }y_{i}\varepsilon ^{i}\right){\Bigg |}_{\varepsilon =0}=y_{2}{\frac {\partial }{\partial y}}f\left(t,y_{0}\right)+{1 \over 2}y_{1}{\frac {\partial ^{2}}{\partial y^{2}}}f\left(t,y_{0}\right).$