Simulating continuous dynamical systems under conditions of ... - ULB

0 downloads 0 Views 894KB Size Report
ous dynamical system and we compare it with consolidated results in probability theory. The analysis illustrates the limits of conventional Monte Carlo ...
       "!#$ !&%'( ) *+,!&%"  -/.' 0! 1 %)2 2( -3 4 )0! 1 5""26 74 11 %)5 890!: ; =@?BADCFE=HGJI@?LKNMOPQ< RHT(S UWV XZY[T(\8T(]^Y`_Qa bc]edgfX[\VY[hji^kLTml@no]LhqprTsX[tZhuY T&S vFhqw@X[T_eT&xyX[k@z@T({q{jTst(l|x"T({q}hjke\ Ts\Vhq{~Z}mw|f]€Y[T6JkL{qwF‚ Vƒ‚ w|T „…‡†(ˆ(‰ ŠL‹mˆmŒ  X[fwWVwehq{jhuYcŽ'V]L_U|fmtZtZhqwLhq{qhuYŽY[LTsfXZŽ_@T6V{y‘`hqY[’_@hq“BTNX[T(]^Y:YcŽ@U|T(t#fdJke]Lƒ(TNXZ” Y•Vhq]€YcŽ^‚^– eToULVU|TsX"_ehqtZƒ(ketZtZT(t"Y[LT)X[f{jTofdBY[LTstZT9d>fX[\V{qhqtZ\8t"hq]8Y[LT`tZhq\ ke{jVY[hqfm]fdQƒ(f]^Y[hj]@” kefmketo_€Ž@]WV\8hqƒ6V{Bt—Žet—Y[Ts\8t ‘`LTsX[TULVX•V\8TNY[TsX[t`V]L_W˜fX9hj]ehuY[hDV{7ƒ(f]W_@hqY[hqf]Lt9VX[T&kL]eƒ(TsXZY•Vhj]Q‚ ™ T_eTNX[hjpmTV]LV{uŽ€Y[hjƒ(V{q{uŽY[LTJT(prfm{qk@Y[hjf]{jV(‘šfd‡U|fmtZtZhqwehj{qhuYŽH_ehqt—YZX[hqwLk@Y[hjf]Lt)hq]Hƒ(f]^Y[hj]€k@” fkLt&_€Že]LV\8hqƒ6V{-t—Ž@t—Y[T(\›V]L_‘"Tœƒ(f\8UWVX[T#huY‘`huY[žƒ(f]LtZfm{qhj_LV Y[T6_ X[TstZkL{uY[thq]žUeX[fwWVwLhq{qhqYcŽ Y[eT(fXZŽ^‚|– LTœV]WV{qŽ@tZhqthq{q{jket—YZX•VY[T(tY[eT#{qhq\8huY[tfd"ƒ(fm]€prTs]€Y[hqf]WV{‡Ÿfm]^Y[T8  V X[{jf5Y[T(ƒ•L]Lhji^kLTst d>fXY[LT tZhq\#kL{jVY[hqfm] fdQ_@Ž@]LV\8hqƒ6V{€t—Ž@t—Y[T(\8t-‘`eTsX[T)Y[LT ke]Lƒ(TsXZY•Vhq]€YcŽ:hjt¡Tsz@U@X[T(tZtZT6_ hq]8VJ]efm] U@X[fmwLVwLhq{qhqt—Y[hqƒ dgfX[\H‚ ™ T#tZef6‘¢Y[WVYoY[LT#hj]^Y[TsX•VƒsY[hqfm]Tsz@hjt—Y[hq]e}w|TsYc‘ T(Ts] Y[eT UeX[fwWVwLhq{qhqYcŽ _@hqt—YZX[hjwekeY[hqfm]HV]W_5Y[LT_€Že]LV\8hqƒ(t)fd‡Y[LTt—Ž@t—Y[T(\£\V(ŽUeX[f€_@kLƒ(TmlLhq]tZfm\8Tƒ(VtZT(t(lBV]Hhq]WVƒs” ƒskeX•VY[T)fkeY[ƒ(f\8T)fd|Y[eT9Ÿfm]^Y[T9  VX[{qf&tZhq\ kL{jV Y[hjf]Flm\V¤€hq]L}oY[eT(]8V_ep€hqt[VwL{qT)Y[LT Vm_efUeY[hqfm] fd¡V]H¥ VULU@X[fmU@X[hDV Y[TU|fmtZtZhqwLhq{qhqt—Y[hqƒ:Y[T(ƒ•L]ehji€keTm‚ {qtZfel€‘ TU@X[fmU|ftZT&V#]LTN‘V{q}mfX[hqY[e\¦dgfX9Y[eTJ]€kL\8TNX[hjƒ(V{BtZhq\ kL{jV Y[hjf]fd¡V _@hu“BTsX[T(]^Y[hjV{ t—Ž@t—Y[T(\Hly‘`LTsX[TY[LT kL]eƒ(TsXZY•Vhj]^YcŽ§fd:ULVX•V\8TNY[TsX[tHV]W_W˜ fXhj]ehuY[hDV{ƒ(fm]L_ehuY[hqfm]etHhqtX[T(U@X[Ts” tZTs]€Y[T(_/w^Ždgke¨(¨sŽ/_ehqt—YZX[hqwLk@Y[hqfm]Lt(‚"©-hq]WV{j{uŽ^l¡‘ T_@hqtZƒ(kLtZtœY[eTHtZhq\ kL{jV Y[hjf]/fd&Vž_€Ž@]WV\8hqƒ6V{ t—Ž@t—Y[T(\ª‘`LTNX[TY[eTT(prfm{qk@Y[hjf]'fd`Y[LTU|fmtZtZhqwehj{qhqt—Y[hqƒHkL]eƒ(TsXZY•Vhj]^YcŽžhqt tZhq\#kL{jVY[T6_w^ŽkLtZhq]L} w|fY[U@X[fmwWVwLhq{qhqt—Y[hjƒ:V]W_U|ftZtZhjwehq{jhqt—Y[hqƒ&Y[T(ƒNe]Lhji^kLT(t(‚ « ¬L­9®|¯B°±#²:³7®B´[°9­ µ ¶ M:?BI@KNI@?I@·E Ie?eKsC|KNI@?’OMK ¶ I|¹B¸Ñ-çZ?’PQ=@¿•KsECBAg=€¿:ÜJM ÜA`·ÁI|E CF¸I@?§K ¶ M ¸N§:™^•í ¼FMHK ¶ M¸NM KI€· PFIeCBOHMI€·9= ¿NMß@I@?ÒHìv^•í:M8¸•O5CBAg=^KsI@?Qß@?BßKNIÂñ:G&ù Ñ çZ?’î ïBÅ %€ðFÜJM:PB¿sI@P7Ie¸NM¹H=@? =€PBPQ¿NIL=@E ¶ ·ÁI@¿oK ¶ M#¿NMm¸•IeADCBKNI@?ÂI€·¡?LKNM¿NÄ^=€A7¹|ADI^Ü?Bß÷ÜJMÜAœ¹|M?BI@KNM'