(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 4.0, MathReader 4.0, or any compatible application. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. ***********************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 41513, 1324]*) (*NotebookOutlinePosition[ 42210, 1348]*) (* CellTagsIndexPosition[ 42166, 1344]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell["Absorbing boundary potential example", "Title"], Cell["\<\ This example shows you how you can use the output to keep track of certain \ aspects of the solution.\ \>", "Text"], Cell["\<\ The default for the AbsorbingBoundary option is to have an absorbing boundary \ potential of height 1 for 1/4 of the length of the interval at each end of \ the computational domain.\ \>", "Text"], Cell[BoxData[{ \(\(xmin\ = \ \(-10\)\ \[Pi];\)\), "\[IndentingNewLine]", \(\(xmax\ = \ 10\ \[Pi];\)\)}], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(Plot[ Im[\(SplitStep`Private`AbsorbingBoundaryConditionPotential[1, \ {xmin, xmax}, {xmin, xmax}\ + \ \((1/4)\)\ \((xmax\ - \ xmin)\)\ {1, \(-1\)}]\)[x]], {x, xmin, xmax}];\)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.5 0.0151576 0.603319 0.588604 [ [.04527 .59082 -9 -9 ] [.04527 .59082 9 0 ] [.19685 .59082 -9 -9 ] [.19685 .59082 9 0 ] [.34842 .59082 -9 -9 ] [.34842 .59082 9 0 ] [.65158 .59082 -6 -9 ] [.65158 .59082 6 0 ] [.80315 .59082 -6 -9 ] [.80315 .59082 6 0 ] [.95473 .59082 -6 -9 ] [.95473 .59082 6 0 ] [.4875 .01472 -12 -4.5 ] [.4875 .01472 0 4.5 ] [.4875 .13244 -24 -4.5 ] [.4875 .13244 0 4.5 ] [.4875 .25016 -24 -4.5 ] [.4875 .25016 0 4.5 ] [.4875 .36788 -24 -4.5 ] [.4875 .36788 0 4.5 ] [.4875 .4856 -24 -4.5 ] [.4875 .4856 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .04527 .60332 m .04527 .60957 L s [(-30)] .04527 .59082 0 1 Mshowa .19685 .60332 m .19685 .60957 L s [(-20)] .19685 .59082 0 1 Mshowa .34842 .60332 m .34842 .60957 L s [(-10)] .34842 .59082 0 1 Mshowa .65158 .60332 m .65158 .60957 L s [(10)] .65158 .59082 0 1 Mshowa .80315 .60332 m .80315 .60957 L s [(20)] .80315 .59082 0 1 Mshowa .95473 .60332 m .95473 .60957 L s [(30)] .95473 .59082 0 1 Mshowa .125 Mabswid .07559 .60332 m .07559 .60707 L s .1059 .60332 m .1059 .60707 L s .13622 .60332 m .13622 .60707 L s .16653 .60332 m .16653 .60707 L s .22716 .60332 m .22716 .60707 L s .25748 .60332 m .25748 .60707 L s .28779 .60332 m .28779 .60707 L s .31811 .60332 m .31811 .60707 L s .37874 .60332 m .37874 .60707 L s .40905 .60332 m .40905 .60707 L s .43937 .60332 m .43937 .60707 L s .46968 .60332 m .46968 .60707 L s .53032 .60332 m .53032 .60707 L s .56063 .60332 m .56063 .60707 L s .59095 .60332 m .59095 .60707 L s .62126 .60332 m .62126 .60707 L s .68189 .60332 m .68189 .60707 L s .71221 .60332 m .71221 .60707 L s .74252 .60332 m .74252 .60707 L s .77284 .60332 m .77284 .60707 L s .83347 .60332 m .83347 .60707 L s .86378 .60332 m .86378 .60707 L s .8941 .60332 m .8941 .60707 L s .92441 .60332 m .92441 .60707 L s .01496 .60332 m .01496 .60707 L s .98504 .60332 m .98504 .60707 L s .25 Mabswid 0 .60332 m 1 .60332 L s .5 .01472 m .50625 .01472 L s [(-1)] .4875 .01472 1 0 Mshowa .5 .13244 m .50625 .13244 L s [(-0.8)] .4875 .13244 1 0 Mshowa .5 .25016 m .50625 .25016 L s [(-0.6)] .4875 .25016 1 0 Mshowa .5 .36788 m .50625 .36788 L s [(-0.4)] .4875 .36788 1 0 Mshowa .5 .4856 m .50625 .4856 L s [(-0.2)] .4875 .4856 1 0 Mshowa .125 Mabswid .5 .04415 m .50375 .04415 L s .5 .07358 m .50375 .07358 L s .5 .10301 m .50375 .10301 L s .5 .16187 m .50375 .16187 L s .5 .1913 m .50375 .1913 L s .5 .22073 m .50375 .22073 L s .5 .27959 m .50375 .27959 L s .5 .30902 m .50375 .30902 L s .5 .33845 m .50375 .33845 L s .5 .39731 m .50375 .39731 L s .5 .42674 m .50375 .42674 L s .5 .45617 m .50375 .45617 L s .5 .51503 m .50375 .51503 L s .5 .54446 m .50375 .54446 L s .5 .57389 m .50375 .57389 L s .25 Mabswid .5 0 m .5 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .5 Mabswid .02381 .01472 m .02499 .01475 L .02605 .01484 L .02729 .01503 L .02846 .01527 L .03053 .01587 L .03279 .01678 L .03527 .01807 L .0379 .01979 L .04262 .02373 L .05205 .03491 L .06244 .05214 L .0842 .10331 L .10458 .16662 L .14335 .31094 L .18458 .46296 L .20517 .52463 L .21519 .54917 L .22428 .5678 L .2338 .58332 L .24277 .59399 L .2477 .59817 L .25224 .60093 L .25461 .60196 L .25597 .60242 L .25722 .60276 L .25853 .60303 L .25925 .60314 L .25994 .60322 L .26115 .6033 L .26184 .60332 L .26248 .60332 L .26365 .60332 L .26489 .60332 L .26711 .60332 L .2721 .60332 L .28114 .60332 L .30174 .60332 L .33985 .60332 L .38041 .60332 L .41946 .60332 L .46096 .60332 L .50094 .60332 L .5394 .60332 L .58032 .60332 L .61973 .60332 L .65761 .60332 L .69795 .60332 L .71653 .60332 L .72634 .60332 L Mistroke .73172 .60332 L .73435 .60332 L .73551 .60332 L .73678 .60332 L .73748 .60332 L .73813 .60332 L .73881 .60331 L .73955 .60326 L .74035 .60319 L .74112 .60309 L .74254 .60281 L .74384 .60247 L .74505 .60208 L .74776 .60093 L .75012 .59962 L .75271 .59787 L .75795 .59328 L .76295 .58764 L .76833 .58021 L .77805 .56335 L .79813 .5157 L .81925 .45012 L .85741 .30799 L .89802 .15787 L .91814 .0969 L .92802 .07218 L .93712 .05297 L .94624 .0374 L .95588 .02522 L .96088 .0207 L .96556 .01761 L .96824 .01633 L .97069 .01549 L .97204 .01516 L .97275 .01502 L .97351 .0149 L .97417 .01482 L .9749 .01476 L .97619 .01472 L Mfstroke % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgOol008Ioo`D000Aoo`03001oogoo08ioo`00R7oo00<007ooOol017oo 00<007ooOol0SWoo0007Ool3001nOol00`00Oomoo`04Ool30024Ool20008Ool000Uoo`8007Eoo`@0 00Aoo`03001oogoo00Aoo`03001oogoo089oo`8000Yoo`002goo00<007ooOol0NWoo00<007ooOol0 17oo00<007ooOol0PGoo00<007ooOol02Woo000Ool0011oo`03001oogoo07aoo`03001oogoo07eoo`03001oogoo00ioo`0047oo00<0 07ooOol0O7oo0P00OGoo00<007ooOol03goo000AOol00`00Oomoo`1kOol00`00Oomoo`1kOol00`00 Oomoo`0@Ool0015oo`03001oogoo07]oo`03001oogoo07]oo`03001oogoo011oo`004Woo00<007oo Ool0NWoo00<007ooOol0NWoo00<007ooOol04Goo000BOol00`00Oomoo`1jOol00`00Oomoo`1iOol0 0`00Oomoo`0BOol001=oo`03001oogoo07Uoo`03001oogoo07Uoo`03001oogoo019oo`004goo00<0 07ooOol0NGoo00<007ooOol0N7oo00<007ooOol04goo000DOol00`00Oomoo`1hOol00`00Oomoo`1h Ool00`00Oomoo`0COol001Aoo`03001oogoo07Qoo`03001oogoo07Moo`03001oogoo01Aoo`005Goo 00<007ooOol0Mgoo0P00N7oo00<007ooOol057oo000EOol00`00Oomoo`1gOol00`00Oomoo`1fOol0 0`00Oomoo`0EOol001Eoo`03001oogoo07Moo`03001oogoo07Ioo`03001oogoo01Eoo`005Woo00<0 07ooOol0MWoo00<007ooOol0MWoo00<007ooOol05Goo000FOol00`00Oomoo`1fOol00`00Oomoo`1f Ool00`00Oomoo`0EOol001Moo`03001oogoo07Eoo`03001oogoo07Eoo`03001oogoo01Ioo`005goo 00<007ooOol0MGoo00<007ooOol0MGoo00<007ooOol05Woo000HOol00`00Oomoo`1dOol00`00Oomo o`1eOol00`00Oomoo`0FOol001Qoo`03001oogoo07Aoo`8007Ioo`03001oogoo01Ioo`0067oo00<0 07ooOol0M7oo00<007ooOol0M7oo00<007ooOol05goo000IOol00`00Oomoo`1cOol00`00Oomoo`1d Ool00`00Oomoo`0GOol001Uoo`03001oogoo07=oo`03001oogoo07Aoo`03001oogoo01Moo`006Goo 00<007ooOol0Lgoo00<007ooOol0Lgoo00<007ooOol067oo000JOol00`00Oomoo`1bOol00`00Oomo o`1cOol00`00Oomoo`0HOol001Yoo`03001oogoo079oo`03001oogoo07=oo`03001oogoo01Qoo`00 6Woo00<007ooOol0GWoo0P0017oo0P0017oo0P001Woo00<007ooOol0LWoo00<007ooOol06Goo000K Ool00`00Oomoo`1LOol01000Oomoo`0027oo00@007ooOol000Eoo`03001oogoo079oo`03001oogoo 01Uoo`006goo00<007ooOol0G7oo00@007ooOol000Qoo`04001oogoo0005Ool3001aOol00`00Oomo o`0JOol001]oo`03001oogoo05Ioo`@0009oo`04001oogoo0009Ool20006Ool00`00Oomoo`1aOol0 0`00Oomoo`0JOol001aoo`03001oogoo05]oo`04001oogoo0008Ool01000Oomoo`001Goo00<007oo Ool0LGoo00<007ooOol06Woo000LOol00`00Oomoo`1LOol2000:Ool20006Ool00`00Oomoo`1`Ool0 0`00Oomoo`0KOol001aoo`03001oogoo071oo`03001oogoo071oo`03001oogoo01]oo`007Goo00<0 07ooOol0Kgoo00<007ooOol0L7oo00<007ooOol06goo000MOol00`00Oomoo`1_Ool00`00Oomoo`1_ Ool00`00Oomoo`0LOol001eoo`03001oogoo06moo`03001oogoo06moo`03001oogoo01aoo`007Woo 00<007ooOol0KWoo0P00L7oo00<007ooOol077oo000NOol00`00Oomoo`1^Ool00`00Oomoo`1^Ool0 0`00Oomoo`0MOol001ioo`03001oogoo06ioo`03001oogoo06ioo`03001oogoo01eoo`007goo00<0 07ooOol0KGoo00<007ooOol0KWoo00<007ooOol07Goo000OOol00`00Oomoo`1]Ool00`00Oomoo`1^ Ool00`00Oomoo`0MOol001moo`03001oogoo06eoo`03001oogoo06eoo`03001oogoo01ioo`007goo 00<007ooOol0KGoo00<007ooOol0KGoo00<007ooOol07Woo000POol00`00Oomoo`1/Ool00`00Oomo o`1]Ool00`00Oomoo`0NOol0021oo`03001oogoo06aoo`03001oogoo06aoo`03001oogoo01moo`00 87oo00<007ooOol0K7oo0P00KGoo00<007ooOol07goo000POol00`00Oomoo`1/Ool00`00Oomoo`1/ Ool00`00Oomoo`0OOol0025oo`03001oogoo06]oo`03001oogoo06aoo`03001oogoo01moo`008Goo 00<007ooOol0Jgoo00<007ooOol0Jgoo00<007ooOol087oo000QOol00`00Oomoo`1[Ool00`00Oomo o`1[Ool00`00Oomoo`0POol0025oo`03001oogoo06]oo`03001oogoo06]oo`03001oogoo021oo`00 8Woo00<007ooOol0JWoo00<007ooOol0Jgoo00<007ooOol087oo000ROol00`00Oomoo`1ZOol00`00 Oomoo`1ZOol00`00Oomoo`0QOol0029oo`03001oogoo06Yoo`8006]oo`03001oogoo025oo`008Woo 00<007ooOol0JWoo00<007ooOol0JWoo00<007ooOol08Goo000SOol00`00Oomoo`1YOol00`00Oomo o`1YOol00`00Oomoo`0ROol002=oo`03001oogoo06Uoo`03001oogoo06Uoo`03001oogoo029oo`00 8goo00<007ooOol0JGoo00<007ooOol0JGoo00<007ooOol08Woo000TOol00`00Oomoo`1XOol00`00 Oomoo`1YOol00`00Oomoo`0ROol002Aoo`03001oogoo05Aoo`8000Aoo`8000=oo`<000Ioo`03001o ogoo06Qoo`03001oogoo02=oo`0097oo00<007ooOol0Dgoo00@007ooOol000Qoo`04001oogoo0005 Ool00`00Oomoo`1XOol00`00Oomoo`0SOol002Aoo`03001oogoo05=oo`04001oogoo0008Ool01000 Oomoo`001Goo0`00J7oo00<007ooOol08goo000UOol00`00Oomoo`17oo00<007ooOol0>Goo10000Woo00@007ooOol0 00Yoo`03001oogoo00Aoo`03001oogoo05Aoo`03001oogoo03Moo`00>7oo00<007ooOol0?goo00@0 07ooOol000Qoo`04001oogoo0005Ool00`00Oomoo`1DOol00`00Oomoo`0gOol003Qoo`03001oogoo 041oo`8000Yoo`8000Ioo`03001oogoo05Aoo`03001oogoo03Moo`00>Goo00<007ooOol0Dgoo00<0 07ooOol0Dgoo00<007ooOol0>7oo000iOol00`00Oomoo`1COol00`00Oomoo`1COol00`00Oomoo`0h Ool003Uoo`03001oogoo05=oo`03001oogoo05=oo`03001oogoo03Qoo`00>Woo00<007ooOol0DWoo 00<007ooOol0DWoo00<007ooOol0>Goo000jOol00`00Oomoo`1BOol2001COol00`00Oomoo`0iOol0 03Yoo`03001oogoo059oo`03001oogoo059oo`03001oogoo03Uoo`00>goo00<007ooOol0DGoo00<0 07ooOol0DGoo00<007ooOol0>Woo000kOol00`00Oomoo`1AOol00`00Oomoo`1AOol00`00Oomoo`0j Ool003]oo`03001oogoo055oo`03001oogoo051oo`03001oogoo03]oo`00?7oo00<007ooOol0D7oo 00<007ooOol0D7oo00<007ooOol0>goo000lOol00`00Oomoo`1@Ool00`00Oomoo`1?Ool00`00Oomo o`0lOol003eoo`03001oogoo04moo`03001oogoo04moo`03001oogoo03aoo`00?Goo00<007ooOol0 Cgoo0P00D7oo00<007ooOol0?7oo000nOol00`00Oomoo`1>Ool00`00Oomoo`1>Ool00`00Oomoo`0m Ool003ioo`03001oogoo04ioo`03001oogoo04ioo`03001oogoo03eoo`00?Woo00<007ooOol0CWoo 00<007ooOol0CGoo00<007ooOol0?Woo000oOol00`00Oomoo`1=Ool00`00Oomoo`1=Ool00`00Oomo o`0nOol003moo`03001oogoo04eoo`03001oogoo04aoo`03001oogoo03moo`0037oo0P0017oo0P00 8Woo10000goo0P0000=oo`00Ool087oo1@000Woo0P0097oo00<007ooOol08goo1@000Woo0P0087oo 00<007ooOol010000goo0P008goo0P0017oo0P002goo000;Ool01000Oomoo`000Woo00@007ooOol0 025oo`03001oogoo00=oo`03001oogoo008002=oo`05001oogooOol00002Ool00`00Oomoo`0QOol0 0`00Oomoo`0UOol01@00Oomoogoo00000Woo00<007ooOol07Goo00@007ooOol000Eoo`04001oogoo 000QOol01000Oomoo`000Woo00@007ooOol000Yoo`003Goo00D007ooOomoo`00009oo`03001oogoo 021oo`03001oogoo009oo`06001oogoo001oo`008Woo00D007ooOomoo`00009oo`03001oogoo025o o`03001oogoo02Eoo`05001oogooOol00002Ool00`00Oomoo`0LOol00`00Oomoo`02Ool00`00Oomo o`02Ool01000Oomoo`008goo00D007ooOomoo`00009oo`03001oogoo00Qoo`001Goo10001Goo00@0 07ooOol0009oo`03001oogoo01Uoo`@000Aoo`05001oogooOol00002Ool00`00Ool0000JOol40004 Ool01@00Oomoogoo00000Woo00<007ooOol08Goo0P009Woo00D007ooOomoo`00009oo`03001oogoo 01aoo`03001oogoo00=oo`05001oogooOol00002Ool00`00Oomoo`0ROol01000Oomoo`000Woo00<0 07ooOol027oo000;Ool01000Oomoo`000Woo00@007ooOol0025oo`04001oogoo0002Ool01000Oomo o`000Woo00<007ooOol07goo00D007ooOomoo`00009oo`03001oogoo025oo`03001oogoo02Eoo`05 001oogooOol00002Ool00`00Oomoo`0KOol00`00Oomoo`02Ool01000Oomoo`000Woo00@007ooOol0 025oo`04001oogoo0002Ool01000Oomoo`002Woo000Ool008moo`03001oogoo08ioo`00\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-33.1234, -1.02675, 0.230825 , 0.00594415}}] }, Open ]], Cell[BoxData[ \(sol\ = \ SplitStepSolve[\(Soliton[1, 1]\)[x, 0], {x, \(-10\)\ \[Pi], \ 10\ \[Pi], \ 180}, {{0, 20, 1}, Automatic}, \ AbsorbingBoundary \[Rule] True, \ Output \[Rule] {Mass, \ {Plot, \ PlotRange \[Rule] All}}]\)], "Input"], Cell["\<\ By specifying the Mass in the output, you can see how this ordinarily \ conserved quantities decays due to the absorbing potential.\ \>", "Text"], Cell[CellGroupData[{ Cell[BoxData[ \(\(ListPlot[{Time, \ Log[10, Mass]}\ /. \ sol];\)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.047619 0.550167 0.176568 [ [.2619 .53767 -3 -9 ] [.2619 .53767 3 0 ] [.5 .53767 -6 -9 ] [.5 .53767 6 0 ] [.7381 .53767 -6 -9 ] [.7381 .53767 6 0 ] [.97619 .53767 -6 -9 ] [.97619 .53767 6 0 ] [.01131 .02046 -12 -4.5 ] [.01131 .02046 0 4.5 ] [.01131 .10875 -24 -4.5 ] [.01131 .10875 0 4.5 ] [.01131 .19703 -12 -4.5 ] [.01131 .19703 0 4.5 ] [.01131 .28531 -24 -4.5 ] [.01131 .28531 0 4.5 ] [.01131 .3736 -12 -4.5 ] [.01131 .3736 0 4.5 ] [.01131 .46188 -24 -4.5 ] [.01131 .46188 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .2619 .55017 m .2619 .55642 L s [(5)] .2619 .53767 0 1 Mshowa .5 .55017 m .5 .55642 L s [(10)] .5 .53767 0 1 Mshowa .7381 .55017 m .7381 .55642 L s [(15)] .7381 .53767 0 1 Mshowa .97619 .55017 m .97619 .55642 L s [(20)] .97619 .53767 0 1 Mshowa .125 Mabswid .07143 .55017 m .07143 .55392 L s .11905 .55017 m .11905 .55392 L s .16667 .55017 m .16667 .55392 L s .21429 .55017 m .21429 .55392 L s .30952 .55017 m .30952 .55392 L s .35714 .55017 m .35714 .55392 L s .40476 .55017 m .40476 .55392 L s .45238 .55017 m .45238 .55392 L s .54762 .55017 m .54762 .55392 L s .59524 .55017 m .59524 .55392 L s .64286 .55017 m .64286 .55392 L s .69048 .55017 m .69048 .55392 L s .78571 .55017 m .78571 .55392 L s .83333 .55017 m .83333 .55392 L s .88095 .55017 m .88095 .55392 L s .92857 .55017 m .92857 .55392 L s .25 Mabswid 0 .55017 m 1 .55017 L s .02381 .02046 m .03006 .02046 L s [(-3)] .01131 .02046 1 0 Mshowa .02381 .10875 m .03006 .10875 L s [(-2.5)] .01131 .10875 1 0 Mshowa .02381 .19703 m .03006 .19703 L s [(-2)] .01131 .19703 1 0 Mshowa .02381 .28531 m .03006 .28531 L s [(-1.5)] .01131 .28531 1 0 Mshowa .02381 .3736 m .03006 .3736 L s [(-1)] .01131 .3736 1 0 Mshowa .02381 .46188 m .03006 .46188 L s [(-0.5)] .01131 .46188 1 0 Mshowa .125 Mabswid .02381 .03812 m .02756 .03812 L s .02381 .05578 m .02756 .05578 L s .02381 .07343 m .02756 .07343 L s .02381 .09109 m .02756 .09109 L s .02381 .1264 m .02756 .1264 L s .02381 .14406 m .02756 .14406 L s .02381 .16172 m .02756 .16172 L s .02381 .17937 m .02756 .17937 L s .02381 .21469 m .02756 .21469 L s .02381 .23234 m .02756 .23234 L s .02381 .25 m .02756 .25 L s .02381 .26766 m .02756 .26766 L s .02381 .30297 m .02756 .30297 L s .02381 .32063 m .02756 .32063 L s .02381 .33829 m .02756 .33829 L s .02381 .35594 m .02756 .35594 L s .02381 .39126 m .02756 .39126 L s .02381 .40891 m .02756 .40891 L s .02381 .42657 m .02756 .42657 L s .02381 .44423 m .02756 .44423 L s .02381 .47954 m .02756 .47954 L s .02381 .4972 m .02756 .4972 L s .02381 .51485 m .02756 .51485 L s .02381 .53251 m .02756 .53251 L s .02381 .00281 m .02756 .00281 L s .02381 .56782 m .02756 .56782 L s .02381 .58548 m .02756 .58548 L s .02381 .60314 m .02756 .60314 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .008 w .02381 .60332 Mdot .07143 .60332 Mdot .11905 .60332 Mdot .16667 .60332 Mdot .21429 .60332 Mdot .2619 .60332 Mdot .30952 .60332 Mdot .35714 .60331 Mdot .40476 .60304 Mdot .45238 .59903 Mdot .5 .58246 Mdot .54762 .54553 Mdot .59524 .48564 Mdot .64286 .40682 Mdot .69048 .31814 Mdot .7381 .2307 Mdot .78571 .15503 Mdot .83333 .09755 Mdot .88095 .05802 Mdot .92857 .03196 Mdot .97619 .01472 Mdot % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg1oo`80021oo`006goo00<007ooOol0ogoo0goo000KOol00`00Oomoo`3oOol3Ool001]o o`800?moo`Aoo`006goo00<007ooOol0ogoo0goo000KOol00`00Oomoo`3oOol3Ool001]oo`03001o ogoo0?moo`=oo`006goo00<007ooOol0ogoo0goo000KOol2003oOol4Ool001]oo`03001oogoo0=Ao o`8002aoo`006goo00<007ooOol0e7oo0P00;7oo0007Ool40003Ool20004Ool20005Ool00`00Oomo o`3oOol3Ool000Moo`03001oogoo00Uoo`04001oogoo0004Ool00`00Oomoo`3oOol3Ool000Qoo`03 001oogoo00]oo`03001oogoo009oo`800?moo`Aoo`0000=oo`0000000P0017oo00<007ooOol027oo 0P001Goo00<007ooOol0ogoo0goo0007Ool01000Oomoo`002Goo00<007ooOol017oo00<007ooOol0 ogoo0goo0008Ool2000:Ool30004Ool00`00Oomoo`3oOol3Ool001]oo`800?moo`Aoo`006goo00<0 07ooOol0ogoo0goo000KOol00`00Oomoo`3oOol3Ool001]oo`03001oogoo0?moo`=oo`006goo00<0 07ooOol0ogoo0goo000KOol2003oOol4Ool001]oo`03001oogoo0?moo`=oo`006goo00<007ooOol0 ogoo0goo000KOol00`00Oomoo`37Ool2000iOol001]oo`03001oogoo07oo00D007ooOomoo`00009oo`03001oogoo005oo`006goo00<007ooOol0 >goo00<007ooOol0>Woo00D007ooOomoo`00009oo`03001oogoo03Ioo`03001oogoo009oo`03001o ogoo03Eoo`04001oogoo0002Ool01000Oomoo`000goo000KOol00`00Oomoo`0kOol3000iOol20004 Ool2000hOol20004Ool3000fOol20004Ool20004Ool001]oo`800?moo`Aoo`006goo00<007ooOol0 ogoo0goo000KOol00`00Oomoo`3oOol3Ool001]oo`03001oogoo0?moo`=oo`006goo00<007ooOol0 R7oo0P00N7oo000DOooo000=000001]oo`03001oogoo00Uoo`03001oogoo00Yoo`03001oogoo00Yo o`03001oogoo00Uoo`03001oogoo00Yoo`03001oogoo00Yoo`03001oogoo00Uoo`03001oogoo00Yo o`03001oogoo00Yoo`03001oogoo00Uoo`03001oogoo00Yoo`03001oogoo00Yoo`03001oogoo00Uo o`03001oogoo00Yoo`03001oogoo00Yoo`03001oogoo00Uoo`03001oogoo00Yoo`03001oogoo00Yo o`03001oogoo00Uoo`03001oogoo00Yoo`03001oogoo00Eoo`006goo00<007ooOol0?7oo00<007oo Ool0?7oo00<007ooOol0?Goo00<007ooOol0?7oo00<007ooOol01Goo000KOol00`00Oomoo`3oOol3 Ool001]oo`03001oogoo0?moo`=oo`006goo0P00ogoo17oo000KOol00`00Oomoo`3oOol3Ool001]o o`03001oogoo0?moo`=oo`006goo00<007ooOol0Ngoo0P00QGoo000KOol00`00Oomoo`1kOol20025 Ool001]oo`800?moo`Aoo`006goo00<007ooOol0ogoo0goo000KOol00`00Oomoo`3oOol3Ool001]o o`03001oogoo06ioo`80099oo`006Woo0`002Woo0P002goo0P002Woo0P002goo0P002goo0P002Woo 0P002goo0P002goo0P002Woo0P00TWoo000JOol2000;Ool2000;Ool2000:Ool2000;Ool2000;Ool2 000:Ool2000;Ool2000;Ool2002NOol001]oo`03001oogoo0?moo`=oo`006goo00<007ooOol0ogoo 0goo000KOol00`00Oomoo`3oOol3Ool00?moob5oo`00ogoo8Goo003oOolQOol00?moob5oo`00ogoo 8Goo003oOolQOol00001\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-2.15729, -3.24937, 0.0789457 , 0.021291}}] }, Open ]], Cell["\<\ You can get a better overall reduction by balancing the reflection and \ transmission. In this example, the computation is stopped when the Mass \ becomes less than 0.0005. The value of the StoppingCondition option is \ assumed to be a function which is applied to the current output.\ \>", "Text"], Cell[BoxData[ \(\(sol\ = \ SplitStepSolve[\(Soliton[1, 1]\)[x, 0], {x, \(-10\)\ \[Pi], \ 10\ \[Pi], \ 180}, {{0, 30, 1}, Automatic}, \ AbsorbingBoundary \[Rule] { .5, 1/3}, \ Output \[Rule] {Mass, \ {Plot, \ PlotRange \[Rule] All}}, \ StoppingCondition \[Rule] Function[\((Mass\ /. \ #)\)\ < \ .0005]];\)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(ListPlot[{Time, \ Log[10, Mass]}\ /. \ sol];\)\)], "Input"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.04329 0.554446 0.162353 [ [.24026 .54195 -3 -9 ] [.24026 .54195 3 0 ] [.45671 .54195 -6 -9 ] [.45671 .54195 6 0 ] [.67316 .54195 -6 -9 ] [.67316 .54195 6 0 ] [.88961 .54195 -6 -9 ] [.88961 .54195 6 0 ] [.01131 .06739 -12 -4.5 ] [.01131 .06739 0 4.5 ] [.01131 .14856 -24 -4.5 ] [.01131 .14856 0 4.5 ] [.01131 .22974 -12 -4.5 ] [.01131 .22974 0 4.5 ] [.01131 .31092 -24 -4.5 ] [.01131 .31092 0 4.5 ] [.01131 .39209 -12 -4.5 ] [.01131 .39209 0 4.5 ] [.01131 .47327 -24 -4.5 ] [.01131 .47327 0 4.5 ] [ 0 0 0 0 ] [ 1 .61803 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 g .25 Mabswid [ ] 0 setdash .24026 .55445 m .24026 .5607 L s [(5)] .24026 .54195 0 1 Mshowa .45671 .55445 m .45671 .5607 L s [(10)] .45671 .54195 0 1 Mshowa .67316 .55445 m .67316 .5607 L s [(15)] .67316 .54195 0 1 Mshowa .88961 .55445 m .88961 .5607 L s [(20)] .88961 .54195 0 1 Mshowa .125 Mabswid .0671 .55445 m .0671 .5582 L s .11039 .55445 m .11039 .5582 L s .15368 .55445 m .15368 .5582 L s .19697 .55445 m .19697 .5582 L s .28355 .55445 m .28355 .5582 L s .32684 .55445 m .32684 .5582 L s .37013 .55445 m .37013 .5582 L s .41342 .55445 m .41342 .5582 L s .5 .55445 m .5 .5582 L s .54329 .55445 m .54329 .5582 L s .58658 .55445 m .58658 .5582 L s .62987 .55445 m .62987 .5582 L s .71645 .55445 m .71645 .5582 L s .75974 .55445 m .75974 .5582 L s .80303 .55445 m .80303 .5582 L s .84632 .55445 m .84632 .5582 L s .9329 .55445 m .9329 .5582 L s .97619 .55445 m .97619 .5582 L s .25 Mabswid 0 .55445 m 1 .55445 L s .02381 .06739 m .03006 .06739 L s [(-3)] .01131 .06739 1 0 Mshowa .02381 .14856 m .03006 .14856 L s [(-2.5)] .01131 .14856 1 0 Mshowa .02381 .22974 m .03006 .22974 L s [(-2)] .01131 .22974 1 0 Mshowa .02381 .31092 m .03006 .31092 L s [(-1.5)] .01131 .31092 1 0 Mshowa .02381 .39209 m .03006 .39209 L s [(-1)] .01131 .39209 1 0 Mshowa .02381 .47327 m .03006 .47327 L s [(-0.5)] .01131 .47327 1 0 Mshowa .125 Mabswid .02381 .08362 m .02756 .08362 L s .02381 .09986 m .02756 .09986 L s .02381 .11609 m .02756 .11609 L s .02381 .13233 m .02756 .13233 L s .02381 .1648 m .02756 .1648 L s .02381 .18103 m .02756 .18103 L s .02381 .19727 m .02756 .19727 L s .02381 .2135 m .02756 .2135 L s .02381 .24597 m .02756 .24597 L s .02381 .26221 m .02756 .26221 L s .02381 .27845 m .02756 .27845 L s .02381 .29468 m .02756 .29468 L s .02381 .32715 m .02756 .32715 L s .02381 .34339 m .02756 .34339 L s .02381 .35962 m .02756 .35962 L s .02381 .37586 m .02756 .37586 L s .02381 .40833 m .02756 .40833 L s .02381 .42456 m .02756 .42456 L s .02381 .4408 m .02756 .4408 L s .02381 .45703 m .02756 .45703 L s .02381 .4895 m .02756 .4895 L s .02381 .50574 m .02756 .50574 L s .02381 .52198 m .02756 .52198 L s .02381 .53821 m .02756 .53821 L s .02381 .05115 m .02756 .05115 L s .02381 .03492 m .02756 .03492 L s .02381 .01868 m .02756 .01868 L s .02381 .00245 m .02756 .00245 L s .02381 .57068 m .02756 .57068 L s .02381 .58692 m .02756 .58692 L s .02381 .60315 m .02756 .60315 L s .25 Mabswid .02381 0 m .02381 .61803 L s 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath .008 w .02381 .60332 Mdot .0671 .60332 Mdot .11039 .60332 Mdot .15368 .60332 Mdot .19697 .60332 Mdot .24026 .6033 Mdot .28355 .60284 Mdot .32684 .59995 Mdot .37013 .59192 Mdot .41342 .57646 Mdot .45671 .55211 Mdot .5 .51845 Mdot .54329 .47615 Mdot .58658 .42674 Mdot .62987 .37222 Mdot .67316 .3147 Mdot .71645 .25621 Mdot .75974 .1988 Mdot .80303 .14483 Mdot .84632 .09721 Mdot .88961 .05901 Mdot .9329 .03186 Mdot .97619 .01472 Mdot % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgioo`80019oo`006goo0P00kgoo0P004Woo000KOol00`00Oomo o`3oOol3Ool001]oo`03001oogoo0?moo`=oo`006goo00<007ooOol0ogoo0goo000KOol2003oOol4 Ool001]oo`03001oogoo0?moo`=oo`006goo00<007ooOol0hgoo0P007Goo000DOol20005Ool00`00 Oomoo`3SOol2000MOol001=oo`04001oogoo0004Ool00`00Oomoo`3oOol3Ool001Eoo`03001oogoo 00=oo`800?moo`Aoo`003Goo10001Goo00<007ooOol00Woo00<007ooOol0ogoo0goo000COol01000 Oomoo`0017oo00<007ooOol0ogoo0goo000DOol20005Ool00`00Oomoo`3oOol3Ool001]oo`800?mo o`Aoo`006goo00<007ooOol0ogoo0goo000KOol00`00Oomoo`3oOol3Ool001]oo`03001oogoo0=Mo o`8002Uoo`006goo0P00f7oo0P00:Goo000KOol00`00Oomoo`3oOol3Ool001]oo`03001oogoo0?mo o`=oo`006goo00<007ooOol0ogoo0goo000KOol00`00Oomoo`3oOol3Ool001]oo`800?moo`Aoo`00 6goo00<007ooOol0ogoo0goo000KOol00`00Oomoo`3oOol3Ool001]oo`03001oogoo0?moo`=oo`00 6goo0P00ogoo17oo000KOol00`00Oomoo`3oOol3Ool000Moo`@000=oo`8000Aoo`8000Eoo`03001o ogoo0?moo`=oo`001goo00<007ooOol02Goo00@007ooOol000Aoo`03001oogoo0"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-2.37302, -3.56022, 0.0868403 , 0.0231552}}] }, Open ]] }, Open ]] }, FrontEndVersion->"4.0 for Microsoft Windows", ScreenRectangle->{{0, 1024}, {0, 695}}, WindowSize->{496, 599}, WindowMargins->{{183, Automatic}, {23, Automatic}}, StyleDefinitions -> "ArticleClassic.nb" ] (*********************************************************************** Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. ***********************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[1739, 51, 53, 0, 111, "Title"], Cell[1795, 53, 125, 3, 44, "Text"], Cell[1923, 58, 206, 4, 61, "Text"], Cell[2132, 64, 124, 2, 51, "Input"], Cell[CellGroupData[{ Cell[2281, 70, 275, 6, 131, "Input"], Cell[2559, 78, 16722, 483, 186, 4690, 330, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[19296, 564, 276, 5, 91, "Input"], Cell[19575, 571, 155, 3, 44, "Text"], Cell[CellGroupData[{ Cell[19755, 578, 82, 1, 31, "Input"], Cell[19840, 581, 10295, 349, 186, 3715, 263, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[30150, 933, 311, 5, 78, "Text"], Cell[30464, 940, 391, 7, 111, "Input"], Cell[CellGroupData[{ Cell[30880, 951, 82, 1, 31, "Input"], Cell[30965, 954, 10520, 366, 186, 3923, 280, "GraphicsData", "PostScript", "Graphics"] }, Open ]] }, Open ]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)