Analyse des algorithmes de résolution d'équations différentielles numériques

-&(M0B#>>(5"(B0&"(),&'-&(@),B,'#,&5<(B(''),&'0-<@,)'5*('1C

C"-+&A-&9())&$-$,$",$#&H9+6--))W&9$-&-&+&H96%A$9"&

>"7"$+

:(<@()<,&&(<?#1'#B(<"(@0<<#?#'=@0<<(&'>(-)<'(B@<Z*

'),-?()"(<M,->(<"()0#<,&<N-#>(<(B@S1/(&'"F0P#)50->#(-"(<(*

1,&1(&')()<-)>(<B,2(&<N-#@,-))0#(&'>(<0#"()ZB(&()ZU#(&

>(-)<*#"=(<C

:(<@()<,&&(<01'#?(<5(>>(<5@>0&#M#(&'>(-)<)S?(<"(M0b,&Z

>(<)=0>#<()1,&1)Q'(B(&'C!>>(<<(",&&(&'>(<B,2(&<"F2@0)?(&#)C

!>>(<@0)'(&'"-@)#&1#@(N-()#(&&(<(@),"-#)0A0B0#<<0&<-&(

)=(>>(#&#'#0'#?(@()<,&&(>>(C!>>(<@0<<(&'",&1Z>F01'#,&<0&<*

@()")("('(B@<C

:(<01'#M<0P#<<(&''0&"#<N-(>(<@0<<#M<<F0U<'#(&&(&'5>0#<<0&'0-

@0<<0P(>(<,@@,)'-&#'=<>(-)=1/0@@()5('>(-)1,&M#0&1(

@()<,&&(>>(<F(MM,&")()#&(J,)0U>(B(&'C

:(<01'#M<<0?(&'N-F#><@(-?(&'@)(&")(>(1,&')c>("(>(-)?#(5*

BSB("0&<>(<<#'-0'#,&<"#MM#1#>(<510)#><1,B@'(&'<-)>(-)@),@)(*

10@01#'=Z<()(>(?()>(<B0&1/(<('Z0P#)C

D

*

[3V7V*\\\C1>-U$("#'(-)$#&'()&('C1,B

*

W>>(Ga(:FW?0&'

!

:(<@0<<#M<,&'>F#B@)(<<#,&"(<-U#)>(-)(J#<'(&1(5',-'>(-)

<(BU>()(>(?()"(>0M0'0>#'=('"-U,&?,->,#)"(1(-JN-#>(<*

(&',-)(&'C

%>(<'=?#"(&'N-("0&<>(M,&"5>0@>-@0)'"(<@()<,&&(<*

@0<<#?(<0#B()0#(&')(A,#&")(>(P),-@("(<@()<,&&(<01'#?(<N-#5*

"(M0#'51,&&0#<<(&'-&(B(#>>(-)()=-<<#'("0&<>(-)?#(*

@),M(<<#,&&(>>(('@()<,&&(>>(C

.,BB(&'M0#)(0>,)<@,-)011="()Z1(P),-@(N-#1,&<')-#'

N-,'#"#(&&(B(&'<,&<-11Q<*e

Y

*

[3V7V*\\\C1>-U$("#'(-)$#&'()&('C1,B

*

W>>(Ga(:FW?0&'

!

<-6-8K&L&)2",,-,

Hh

-0&"AF0-)0#1/0&P="(')0?0#>5AF01/Q'()0#B0*

B0#<,&I5HN-0&"A(<()0#@),@)#='0#)(5A(

"(B0&"()0#>0B0#&"(B01,@#&(I5HN-0&"A(

<()0#B0)#=5A(?#<#'()0#>FW-<')0>#(I5*('1C

:0?#((<'-&(<-#'(@,'(&'#(>>("(HN-0&"I5('1F(<'U#(&<f)1(*

N-#(&M0#'',-'>(<(>_<#?,-<&F0?#(G@>-<0-1-&(M,)B("(

M)-<')0'#,&5?,-<&F0-)#(G(&?#("()#(&5('',-'?,-<<(BU>()0#'U#(&*

')#<'(*R

g(&()(B('<@0<(&10-<(>F#&'=)S'"(>F0''(&'(N-#@(-'@0)M,#<

S')(0-<<#"=>#1#(-<(N-(>F,UA(1'#M"(>F0''(&'(>-#$BSB(5('#>(<'*

@0)M0#'(B(&'>,P#N-(N-(1()'0#&(<01'#,&<"(?,')(?#(&=1(<<#'(&'*

-&@(-"(*@0'#(&1(C

:(?)0#@),U>QB("(>0@0<<#?#'=<F#&<'0>>(N-0&"?,')(?#(

(&'#Q)("(?#(&'-&(>,&P-(<-#'("F0''(&'(<<-11(<<#?(<5('N-(?,-<

@)(&(G>F/0U#'-"("F(&M0#)(?,')(M,&1'#,&&(B(&'P=&=)0>C

42",,-%K&H"+&)&X767-,&$%(")Y&8"9&$)&-2Z$+,&H"+

%>(<'U,&"(@0'#(&'()N-(>N-(<B,#<@,-)"(B0&"()-&(

@),B,'#,&50M#&"(B,&')()@(&"0&'1('(B@<$>Z1(N-F,&?0-'"0&<*

<,&')0?0#>5B0#<(<'$#>U,&"F0''(&")(@(&"0&'"(<0&&=(<5(&

(<@=)0&'N-(N-(>N-F-&<0-)05Z?,')(@>01(5'(&#)1,B@'("(1(N-(

?,-<"=<#)(G<(1)Q'(B(&'*e

;#?,-<?,->(G=@,-<()?,')(@0)'(&0#)(5(&<01/0&'

#&'=)#(-)(B(&'N-F#><F0P#'?)0#B(&'"(>0U,&&(@()<,&&(50?(G$?,-<

U(<,#&"F0''(&")(N-(',-<>(<M01'(-)<#"=0-J<,#(&'0U<,>-B(&'*

)=-&#<@,-)?,-<"=1>0)()50>,)<N-(?,-<>0#<<(G@(&"0&'1('(B@<$

>Z@0<<()-&'(B@<@)=1#(-J*e

:F0''(&'((<'Z>F,@@,<="(>F01'#,&5('#>M0-'@0)M,#<011(@'()"(

@)(&")(>(<1/,<(<1,BB((>>(<?#(&&(&'5(&)(M-<0&'"(?,-<

<,-B('')(Z>0&,'#,&"(HB,B(&'@0)M0#'I510)1(B,B(&'$>Z

&F(J#<'(A0B0#<?)0#B(&'C

:F0&'#1#@0'#,&(<'1()'(<-&(U,&&(1/,<(5B0#<,&&(B(<-)(

@0<>0)=-<<#'("(N-(>N-F-&<(>,&>010@01#'=N-F#>0Z)=P>()',-<>(<

K

*

[3V7V*\\\C1>-U$("#'(-)$#&'()&('C1,B

*

W>>(Ga(:FW?0&'

!

@),U>QB(<N-F#>@,-))0#'@,'(&'#(>>(B(&')(&1,&')()5@0)0?0&1(C

W>>()"(>F0?0&'(')=P>()>(<"#MM#1->'=<0-M#>"-'(B@<(<'>0*

B0)N-("(<@()<,&&(<01'#?(<N-#<,&'1,&<1#(&'(<N-F(>>(<?()),&'*

1()'0#&(B(&'"(<,U<'01>(<<-)>(-)),-'(5B0#<N-#&F(&M,&'@0<*

',-'(-&(/#<',#)(('<F(&<()?(&'BSB(0-1,&')0#)(1,BB("F-&(*

<,-)1("(B,'#?0'#,&C

O-&+"A,&%"-+&)&#$%&+,&H"9G6$+&)"&7$))A9&,8[-$PA

;#?,-<S'(<?,-<$BSB(/0U#'-=Z)(B('')((&@()B0&(&1(',-<

?,<@>0&<Z@>-<'0)"5&(@()<#<'(G@0<"0&<1('2@("(

1,B@,)'(B(&'10)>0)=-<<#'(&(M()0N-(?,-<M-#)C

W",@'(G",&1>(<0''#'-"(<<-#?0&'(<50M#&"()(&,&1()Z*

>F0''(&'(*_

• d-U>#(G>0&,'#,&"()#<N-(V_?,-<&(@,-))(GA0B0#<S')(<f)*

N-F-& @),A(' &( 1,B@,)'( 0-1-& )#<N-(5 ",&1 011(@'(G 1(*

@,<'->0'"Q<B0#&'(&0&'C

;#5@0)(J(B@>(5?,-<0?(G(&?#("(?,-<B0)#()5B0#<N-(*

?,-<0''(&"(GH>(U,&B,B(&'I51(>-#,i?,-<<()(G<f)"(&(

A0B0#<"#?,)1()5?,-<S'(<(&')0#&"(B('')(-&(1,&"#'#,&

#B@,<<#U>(Z>0)=0>#<0'#,&"(?,')(B0)#0P(*R

`,-<@,-?(G)="-#)(>(<M01'(-)<"()#<N-(5@0)(J(B@>(#1#5(&*

(J@>#N-0&'?,<1)0#&'(<Z?,')(@0)'(&0#)(5('(&1,&M),&'0&'*

?,<P)0&"<@),A('<"(?#(@,-)(&B(<-)()>F0"=N-0'#,&5B0#<

?,-<&(@,-))(GA0B0#<S')(<f)"(*)#(&C

L&<0-'"0&<>(?#"((<'@0)M,#<>0B(#>>(-)('(1/&#N-(@,-)&(

@0<)=M>=1/#)@(&"0&'"(<0&&=(<0-J)#<N-(<N-('(>>(,-'(>>(*

"=1#<#,&?,-<M()0#'*@)(&")(C

• +=P>(G?,<@),U>QB(<0-M-)('ZB(<-)(N-F#><<-)?#(&&(&'5

<0&<0''(&")(&,&@>-<N-F#><@)(&&(&'"(>F0B@>(-)5B0#<<0&<

>(<0&'#1#@()0-@,#&'"(&(@>-<@(&<()N-FZ*1(>0C

E(@0<<(G@0<',-'?,')('(B@<Z#B0P#&()N-('(>,-'(><,-1#*

@,-))0#'<-)?(&#)5B0#<?#?(G(&020&'1,&M#0&1((&?,-<$BSB(

('(&?,')(0@'#'-"(ZM0#)(M01(>(B,B(&'*?(&-C

O

*

[3V7V*\\\C1>-U$("#'(-)$#&'()&('C1,B

*

W>>(Ga(:FW?0&'

!

4&H9%K&H"+&%&,7H+W&"\$++K&%U+&7"$-,*-"-,&]

a(M0b,&)=<-B=(5>F0''(&'(1,BB(B,"("(?#((<'>F(&&(B#(

&0'-)(>>("(>0)=0>#<0'#,&"(<,#('@(-'?,-<(B@,#<,&&()

>F(J#<'(&1(5A-<N-FZ)=0>#<()5Z>F]P("(4V0&<5N-(?,-<0?(G0''(&"-

',-'(?,')(?#(@,-))#(&*R

;,2(G51()'(<5)0#<,&&0U>(@,-)&(@0<M0#)(',-'1(N-#?,-<*

@0<<((&@()B0&(&1(@0)>0'S'(('>()=0>#<()&F#B@,)'(1,BB(&'5*

B0#<011(@'(G5@0)1,&')(5"(H?,-<>]1/()I<-)1(<@),A('<N-#?,-<

'#(&&(&'Z1j-)"(@-#<<#>,&P'(B@<5<#?,-<&(?,->(G@0<M#&#)@0)*

?,#)?,')(?#(<("=),->()',-'(<(->("(?0&'?,<2(-J5<0&<N-(

?,-<<,2(G)=(>>(B(&'0-J1,BB0&"(<C

E(?,-<1,&'(&'(G@0<"F0''(&")(5B0#<0P#<<(G"Q<0-A,-)"F/-#C

4

*

[3V7V*\\\C1>-U$("#'(-)$#&'()&('C1,B

*

W>>(Ga(:FW?0&'

!

42-,99K&H"+&#6+&$%(+

g

(1,&&0#<@>-<#(-)<@()<,&&(<0-',-)"(B,#N-#,&')=-<<#

ZM0#)("#<@0)0k')((>>(<$BSB(<"(')Q<U,&&(<#"=(<N-F(>>(<*

0?0#(&'(-(<5(&)=-<<#<<0&'Z<(@()<-0"()N-F(>>(<&F(&*

?0>0#(&'@0<>0*@(#&(C

h-F#><F0P#<<("F=1)#)(-&>#?)(5"(B,&'()<0@),@)((&')(@)#<(5

"()=&,?()<0B0#<,&5('1C5?,-<@,-?(GS')(Z>0<,-)1("F-&(')Q<*

U,&&(#"=(0-@,'(&'#(>1()'0#&<#?,-<&(U0#<<(G@0<>(<U)0<0-*

"()&#()B,B(&'@0)1(N-(?,-<0?(G<,-"0#&@0<<='),@"('(B@<Z

)(P0)"()"(@)Q<>(<"='0#><"(<0*)=0>#<0'#,&C

;#?,-<?,-<1,&1(&')(G'),@>,&P-(B(&'<-)',-<>(<"='0#><

"F-&(,@=)0'#,&5#>(<'M,)'@),U0U>(N-(?,-<02(G(&?#("(U0#<<()>(<*

U)0<5('N-(?,-<0))#?#(G0>,)<Z)(B('')((&10-<(?,')(#"=(

"F,)#P#&(*(>>($BSB(C

a(>0BSB(M0b,&N-F,&1,&<(#>>(Z-&(&M0&'N-#(<'',BU="(*

?=>, "( )(B,&'() "(<<-< #BB="#0'(B(&' @,-) &( @0< <0&< 1(<<(

H)(?#?)(I <0 1/-'( B(&'0>(B(&'5 &,'(G N-( >F0''(&'( (' >F#&"=1#<#,&*

0>#B(&'(&'(&@()B0&(&1(>0@(-)"F=1/,-()('"(&(@0<)=-<<#)Z

011,B@>#)1(N-(>F,&0(&'S'(C

O-&$%(&+,&G"$,&H6A9&^,9&+A$#$&%2"8,$6-+

`,<#"=(<"(?)0#(&'',-A,-)<S')(<-#?#(<"F01'#,&<5@,-)&(@0<

',BU()"0&<>(",B0#&("-M0&'0<B(('"(>F#))=(>5@0)&0'-)(*

#&0''(#P&0U>(*_

• `,-<0?(G-&(#"=("(>#?)(eEF0''(&"(G@0<('1,&1(?(G*

"Q<0-A,-)"F/-#<,&@>0&(')="#P(G<,&#&'),"-1'#,&_>(M0#'"(

<0-'()>(@0<5(&020&'-&(01'#,&1,&1)Q'("0&<>01,&'#&-#'=*

"(?,')(#"=(5&(@,-))0N-(?,-<",&&()(&?#("(1,&'#&-()C*

;#5@0)1,&')(5?,-<0''(&"(G'),@>,&P'(B@<5>('/QB("-

>#?)(5@-#<>F#"=("->#?)((>>($BSB(('?,')(@),@)(10@01#'=Z*

>F=1)#)(<(),&'(&)0@#"(B(&')(B#<(&10-<(@0)?,')(*1()?(0-C

• `,-<0?(G(&?#("(1)=()?,')((&')(@)#<(eE(>0#<<(G@0<?,')(

#"=(?,-<=1/0@@()B0#<50-B,#&<5B(''(G(&@>01(*

T

*

[3V7V*\\\C1>-U$("#'(-)$#&'()&('C1,B

*

W>>(Ga(:FW?0&'

!

)0@#"(B(&'-&(='-"("(B0)1/=N-#?,-<",&&()0>0B(<-)(

1,&1)Q'("(>0M0#<0U#>#'="(?,')(#"=(C;#?,-<0''(&"(G'),@*

>,&P'(B@<5?,')(10@01#'=Z"(?(&#)?,')(@),@)(@0'),&?,-<*

<(BU>()0','0>(B(&'#&(J#<'0&'(5BSB(<#>(1,&1(@'N-(?,-<*

0?(G(&'S'((<'?)0#B(&'*?0>0U>(C

• `,-<0?(G(&?#("()=&,?()?,')(B0#<,&e9)(&(G)(&"(G$

?,-<"Q<B0#&'(&0&'0?(1-&@),M(<<#,&&(>N-#@,-))0?,-<*

P-#"()5,-01/('(G>(<B0'=)#0-J&=1(<<0#)(<@,-)?,-<>0&1()C;#

?,-<0''(&"(G'),@>,&P'(B@<5?,-<&(M()(GN-(?,-<",&&()

"(<)0#<,&<"FS')(M0'#P-=Z>F0?0&1(5('?,-<&(1/0&P()(G

A0B0#<)#(&Z?,')(#&'=)#(-)C

/+&$%(+&+6-,&H9(8$A++&,&HA#-,&^,9&+6A98+&%&

#(9$,";)+&9(A++$,+

`,<#"=(<<,&'@)=1#(-<(<('@(-?(&'S')(Z>0<,-)1("(*

?=)#'0U>(<)=-<<#'(<5@,-)@(-N-(?,-<<,2(G10@0U>("(M)0&1/#)1(*

',-'@('#'@0<N-#<=@0)(>F#&01'#,&"(*>F01'#,&C

l(''(G$?,-<"0&<>0@(0-"F-&@0)01/-'#<'(_>(<0-'(&>-#$*

BSB(@(-'S')(B()?(#>>(-J<F#>(<'011,B@>#B0#<&F0-)0A0B0#<>#(-*

<#?,-<)(<'(G"(</(-)(<(&>F0#)Z)(<<0<<()>0)0#<,&"FS')("-<0-'

('>(<@(-)<N-F#>?,-<*#&<@#)(C

E(1)0#P&(G@0<"(1,&M#()?,')(#"=(ZN-(>N-F-&N-#?,-<*

1,&&0k'U#(&('N-#&F/=<#'()0@0<Z?,-<",&&()>(@('#'1,-@"(@#("

0-JM(<<(<@0)M,#<&=1(<<0#)(@,-)(MM(1'-()1(<0-'"(*>0&1(B(&'C

E("(?(&(G@0<1(<@()<,&&(<N-(>F,&(&'(&"<#<,-?(&'"#)(*

H0/5<#<(->(B(&'AF0?0#<<-#?#B,&#"=(50-A,-)"F/-#B,&(&')(@)#<(

0-)0#'"=1,>>=I5,-H<#<(->(B(&'AF0?0#<<0-'=>(@0<5A(&F0-)0#<@0<

0-'0&'"()(P)('<ZB,&*]P(IC

=)"-$G$K&)&H)A+&9"H$%7-,&H6++$;)&%+&"8,$6-+&

86-89U,*+

L&(#"=((<'M0#'(@,-)S')(<-#?#("F01'#,&<('&,&@0<@,-)S')(*

(&'())=(0-M,&""F-&(U,k'(Z10-<("(<@(-)<N-F(>>(P=&Q)(C

7V

*

[3V7V*\\\C1>-U$("#'(-)$#&'()&('C1,B

*

W>>(Ga(:FW?0&'

!

;-#?(G",&1>(<?c')(<5@>0&#M#(G>(@>-<)0@#"(B(&'@,<<#U>(*

"(<01'#,&<1,&1)Q'(<5('1,&1(&')(G$?,-<-&#N-(B(&'<-)',-<>(<*

U=&=M#1(<N-#@,-))0#(&'"=1,->()"(1(''(#"=(C

:F01'#,&(<'>(<(->)(BQ"((MM#101(1,&')(>0@(-)50>,)<50P#<<(G5*

M0#'(<<F=?0&,-#)?,<@(-)<(')=1,>'(G>(<M)-#'<"(?,<#"=(<*R

77

*

[3V7V*\\\C1>-U$("#'(-)$#&'()&('C1,B

*

W>>(Ga(:FW?0&'

!

O,$)$+K&)2"8,$6-&7(8"-$PA

#A(?,-<"#<N-(@,-)<F2B('')(5#>M0-'<F2B('')(5?,-<*

0>>(G1()'0#&(B(&'@(&<()N-(A(?,-<<,)<-&(U(>>(*

=?#"(&1(R!'@,-)'0&'CCC

;

:(@>-<<,-?(&'5',-'(<>(<U(>>(<'/=,)#(<@<21/,>,P#N-(<&(*

?,-<0?0&1(),&'@0<@>-<N-(1(''(<#B@>(>,P#N-(N-#(<'N-(>F01'#,&

(&')0k&(*>F01'#,&C

E,')(B,&"((<'"(?(&-)(>0'#?(B(&'B,-Z>FQ)(B,"()&(5*

10)&,-<<,BB(<(&',-)=<"(B01/#&(<5"(P0"P('<('"(<2<'QB(<*

N-#M,&'>(')0?0#>Z&,')(@>01(Ca0&<-&(B0A,)#'="(10<51(<0#"(<

<,&'>(<U#(&?(&-(<10)(>>(<&,-<M,&'P0P&()"-'(B@<('&,-<

@()B(''(&'"(&,-<1,&1(&')()<-)>F(<<(&'#(>5B0#<(>>(<&,-<,&'5

"0&<>F(&<(BU>(5)(&"-&0'-)(>>(B(&'@0)(<<(-JC

D-$,$K&A-&"8,$6-&,&))&+&H6A9+A$#9"&%2))B7^7&"++K&

G"8$)7-,

E,-<0?,&<@()"->F/0U#'-"("(M0#)("(<(MM,)'<1,&1)('<('

)=@='=<@,-),U'(&#)1(N-(&,-<"=<#),&<5('@0)M,#<5&,-<&,-<*

<(&',&<"=@0<<=<Z>F#"=("("(?,#)&,-<B('')((&*B,-?(B(&'C

9,-)'0&'5#><-MM#'<,-?(&'"F#&#'#()-&(01'#,&@,-)N-F(>>(<(

@,-)<-#?("F(>>($BSB(0<<(GM01#>(B(&'5('',-'(>0"#MM=)(&1((&')(*

>F01'#M('>(@0<<#M?01,&<#<'()(&1(7()P(<'("(>0&1(B(&'5N-#*

0B,)1(>F01'#,&B=10&#N-(*_

• ;#?,-<N-#''(G>0'0U>(<0&<>0"=U0))0<<()('N-(?,-<@0)'(G

M0#)(0-')(1/,<((&@(&<0&'0-B,B(&'@=&#U>(,i?,-<

"(?)(G)(?(&#)@,-)&('',2()5?,-<&(M()(GN-()(@,-<<()*

>F=1/=0&1(('?,-<&(@),M#'()(G@0<?)0#B(&'"(1(N-(?,-<*

S'(<(&')0#&"(M0#)(CW>,)<>(?(G$?,-<(&@)(&0&'*

#BB="#0'(B(&'N-(>N-(<0<<#(''(<510)-&(M,#<>0B=10&#N-(

>0&1=(5?,-<M()(G&0'-)(>>(B(&'>(*)(<'(C

• ;#,&?,-<)(B('-&",<<#()Z')0#'()('N-(?,-<?,-<

1,&'(&'(G"(>F,U<()?()@,<=<-)?,')(U-)(0-5(&#B0P#&0&'Z

N-(>@,#&'1(>0?0S')("-)"(?,-<(&,11-@()5?,-<)(&"(G>(<*

73

*

[3V7V*\\\C1>-U$("#'(-)$#&'()&('C1,B

*

W>>(Ga(:FW?0&'

!

1/,<(<"(B,#&<(&B,#&<0U,)"0U>(<CW>,)<,-?)(G$>(<0&<*

)=M>=1/#)51,BB(&1(GZ>(>#)((&@)(&0&'"(<&,'(<('?,-<*(&

0-)(G)0@#"(B(&'@0)1,-)->0B,#'#=<0&<?,-<(&)(&")(*

1,B@'(C

/2"8,$6-&%6$,&+&H96%A$9&X7"$-,-"-,YW&,&-6-&H)A+&,"9%

:0M0b,&>0@>-<<#B@>("F-'#>#<()>F01'#,&B=10&#N-((<'"(

U0&&#)-&1()'0#&&,BU)("(B,'<"(?,')(?,10U->0#)(5>,)<N-(?,-<*

=?,N-(G>F01'#,&_H-&A,-)I5H"(B0#&I5H>0<(B0#&(@),1/0#&(I5('1C

<,&'0-'0&'"(B,'<N-#<,&'#&1,B@0'#U>(<0?(1>F01'#,&@-#<N-F#><

)(>Q?(&'"(>F#&'(&'#,&C

:F01'#,&<(@),"-#'HB0#&'(&0&'I5('-&#N-(B(&'B0#&'(&0&'5('

?,-<"(?(G0))S'...