(*********************************************************************** Mathematica-Compatible Notebook This notebook can be used on any computer system with Mathematica 3.0, MathReader 3.0, or any compatible application. The data for the notebook starts with the line of 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[ 132232, 3656]*) (*NotebookOutlinePosition[ 133006, 3683]*) (* CellTagsIndexPosition[ 132962, 3679]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell["Numerical Methods for Conservation Laws.", "Title"], Cell[CellGroupData[{ Cell["Page 29", "Section"], Cell[BoxData[ \(\(viscousShockPlot[wl_, wr_, \[Epsilon]_] := Plot[wr + \(\((wl - wr)\)\/2\) \((1 - Tanh[\(\((wl - wr)\)\/\(4\ \[Epsilon]\)\) \[Theta]\ ])\), { \[Theta], \(-3\), 3}]; \)\)], "Input"], Cell[BoxData[ \(\(Table[viscousShockPlot[2, 0, 1\/10\^i], {i, 0, 3}]; \)\)], "Input"] }, Closed]], Cell[CellGroupData[{ Cell["Page 72", "Section"], Cell[CellGroupData[{ Cell[BoxData[ \(Solve[{m\&~ - m\&^ == s \((\[Rho]\&~ - \[Rho]\&^)\), \n\t\t \((m\&~\^2\/\[Rho]\&~ + \(a\^2\) \[Rho]\&~)\) - \((m\&^\^2\/\[Rho]\&^ + \(a\^2\) \[Rho]\&^)\) == s \((m\&~ - m\&^)\)}, \n\t{m\&~, s}] // Simplify\)], "Input"], Cell[BoxData[ \({{m\&~ \[Rule] \(a\ \[Rho]\&^\^\(3/2\)\ \@\[Rho]\&~ + m\&^\ \[Rho]\&~ - a\ \@\[Rho]\&^\ \[Rho]\&~\^\(3/2\)\)\/\[Rho]\&^, s \[Rule] \(m\&^ - a\ \@\[Rho]\&^\ \@\[Rho]\&~\)\/\[Rho]\&^}, { m\&~ \[Rule] \(\(-a\)\ \[Rho]\&^\^\(3/2\)\ \@\[Rho]\&~ + m\&^\ \[Rho]\&~ + a\ \@\[Rho]\&^\ \[Rho]\&~\^\(3/2\)\)\/\[Rho]\&^, s \[Rule] \(m\&^ + a\ \@\[Rho]\&^\ \@\[Rho]\&~\)\/\[Rho]\&^}}\)], "Output"] }, Open ]], Cell[BoxData[ \({{m\&~ \[Rule] \(m\&^\ \[Rho]\&~\)\/\[Rho]\&^ + \((a\ \((\[Rho]\&^ - \[Rho]\&~)\)\ \@\[Rho]\&~)\)\/\@\[Rho]\&^, s \[Rule] m\&^\/\[Rho]\&^ - \(a\ \@\[Rho]\&~\)\/\@\[Rho]\&^}, { m\&~ \[Rule] \(m\&^\ \[Rho]\&~\)\/\[Rho]\&^ - \((a\ \((\[Rho]\&^ - \[Rho]\&~)\)\ \@\[Rho]\&~)\)\/\@\[Rho]\&^, s \[Rule] m\&^\/\[Rho]\&^ + \(a\ \@\[Rho]\&~\)\/\@\[Rho]\&^}}\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(Eigensystem[{{v, 1}, {\[Phi], v}}]\)], "Input"], Cell[BoxData[ \({{v - \@\[Phi], v + \@\[Phi]}, {{\(-\(1\/\@\[Phi]\)\), 1}, {1\/\@\[Phi], 1}}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Solve[\n \t{\(v\&^\)\^2\/2 + \(\[Phi]\&^\) + \(-\(v\&~\)\^2\)\/2 - \(\[Phi]\&~\) == s \((\(v\&^\) - \(v\&~\))\), \n \(v\&^\)\ \(\[Phi]\&^\) - \(v\&~\)\ \(\[Phi]\&~\) == s \((\(\[Phi]\&^\) - \(\[Phi]\&~\))\)}, \n{\(v\&~\), s}] // Simplify\)], "Input"], Cell[BoxData[ \({{\(v\&~\) \[Rule] \(v\&^\) + \(\@2\ \((\(\[Phi]\&^\) - \(\[Phi]\&~\))\)\)\/\@\(\(\[Phi]\&^\) + \(\[Phi]\&~\)\), s \[Rule] \(v\&^\) - \(\@2\ \(\[Phi]\&~\)\)\/\@\(\(\[Phi]\&^\) + \(\[Phi]\&~\)\)}, { \(v\&~\) \[Rule] \(v\&^\) + \(\@2\ \((\(-\(\[Phi]\&^\)\) + \(\[Phi]\&~\)) \)\)\/\@\(\(\[Phi]\&^\) + \(\[Phi]\&~\)\), s \[Rule] \(v\&^\) + \(\@2\ \(\[Phi]\&~\)\)\/\@\(\(\[Phi]\&^\) + \(\[Phi]\&~\)\)}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(Out[2] /. {\(\[Phi]\&~\) -> \(\[Phi]\&^\) \((1 + \[Xi])\)} // Simplify \) // PowerExpand\)], "Input"], Cell[BoxData[ \({{\(v\&~\) \[Rule] \(v\&^\) - \(\@2\ \[Xi]\ \@\(\[Phi]\&^\)\)\/\@\(2 + \[Xi]\), s \[Rule] \(v\&^\) - \(\@2\ \((1 + \[Xi])\)\ \@\(\[Phi]\&^\)\)\/\@\(2 + \[Xi]\)}, { \(v\&~\) \[Rule] \(v\&^\) + \(\@2\ \[Xi]\ \@\(\[Phi]\&^\)\)\/\@\(2 + \[Xi]\), s \[Rule] \(v\&^\) + \(\@2\ \((1 + \[Xi])\)\ \@\(\[Phi]\&^\)\)\/\@\(2 + \[Xi]\)}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(Plot[\[Rho] + \((\[Rho] - 1)\) \@\[Rho], {\[Rho], 0, 3}]; \)\)], "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 /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.31746 0.0311077 0.0885214 [ [.18254 .01861 -9 -9 ] [.18254 .01861 9 0 ] [.34127 .01861 -3 -9 ] [.34127 .01861 3 0 ] [.5 .01861 -9 -9 ] [.5 .01861 9 0 ] [.65873 .01861 -3 -9 ] [.65873 .01861 3 0 ] [.81746 .01861 -9 -9 ] [.81746 .01861 9 0 ] [.97619 .01861 -3 -9 ] [.97619 .01861 3 0 ] [.01131 .11963 -6 -4.5 ] [.01131 .11963 0 4.5 ] [.01131 .20815 -6 -4.5 ] [.01131 .20815 0 4.5 ] [.01131 .29667 -6 -4.5 ] [.01131 .29667 0 4.5 ] [.01131 .38519 -6 -4.5 ] [.01131 .38519 0 4.5 ] [.01131 .47371 -6 -4.5 ] [.01131 .47371 0 4.5 ] [.01131 .56224 -6 -4.5 ] [.01131 .56224 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 .18254 .03111 m .18254 .03736 L s [(0.5)] .18254 .01861 0 1 Mshowa .34127 .03111 m .34127 .03736 L s [(1)] .34127 .01861 0 1 Mshowa .5 .03111 m .5 .03736 L s [(1.5)] .5 .01861 0 1 Mshowa .65873 .03111 m .65873 .03736 L s [(2)] .65873 .01861 0 1 Mshowa .81746 .03111 m .81746 .03736 L s [(2.5)] .81746 .01861 0 1 Mshowa .97619 .03111 m .97619 .03736 L s [(3)] .97619 .01861 0 1 Mshowa .125 Mabswid .05556 .03111 m .05556 .03486 L s .0873 .03111 m .0873 .03486 L s .11905 .03111 m .11905 .03486 L s .15079 .03111 m .15079 .03486 L s .21429 .03111 m .21429 .03486 L s .24603 .03111 m .24603 .03486 L s .27778 .03111 m .27778 .03486 L s .30952 .03111 m .30952 .03486 L s .37302 .03111 m .37302 .03486 L s .40476 .03111 m .40476 .03486 L s .43651 .03111 m .43651 .03486 L s .46825 .03111 m .46825 .03486 L s .53175 .03111 m .53175 .03486 L s .56349 .03111 m .56349 .03486 L s .59524 .03111 m .59524 .03486 L s .62698 .03111 m .62698 .03486 L s .69048 .03111 m .69048 .03486 L s .72222 .03111 m .72222 .03486 L s .75397 .03111 m .75397 .03486 L s .78571 .03111 m .78571 .03486 L s .84921 .03111 m .84921 .03486 L s .88095 .03111 m .88095 .03486 L s .9127 .03111 m .9127 .03486 L s .94444 .03111 m .94444 .03486 L s .25 Mabswid 0 .03111 m 1 .03111 L s .02381 .11963 m .03006 .11963 L s [(1)] .01131 .11963 1 0 Mshowa .02381 .20815 m .03006 .20815 L s [(2)] .01131 .20815 1 0 Mshowa .02381 .29667 m .03006 .29667 L s [(3)] .01131 .29667 1 0 Mshowa .02381 .38519 m .03006 .38519 L s [(4)] .01131 .38519 1 0 Mshowa .02381 .47371 m .03006 .47371 L s [(5)] .01131 .47371 1 0 Mshowa .02381 .56224 m .03006 .56224 L s [(6)] .01131 .56224 1 0 Mshowa .125 Mabswid .02381 .04881 m .02756 .04881 L s .02381 .06652 m .02756 .06652 L s .02381 .08422 m .02756 .08422 L s .02381 .10192 m .02756 .10192 L s .02381 .13733 m .02756 .13733 L s .02381 .15504 m .02756 .15504 L s .02381 .17274 m .02756 .17274 L s .02381 .19045 m .02756 .19045 L s .02381 .22585 m .02756 .22585 L s .02381 .24356 m .02756 .24356 L s .02381 .26126 m .02756 .26126 L s .02381 .27897 m .02756 .27897 L s .02381 .31438 m .02756 .31438 L s .02381 .33208 m .02756 .33208 L s .02381 .34978 m .02756 .34978 L s .02381 .36749 m .02756 .36749 L s .02381 .4029 m .02756 .4029 L s .02381 .4206 m .02756 .4206 L s .02381 .43831 m .02756 .43831 L s .02381 .45601 m .02756 .45601 L s .02381 .49142 m .02756 .49142 L s .02381 .50912 m .02756 .50912 L s .02381 .52683 m .02756 .52683 L s .02381 .54453 m .02756 .54453 L s .02381 .0134 m .02756 .0134 L s .02381 .57994 m .02756 .57994 L s .02381 .59764 m .02756 .59764 L s .02381 .61535 m .02756 .61535 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 .5 Mabswid .02381 .03108 m .02499 .02607 L .02605 .02434 L .02729 .02291 L .02846 .02185 L .03053 .02037 L .03279 .01915 L .03527 .01809 L .0379 .01722 L .04262 .01608 L .04496 .01568 L .04749 .01534 L .04987 .01509 L .05205 .01493 L .05337 .01485 L .0546 .0148 L .05595 .01476 L .0567 .01474 L .05738 .01473 L .05863 .01472 L .05934 .01472 L .05999 .01472 L .06118 .01473 L .06244 .01476 L .06309 .01478 L .06379 .0148 L .06507 .01485 L .06746 .01497 L .07025 .01515 L .07283 .01536 L .08253 .01645 L .09322 .01812 L .10458 .02034 L .14233 .03026 L .18253 .04407 L .22122 .05976 L .26236 .07855 L .30199 .09842 L .34009 .11897 L .38065 .14225 L .4197 .16592 L .4612 .19232 L .50118 .2189 L .53965 .24546 L .58056 .27472 L .61997 .30384 L .65785 .33266 L .69819 .36421 L .73702 .39538 L .7783 .42936 L Mistroke .81806 .46287 L .8563 .4958 L .897 .53159 L .93618 .56674 L .97619 .60332 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]HGAYHf4PAg9QL6QYHgc/lC4a000<`00IP00V@00c000o`0c000c <`0cIP0cV@0cc00co`1V001V<`1VIP1VV@1Vc01Vo`2I002I<`2IIP2IV@2Ic02Io`3<003<<`3I03>IIIS>IVC>Ic3>Ioc?<03?<[ZjWAdM:NWYm?Cdl30`:0Pl:j^[WMgMolG0?l_0?m70?mO0?mg0?n?0?nW 0?no0?oG0?o_0?Oo0=oo0@000`40000g00030@0000/0104j00030@0000X0 000A00030@0000@01`4200/170000`40000200030@0000X000<10000>P000`40000g00030@0000`0 00<10000>P000`40000g00030@0000/000<10000>`000`40000:00004@000`4000040A@01@4H0004 0@000@d000<10000=`040CH0104>00030@0003D000@100000P4d00040@0000813@000`40000g0005 0@0000403000014000<100407@040AD00P4>00D1=`020CP00P4>00D1=0040CH0104>00D1=P030@d0 000A00818P040Ol0N00001400P4V00@1o`1d00004@000`40000Y00@1o`1`00001`3o0JL10`000140 00<100002P000`40000:00030@0000/000<100002@020@`000<100002`000`40000:00030@0000/0 00<100002P000`40000:00030@0000/000<100002P000`40000:00030@0000/000<100002P000`40 000:00030@0000/000<100002P000`40000;00030@0000X000<100002P000`40000;00030@0000X0 00<100002P000`40000;00030@0000X000<100002`000`40000:00030@0000X000<100002`000`40 000;00004@000`40000:00030@0000X000<100002`000`40000:00@12@000`40000;00030@0000X0 00<100002`000`40000:00030@0000X000<100002`000`40000:00030@0000X000<100002`000`40 000:00030@0000X000<100002`000`40000:00030@0000/000<100002P000`40000:00030@0000/0 00<100002P000`40000:00030@0000/000<100002P000`40000;00030@0000X000<100002P000`40 000;00030@0000/0000A00030@0003H00`4600030@00040000<10000@0000`40001000030@000400 00<10000@0000`40000;00004@000`40000i0081o`1R00004@000`40000k00<1o`1O00004@000`40 000n00<1o`1L00004@000`40001100<1o`1I00004@020DD00`7o05H0000A00030@0004L00P7o05@0 000A00030@0004T00`7o0540000A00030@0004`00`7o04h0000A00030@0004l00`7o04/0000A0003 0@0005800P7o04T0000A00030@0005@00P7o04L0000A0081E`020Ol0A@00014000<10000F0020Ol0 @`00014000<10000FP020Ol0@@00014000<10000G0020Ol0?`00014000<10000GP020Ol0?@000140 00<10000H0020Ol0>`00014000<10000HP020Ol0>@00014000<10000I0020Ol0=`0001400P5W0081 o`0e00004@000`40001X0081o`0c00004@000`40001Z0081o`0a00004@000`40001/0081o`0_0000 4@000`40001^0081o`0]00004@000`40001`0081o`0[00004@000`40001b00030@000?l0:0000140 0P5d0081o`0X00004@000`40001e0081o`0V00004@000`40001g0081o`0T00004@000`40001i0081 o`0R00004@000`40001k0081o`0P000010060@L000<10000O@020Ol07P0000L000<100001`000`40 001o0081o`0L00001`000`40000700030@0008400P7o01X0000700030@0000L00`630081o`0H0000 1`000`40000700030@0008D00P7o01H0000700030@0000L000<10000Q`000`40003o01<000070003 0@0000L000<10000R0020Ol04`0000L000<100001`000`40002:0081o`0A000010040@T000<10000 S0020Ol03`0000H00P4900030@0008h000<10000o`0<00004@020I000P7o00`0000A00030@000940 0P7o00X0000A00030@0009<000<10000o`0700004@000`40002D0081o`0700004@000`40002F0003 0@000?l01000014000<10000U`020Ol01000014000<10000V@020Ol00P00014000<10000V`000`40 003n00004@020Id00P7n00004@000`40002N00030@000?/0000A00030@0009l00P7k00004@000`40 002Q0081n@00014000<10000X`000`40003f00004@000`40002T0081mP00014000<10000YP020O@0 000A0081Z@000`40003a00004@000`40002Y0081l@00014000<10000Z`000`40003^00004@000`40 002/0081kP00014000<10000[P020N`0000A00030@000;0000<10000j@00014000<10000/@020NT0 000A0081]0020NL0000A00030@000;D000<10000i000014000<10000]P020N@0000A00030@000;P0 00<10000h@00014000<10000^@020N40000300H120000`40002k00030@000=h0000300810`000`40 000600030@000;`00P7N000010020@/000<10000_P020M`0000500812P030L0000<10000f@0000H0 0P4900030@000<400P7I00001`020@P000<10000``000`40003F000020000`40000600030@000<@0 0P7F000020000`40000600030@0000000 4@000`40003<0081cP00014000<10000cP000`40003;00004@000`40003?00030@0000000<10000^@00014000<10000h@020KT0000A00030@000><000<10000 ]P00014000<10000i0000`40002e00004@000`40003U0081]@00014000<10000i`000`40002b0000 4@000`40003X0081/P0001400P7[00030@000:l0000A00030@000>/000<10000[P00014000<10000 k0020Jh0000A00030@000>h000<10000Z`00014000<10000k`020J/0000400@12@000`40003a0003 0@000:P0000300040@00008120000`40003b0081Z00000P000<100001P030O@000<10000Y@0000P0 00<100001P000`40003e00030@000:@00007008120000`40003f0081Y00000D00`4900030@000?P0 00<10000X@0000P000<100001P000`40003i00030@000:00000800030@0000H000<10000nP020J00 000400050@00004020000`40003l00030@0009d0000500<12@000`40003m00030@0009`0000A0081 o`020I`0000A00030@000?l00@000`40002I00004@000`40003o008000<10000V000014000<10000 o`030081V000014000<10000o`0500030@0009D0000A00030@000?l01P000`40002D00004@000`40 003o00L00P6D00004@020Ol02P000`40002A00004@000`40003o00X000<10000T000014000<10000 o`0;0081T000014000<10000o`0=00030@0008d0000A00030@000?l03P000`40002<00004@000`40 003o00l000<10000R`00014000<10000o`0@00030@0008X0000A00030@000?l04@020HX0000A0081 o`0D00030@0008L0000A00030@000?l050000`40002600004@000`40003o01D000<10000Q@000140 00<10000o`0F00030@0008@0000A00030@000?l05`020H@0000A00030@000?l06@000`4000210000 4@000`40003o01X000<10000P00001400P7o01`000<10000O`00014000<10000o`0L0081O`000140 00<10000o`0N00030@0007`0000A00030@000?l07`000`40001k00004@000`40003o020000<10000 NP0000H00`4800030@000?l08@000`40001i00001`000`40000700030@000?l08P020GT0000300H1 20030Ol090000`40001f00000`001@40000100T000<10000o`0U00030@0007D0000400040@000@T0 00<10000o`0V00030@0007@0000400040@000@T000<10000o`0W0081M00000D000<100402@000`40 003o02T000<10000L@0000D000<100402@000`40003o02X000<10000L00000H00P4900030@000?l0 :`000`40001_00001P020@T000<10000o`0/0081K`0001400P7o02l000<10000K000014000<10000 o`0_00030@0006/0000A00030@000?l0<0000`40001Z00004@000`40003o034000<10000J@000140 00<10000o`0b0081J@00014000<10000o`0d00030@0006H0000A00030@000?l0=@000`40001U0000 4@020Ol0=`000`40001T00004@000`40003o03L000<10000H`00014000<10000o`0h0081H`000140 00<10000o`0j00030@000600000A00030@000?l0>`000`40001O00004@000`40003o03`000<10000 GP00014000<10000o`0m0081GP00014000<10000o`0o00030@0005/0000A0081o`1100030@0005X0 000A00030@000?l0@@000`40001I00004@000`40003o048000<10000F000014000<10000o`130081 F000014000<10000o`1500030@0005D0000A00030@000?l0AP000`40001D00004@000`40003o04L0 00<10000D`0001400P7o04T000<10000DP00014000<10000o`190081DP00014000<10000o`1;0003 0@0004l0000A00030@000?l0C0000`40001>00004@000`40003o04d000<10000C@0000@010490003 0@000?l0CP000`40001<000010001@40000100P000<10000o`1?00030@0004/0000900030@0000D0 00<10000o`1@0081B`0000T000<100001@030Ol0DP000`40001800002@000`40000500030@000?l0 D`000`400017000010001@40000100P000<10000o`1D00030@0004H0000400@12@000`40003o05D0 00<10000A@0000@000<100002P000`40003o05H000<10000A00000@000<100002P000`40003o05L0 00<10000@`0000@01@4800030@000?l0F0020D<0000A0081o`1K00030@000400000A00030@000?l0 F`000`40000o00004@000`40003o05`000<10000?P00014000<10000o`1M00030@0003d0000A0003 0@000?l0GP000`40000l00004@000`40003o05l000<10000>`00014000<10000o`1P0081>`000140 00<10000o`1R00030@0003P0000A0081o`1T00030@0003L0000A00030@000?l0I0000`40000f0000 4@000`40003o06D000<10000=@00014000<10000o`1V00030@0003@0000A00030@000?l0I`000`40 000c00004@000`40003o06P00P4c00004@000`40003o06X000<10000<00001400P7o06`000<10000 ;`00014000<10000o`1/00030@0002h0000A00030@000?l0K@000`40000]00004@000`40003o06h0 00<10000;000014000<10000o`1_00030@0002/0000A00030@000?l0L0000`40000Z00004@000`40 003o07400P4Z00004@000`40003o07<000<100009`0001400P7o07D000<100009P00014000<10000 o`1e00030@0002D0000A00030@000?l0MP000`40000T00004@000`40003o07L000<100008`000140 00<10000o`1h00030@000280000600<120000`40003o07T00P4R00001@000`4000020@L000<10000 o`1k00030@0001l0000400030@00008000<100001@030Ol0O0000`40000N000010000`4000020003 0@0000D000<10000o`1m00030@0001d0000400810P020@L000<10000o`1n00030@0001`000040003 0@01008120000`40003o07l000<100006`0000@000<100002P000`40003o080000<100006P0000D0 00<100002@000`40003o084000<100006@0000D00P4:00030@000?l0PP020AT0000700<11`020Ol0 Q@000`40000F00004@000`40003o08D000<100005@00014000<10000o`2600030@0001@0000A0003 0@000?l0Q`000`40000C00004@000`40003o08P000<100004P00014000<10000o`2900030@000140 000A00030@000?l0RP020A40000A00030@000?l0S0000`40000>00004@020Ol0SP000`40000=0000 4@000`40003o08h000<100003000014000<10000o`2?00030@0000/0000A00030@000?l0W@000140 00<10000o`2M00004@000`40003o09d0000A00030@000?l0W@0001400P7o09h0000A00030@000?l0 W@000001\ \>"], ImageRangeCache->{{{0, 431}, {265.875, 0}} -> {-0.128159, -0.503439, 0.0112225, 0.0402469}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(Plot[{\[Rho] - \[Rho]\ Log[\[Rho]], \n\t\t \[Rho] + \[Rho]\ Log[\[Rho]]}, {\[Rho], 10\^\(-2\), 3}]; \)\)], "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 /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.31746 0.0729201 0.196747 [ [.18254 .06042 -9 -9 ] [.18254 .06042 9 0 ] [.34127 .06042 -3 -9 ] [.34127 .06042 3 0 ] [.5 .06042 -9 -9 ] [.5 .06042 9 0 ] [.65873 .06042 -3 -9 ] [.65873 .06042 3 0 ] [.81746 .06042 -9 -9 ] [.81746 .06042 9 0 ] [.97619 .06042 -3 -9 ] [.97619 .06042 3 0 ] [.01131 .17129 -18 -4.5 ] [.01131 .17129 0 4.5 ] [.01131 .26967 -6 -4.5 ] [.01131 .26967 0 4.5 ] [.01131 .36804 -18 -4.5 ] [.01131 .36804 0 4.5 ] [.01131 .46641 -6 -4.5 ] [.01131 .46641 0 4.5 ] [.01131 .56479 -18 -4.5 ] [.01131 .56479 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 .18254 .07292 m .18254 .07917 L s [(0.5)] .18254 .06042 0 1 Mshowa .34127 .07292 m .34127 .07917 L s [(1)] .34127 .06042 0 1 Mshowa .5 .07292 m .5 .07917 L s [(1.5)] .5 .06042 0 1 Mshowa .65873 .07292 m .65873 .07917 L s [(2)] .65873 .06042 0 1 Mshowa .81746 .07292 m .81746 .07917 L s [(2.5)] .81746 .06042 0 1 Mshowa .97619 .07292 m .97619 .07917 L s [(3)] .97619 .06042 0 1 Mshowa .125 Mabswid .05556 .07292 m .05556 .07667 L s .0873 .07292 m .0873 .07667 L s .11905 .07292 m .11905 .07667 L s .15079 .07292 m .15079 .07667 L s .21429 .07292 m .21429 .07667 L s .24603 .07292 m .24603 .07667 L s .27778 .07292 m .27778 .07667 L s .30952 .07292 m .30952 .07667 L s .37302 .07292 m .37302 .07667 L s .40476 .07292 m .40476 .07667 L s .43651 .07292 m .43651 .07667 L s .46825 .07292 m .46825 .07667 L s .53175 .07292 m .53175 .07667 L s .56349 .07292 m .56349 .07667 L s .59524 .07292 m .59524 .07667 L s .62698 .07292 m .62698 .07667 L s .69048 .07292 m .69048 .07667 L s .72222 .07292 m .72222 .07667 L s .75397 .07292 m .75397 .07667 L s .78571 .07292 m .78571 .07667 L s .84921 .07292 m .84921 .07667 L s .88095 .07292 m .88095 .07667 L s .9127 .07292 m .9127 .07667 L s .94444 .07292 m .94444 .07667 L s .25 Mabswid 0 .07292 m 1 .07292 L s .02381 .17129 m .03006 .17129 L s [(0.5)] .01131 .17129 1 0 Mshowa .02381 .26967 m .03006 .26967 L s [(1)] .01131 .26967 1 0 Mshowa .02381 .36804 m .03006 .36804 L s [(1.5)] .01131 .36804 1 0 Mshowa .02381 .46641 m .03006 .46641 L s [(2)] .01131 .46641 1 0 Mshowa .02381 .56479 m .03006 .56479 L s [(2.5)] .01131 .56479 1 0 Mshowa .125 Mabswid .02381 .09259 m .02756 .09259 L s .02381 .11227 m .02756 .11227 L s .02381 .13194 m .02756 .13194 L s .02381 .15162 m .02756 .15162 L s .02381 .19097 m .02756 .19097 L s .02381 .21064 m .02756 .21064 L s .02381 .23032 m .02756 .23032 L s .02381 .24999 m .02756 .24999 L s .02381 .28934 m .02756 .28934 L s .02381 .30902 m .02756 .30902 L s .02381 .32869 m .02756 .32869 L s .02381 .34837 m .02756 .34837 L s .02381 .38772 m .02756 .38772 L s .02381 .40739 m .02756 .40739 L s .02381 .42707 m .02756 .42707 L s .02381 .44674 m .02756 .44674 L s .02381 .48609 m .02756 .48609 L s .02381 .50576 m .02756 .50576 L s .02381 .52544 m .02756 .52544 L s .02381 .54511 m .02756 .54511 L s .02381 .05325 m .02756 .05325 L s .02381 .03357 m .02756 .03357 L s .02381 .0139 m .02756 .0139 L s .02381 .58446 m .02756 .58446 L s .02381 .60414 m .02756 .60414 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 .5 Mabswid .02698 .08395 m .03593 .10496 L .04573 .12282 L .06549 .1512 L .08578 .17407 L .10749 .19393 L .14786 .22204 L .16909 .23334 L .18917 .24224 L .22745 .25516 L .24722 .26003 L .26818 .264 L .27756 .26541 L .28765 .26669 L .29711 .26767 L .30589 .2684 L .31596 .26902 L .32128 .26927 L .32401 .26937 L .32696 .26946 L .32977 .26954 L .33236 .26959 L .33484 .26963 L .33609 .26964 L .33746 .26965 L .3387 .26966 L .33983 .26967 L .34114 .26967 L .34238 .26967 L .34353 .26966 L .34458 .26966 L .34574 .26965 L .34698 .26964 L .34935 .2696 L .35185 .26956 L .35632 .26945 L .36113 .26929 L .36638 .26907 L .37582 .26854 L .3846 .26791 L .40499 .26595 L .42352 .26357 L .46337 .25669 L .50171 .24795 L .54249 .23656 L .58176 .22371 L .61952 .20974 L .65973 .19324 L .69842 .17586 L .73956 .15588 L Mistroke .77919 .13525 L .81731 .11418 L .85787 .09052 L .89692 .06659 L .93841 .03997 L .97619 .01472 L Mfstroke .02698 .06583 m .02922 .06262 L .03162 .05983 L .03593 .0559 L .04106 .05247 L .04573 .05019 L .0506 .04847 L .05335 .04775 L .05589 .04723 L .05811 .04688 L .06055 .04659 L .06178 .04648 L .06246 .04643 L .0631 .04639 L .06435 .04634 L .06549 .04631 L .06614 .0463 L .06685 .04629 L .06813 .04631 L .06938 .04634 L .07053 .04639 L .07178 .04647 L .07314 .04657 L .0759 .04686 L .07832 .04718 L .0809 .0476 L .08557 .04853 L .09703 .05173 L .10749 .05563 L .12651 .06474 L .1472 .07712 L .18784 .10746 L .22698 .14264 L .2646 .1809 L .30467 .22567 L .34323 .2721 L .38423 .32465 L .42373 .378 L .46171 .4316 L .50213 .49089 L .54105 .54996 L .58241 .61474 L s .58241 .61474 m .58445 .61803 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgc/lC4a000<`00IP00V@00c000o`0c000c <`0cIP0cV@0cc00co`1V001V<`1VIP1VV@1Vc01Vo`2I002I<`2IIP2IV@2Ic02Io`3<003<<`3I03>IIIS>IVC>Ic3>Ioc?<03?<[ZjWAdM:NWYm?Cdl30`:0Pl:j^[WMgMolG0?l_0?m70?mO0?mg0?n?0?nW 0?no0?oG0?o_0?Oo0=oo0000080000`40003o08h0000P00030@000?l0SP00020000<10000 o`2>000080000`40003o08h0000P00030@000?l0SP0002000P7o08@000<100002000020000<10000 o`2100812`00020000<10000o`2000030@0000/0000P00030@000?l0OP020@h0000P00030@000?l0 O@000`40000>000080000`40003o07/00P4A000080000`40003o07X000<100004@00020000<10000 o`1i00030@000180000P0081o`1h00815@00020000<10000o`1f00030@0001D0000P00030@000?l0 M0020AP0000P00030@0003@00P4600<11@040C@01P4b00H110030@D0104c00H1<@060@D00`4500@1 9`000`40000:00@12P00020000<10000<`001040004500<11@001@40000103H000<10000=@000`40 000400<11@001@40000103800P4300030@0002l00P4300030@0000<00`4500050@00004090020@`0 00@100000P49000080000`40000;00P17`000`40000200030@0000800`4:00030@0003<000<10000 =@000`40000400<12P000`40000`0081=@020@P00`4:00030@00020000<100004@000`4000070000 80000`40000700@120060AT000<100000P000`40000?00030@0003<000<10000=@000`40000A0003 0@0003400P4e008150000`40000N008150000`400007000080000`40000500814P030AH000<10000 0P000`40000?00030@0003<000<10000=@000`40000A00030@0003800P4e00814`000`40000M0003 0@0001<00P49000080020@@00P4G008150000`40000200030@0000X000D100000@0f00030@0003D0 00<1000030001@40000103H00P4e00813@001@40000101h00P4D00<12P00020000<100000P000`40 000I00814P000`40000200030@0000X0104g00030@0003D000<1000030040CP000<10000=0000`40 000;00@17@020AT000<100001`00020000D100000@0N008140000`40000200030@0000X000<10000 >0000`40000e00030@0000`000<10000>@000`40000d00030@0000/000<100007@000`40000I0003 0@0000L0000P00040@000B400P4?00040@000@d000<10000=@040C@0104>00030@0003@000@10000 0P4a00040@0000813@000`40000K008160001@40000100T0000P00030@0102@00P4>00813P050CD0 0P4f00813P050C<0104c00@13P050AP000<100006@030@X0000P00819`020Ol0>`020BT0000P0003 0@0002P00P7o03P000<10000:@00020000<10000:P000`40003o03<00P4/00005P3o0IX10@000200 00<100002P000`40000:00030@0000X000<100001P020@03004000/000<100002P000`40000:0003 0@0000X000<100002P000`40000:00030@0000X000<100002P000`40000:00030@0000X000<10000 2P000`40000:00030@0000X000<100002P000`40000:00030@0000X000<100002P000`40000:0003 0@0000X000<100002P000`40000:00030@0000X000<100002P020@/000<100002P000`40000:0003 0@0000P0000P00030@0000X000<100002P000`40000:00030@0000P00P4<00030@0000X000<10000 2P000`40000:00030@0000X000<100002P000`40000:00030@0000X000<100002P000`40000:0003 0@0000X000<100002P000`40000:00030@0000X000<100002P000`40000:00030@0000X000<10000 2P000`40000:00030@0000X000<100002P000`40000:00030@0000P00`4<00030@0000X000<10000 2P000`400008000080000`40000`00812`000`40000n00030@0003h000<10000?P000`40000n0003 0@0001@000<100009`000`400008000080000`40000b00030@000?l08P020CD0000P0081=0000`40 003o01l00P4g000080000`400@0d0081o`0M0081>@00020000<10040=P000`40003o01T000<10000 >@0002000P400`01000e0081o`0G0081?000020000@10001>0000`40003o01800P4n000080001@40 000103P00P7o014000<10000?P00020000D100000@0j00030@000?l030020D40000P00030@000080 00<10000>0000`40003o00X000<10000@@00020000<100000P000`40000i0081o`080081A0000200 00<100000`000`40000j00030@000?l00`020DH0000P00030@0000<000<10000>`000`40003o0040 00<10000AP0002000P4500030@0003/000<10000o@020DT0000P00030@0000@000<10000?0000`40 003k00030@0004T0000P00030@0000D000<10000?0000`40003h0081C000020000<100001@000`40 000m00030@000?D00P5>000080000`40000600030@0003d000<10000lP020E00000P00030@0000L0 00<10000?@000`40003_0081DP00020000<100001`000`40000n0081kP000`40001B000080000`40 000800030@0003l000<10000j@020ED0000P00812P000`40000o00030@000>H00P5G000080000`40 000900030@00040000<10000h`020ET0000P00030@0000X000<10000@0000`40003P0081F`000200 00<100002`000`40001000030@000=d00P5M000080000`40000;00030@00044000<10000fP020El0 000P00030@0000`000<10000@@000`40003G0081H@00020000<100003@000`40001100030@000=@0 0P5S000080000`40000=00030@00048000<10000d@020FD0000P00813`000`40001200030@00000D11`000`40 000L00030@0004<000<10000W`020HL0000P00817P000`40001200030@0009d00P69000080000`40 000N00030@00048000<10000V@030H/0000P00030@0001l000<10000@P000`40002F0081SP000200 00<1000080020D<000<10000T`020I00000P00030@00028000<10000@@000`40002?00<1TP000200 00<100008`020D8000<10000S0020ID0000P00030@0002D000<10000?`000`40002:0081U`000200 00<100009P000`40000o00030@0008L00P6I000080020BP00P5000030@0008<00`6K000080000`40 000Y00030@0003h000<10000P0020Ih0000P00030@0002X00P4o00030@0007`00`6P000080000`40 000/00030@0003d000<10000N@020J<0000P00030@0002d00P4m00030@0007H00`6U000080000`40 000_00030@0003/000<10000L`020JP0000P00030@0003000P4l00030@0006l00`6Z000080000`40 000b0081>`000`40001[00<1[@0002000P4e0081>P000`40001V00@1/000020000<10000=P020CT0 00<10000HP030K@0000P00030@0003P00P4h00030@0005h00`6g000080000`40000j0081=`000`40 001J00<1^P00020000<10000?0030C@000<10000EP040Kd0000P00030@0003l00`4b00030@000580 0`71000080000`4000120081<@000`40001>00<1a000020000<10000A0030Bl000<10000B@040LL0 000P0081B0030Bd000<10000A0040L/0000P00030@0004X00P4/00030@0003d01P7?000080000`40 001<00<1:P000`40000g00D1e@00020000<10000C`030BP000<10000<@050MX0000P00030@000580 104U00030@0002X01P7O000080000`40001F00H17`000`40000U00D1i@0001@01P4600030@0005`0 1P4J00030@0001d01`7Z00005`000`40000600030@000680204C00030@0001@0207a00005`000`40 000600<1JP0Z0OT0000G00030@0000H000<10000O`000`40003o00`0000G00030@0000H000<10000 P0000`40003o00/0000G00030@0000H000<10000P@000`40003o00X0000G00030@0000H000<10000 P@000`40003o00X0000G00030@0000H000<10000PP000`40003o00T0000D00@120000`4000230003 0@000?l0200001H00P4800030@0008@000<10000o`07000080000`40002400030@000?l01`000200 0P6600030@000?l01P00020000<10000QP000`40003o00D0000P00030@0008L000<10000o`040000 80000`40002800030@000?l00`00020000<10000R0000`40003o00<0000P00030@0008T000<10000 o`02000080000`40002:00030@000?l00@00020000<10000R`000`40003o000080020H`000<10000 o`00020000<10000S0000`40003n000080000`40002=00030@000?d0000P00030@0008h000<10000 o000020000<10000SP000`40003l000080000`40002?00030@000?/0000P00030@00090000<10000 nP00020000<10000T@000`40003i000080020I8000<10000n@00020000<10000TP000`40003h0000 80000`40002C00030@000?L0000P00030@0009@000<10000mP00020000<10000U0000`40003f0000 80000`40002E00030@000?D0000P00030@0009H000<10000m000020000<10000U`000`40003c0000 80020IP000<10000l`00020000<10000V0000`40003b000080000`40002I00030@000?40000P0003 0@0009T000<10000l@00020000<10000VP000`40003`00000P060@@00`4500@120000`40002K0003 0@000>l0000500030@0000@00`4500050@0000401`000`40002L00030@000>h0000500030@0000@0 0`4:00030@0000@000<10000W0000`40003^00001@000`40000A00030@0000@00`6M00030@000>d0 000500030@00014000<1000010000`40002N00030@000>`0000500030@0000`000D100000@070003 0@0009l000<10000j`0000D000<1000030040@P000<10000W`000`40003[00001@000`40000<0003 0@0000T000<10000X0000`40003Z00000P040@h000<100002@000`40002Q00030@000>T000040081 3P050@L000<10000X@000`40003Y000080000`40002R00030@000>P0000P0081Y0000`40003W0000 80000`40002T00030@000>H0000P00030@000:@000<10000iP00020000<10000Y@000`40003U0000 80000`40002V00030@000>@0000P00030@000:L000<10000h`00020000<10000Y`000`40003S0000 80000`40002X00030@000>80000P0081ZP000`40003Q000080000`40002Y00030@000>40000P0003 0@000:X000<10000h000020000<10000Z`000`40003O000080000`40002/00030@000=h0000P0003 0@000:`000<10000gP00020000<10000[@000`40003M000080000`40002^00030@000=`0000P0081 /0000`40003K000080000`40002_00030@000=/0000P00030@000;0000<10000fP00020000<10000 /@000`40003I000080000`40002a00030@000=T0000P00030@000;8000<10000f000020000<10000 /`000`40003G000080000`40002c00030@000=L0000P0081]@000`40003F000080000`40002e0003 0@000=D0000P00030@000;D000<10000e@00020000<10000]P000`40003D000080000`40002g0003 0@000=<0000C00H11`000`40002g00030@000=<0000C00810`000`40000500030@000;P000<10000 dP0001@00P4:00030@000;T000<10000d@0001D00P4900<1^@000`40003A00005P020@P000<10000 ^P000`40003@00005`020@L000<10000^`000`40003?000060000`40000500030@000;/000<10000 c`0001P000<100001@000`40002l00030@00000 00<10000ZP00020000<10000h@000`40002Y000080000`40003R00030@000:P0000P00030@000>80 00<10000Z0000?l0/@000?l0/@000?l0/@000?l0/@000?l0/@000?l0/@000001\ \>"], ImageRangeCache->{{{0, 431}, {265.875, 0}} -> {-0.24753, -0.453577, 0.0115634, 0.018658}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(DSolve[{\(\[Rho]'\)[\[Xi]] == \[Rho][\[Xi]]\/a, \(m'\)[\[Xi]] == m[\[Xi]]\/a + \[Rho][\[Xi]], \n\t\t \[Rho][\[Xi]\_1] == rowl, m[\[Xi]\_1] == mleft}, {\[Rho][\[Xi]], m[\[Xi]]}, \[Xi]]\)], "Input"], Cell[BoxData[ \(DSolve::"nvld" \( : \ \) "The description of the equations appears to be ambiguous or \ invalid."\)], "Message"], Cell[BoxData[ \({{\[Rho][\[Xi]] \[Rule] E\^\(\[Xi]\/a - \[Xi]\_1\/a\)\ rowl, m[\[Xi]] \[Rule] E\^\(\[Xi]\/a\)\ \(( E\^\(-\(\[Xi]\_1\/a\)\)\ rowl\ \[Xi] + E\^\(-\(\[Xi]\_1\/a\)\)\ \((mleft - rowl\ \[Xi]\_1)\))\)}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Solve[ \(-a\)\ Log[\[Rho]\_m\/\[Rho]\_l] + m\_l\/\[Rho]\_l == a\ Log[\[Rho]\_m\/\[Rho]\_r] + m\_r\/\[Rho]\_r, \[Rho]\_m] // Simplify\)], "Input"], Cell[BoxData[ \({{\[Rho]\_m \[Rule] \(-E\^\(\(m\_l\/\[Rho]\_l - m\_r\/\[Rho]\_r\)\/\(2\ a\)\)\)\ \@\[Rho]\_l\ \@\[Rho]\_r}, { \[Rho]\_m \[Rule] E\^\(\(m\_l\/\[Rho]\_l - m\_r\/\[Rho]\_r\)\/\(2\ a\)\)\ \@\[Rho]\_l\ \@\[Rho]\_r}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(m = {{0, 1, 0}, \n \t\t{\(1\/2\) \((\[Gamma] - 3)\) v\^2, \((3 - \[Gamma])\) v, \((\[Gamma] - 1)\)}, \n \t\t{\(1\/2\) \((\[Gamma] - 1)\) v\^3 - v \((e + p)\)\/\[Rho], \((e + p)\)\/\[Rho] - \((\[Gamma] - 1)\) v\^2, \[Gamma]\ v}} /. {e -> p\/\(\[Gamma] - 1\) + \(1\/2\) \[Rho]\ v\^2} // Simplify; \nEigenvectors[m] // Simplify\)], "Input"], Cell[BoxData[ \({{2\/v\^2, 2\/v, 1}, { \(-\(\(2\ \((\(-1\) + \[Gamma])\)\ \[Rho]\)\/\(\(-2\)\ p\ \[Gamma] + 2\ \@p\ v\ \((\(-1\) + \[Gamma])\)\ \@\[Gamma]\ \@\[Rho] - v\^2\ \((\(-1\) + \[Gamma])\)\ \[Rho]\)\)\), \(2\ \((\(-1\) + \[Gamma])\)\ \((\@p\ \@\[Gamma] - v\ \@\[Rho])\)\ \@\[Rho]\)\/\(\(-2\)\ p\ \[Gamma] + 2\ \@p\ v\ \((\(-1\) + \[Gamma])\)\ \@\[Gamma]\ \@\[Rho] - v\^2\ \((\(-1\) + \[Gamma])\)\ \[Rho]\), 1}, { \(2\ \((\(-1\) + \[Gamma])\)\ \[Rho]\)\/\(2\ p\ \[Gamma] + 2\ \@p\ v\ \((\(-1\) + \[Gamma])\)\ \@\[Gamma]\ \@\[Rho] + v\^2\ \((\(-1\) + \[Gamma])\)\ \[Rho]\), \(2\ \((\(-1\) + \[Gamma])\)\ \((\@p\ \@\[Gamma] + v\ \@\[Rho])\)\ \@\[Rho]\)\/\(2\ p\ \[Gamma] + 2\ \@p\ v\ \((\(-1\) + \[Gamma])\)\ \@\[Gamma]\ \@\[Rho] + v\^2\ \((\(-1\) + \[Gamma])\)\ \[Rho]\), 1}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Out[30] /. {e -> p\/\(\[Gamma] - 1\) + \(1\/2\) \[Rho]\ v\^2} // Simplify\)], "Input"], Cell[BoxData[ \({v, v - \(p\ \[Gamma]\)\/\@\(p\ \[Gamma]\ \[Rho]\), v + \(p\ \[Gamma]\)\/\@\(p\ \[Gamma]\ \[Rho]\)}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(PowerExpand[Out[32]] // Simplify\)], "Input"], Cell[BoxData[ \({v, v - \(\@p\ \@\[Gamma]\)\/\@\[Rho], v + \(\@p\ \@\[Gamma]\)\/\@\[Rho]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Solve[ \((\(\[Rho]\_m\) m\_l\/\[Rho]\_l - \(\@\(\[Rho]\_m\/\[Rho]\_l\)\) \((\[Rho]\_m - \[Rho]\_l)\) == \(\[Rho]\_m\) m\_r\/\[Rho]\_r + \(\@\(\[Rho]\_m\/\[Rho]\_r\)\) \((\[Rho]\_m - \[Rho]\_r)\)) \) /. {\[Rho]\_l -> .28260, m\_l -> .098185, \[Rho]\_r -> .2, m\_r -> 0}, {\[Rho]\_m}]\)], "Input"], Cell[BoxData[ \({{\[Rho]\_m \[Rule] 0.`}, {\[Rho]\_m \[Rule] 0.282599328618680765`}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\[Rho]\_m\) m\_l\/\[Rho]\_l - \(\@\(\[Rho]\_m\/\[Rho]\_l\)\) \((\[Rho]\_m - \[Rho]\_l)\) /. Join[{\[Rho]\_l -> .28260, m\_l -> .098185, \[Rho]\_r -> .2, m\_r -> 0}, \(Out[2]\)[\([2]\)]]\)], "Input"], Cell[BoxData[ \(0.0981854381194643188`\)], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Page 79", "Section"], Cell[CellGroupData[{ Cell[BoxData[ \(Eigensystem[{{v, 1}, {\[Phi], v}}]\)], "Input"], Cell[BoxData[ \({{v - \@\[Phi], v + \@\[Phi]}, {{\(-\(1\/\@\[Phi]\)\), 1}, {1\/\@\[Phi], 1}}}\)], "Output"] }, Open ]], Cell[BoxData[ \(\(eqs = { \(v\&~\)\^2\/2 + \(\[Phi]\&~\) - \(v\&^\)\^2\/2 - \(\[Phi]\&^\) == s \((\(v\&~\) - \(v\&^\))\), \(\[Phi]\&~\)\ \(v\&~\) - \(\[Phi]\&^\)\ \(v\&^\) == s \((\(\[Phi]\&~\) - \(\[Phi]\&^\))\)}; \)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(Solve[eqs, {\(v\&~\), \(\[Phi]\&~\)}] // Simplify\)], "Input"], Cell[BoxData[ \({{\(\[Phi]\&~\) \[Rule] \(\[Phi]\&^\), \(v\&~\) \[Rule] \(v\&^\)}, { \(\[Phi]\&~\) \[Rule] 1\/4\ \((s - \(v\&^\))\)\ \((s - \(v\&^\) + \@\(s\^2 - 2\ s\ \(v\&^\) + \(v\&^\)\^2 + 8\ \(\[Phi]\&^\)\)) \), \(v\&~\) \[Rule] 1\/2\ \(( 3\ s - \(v\&^\) - \@\(s\^2 - 2\ s\ \(v\&^\) + \(v\&^\)\^2 + 8\ \(\[Phi]\&^\)\)) \)}, {\(\[Phi]\&~\) \[Rule] \(-\(1\/4\)\)\ \((s - \(v\&^\))\)\ \((\(-s\) + \(v\&^\) + \@\(s\^2 - 2\ s\ \(v\&^\) + \(v\&^\)\^2 + 8\ \(\[Phi]\&^\)\)) \), \(v\&~\) \[Rule] 1\/2\ \(( 3\ s - \(v\&^\) + \@\(s\^2 - 2\ s\ \(v\&^\) + \(v\&^\)\^2 + 8\ \(\[Phi]\&^\)\)) \)}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Solve[eqs, {\(v\&~\), s}] // Simplify\)], "Input"], Cell[BoxData[ \({{\(v\&~\) \[Rule] \(v\&^\) + \(\@2\ \((\(\[Phi]\&^\) - \(\[Phi]\&~\))\)\)\/\@\(\(\[Phi]\&^\) + \(\[Phi]\&~\)\), s \[Rule] \(v\&^\) - \(\@2\ \(\[Phi]\&~\)\)\/\@\(\(\[Phi]\&^\) + \(\[Phi]\&~\)\)}, { \(v\&~\) \[Rule] \(v\&^\) + \(\@2\ \((\(-\(\[Phi]\&^\)\) + \(\[Phi]\&~\)) \)\)\/\@\(\(\[Phi]\&^\) + \(\[Phi]\&~\)\), s \[Rule] \(v\&^\) + \(\@2\ \(\[Phi]\&~\)\)\/\@\(\(\[Phi]\&^\) + \(\[Phi]\&~\)\)}}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Solve[eqs, {\(\[Phi]\&~\), s}] // Simplify\)], "Input"], Cell[BoxData[ \({{\(\[Phi]\&~\) \[Rule] 1\/4\ \(( \(v\&^\)\^2 + 4\ \(\[Phi]\&^\) + \(v\&~\)\ \((\(v\&~\) - \@\(\(v\&^\)\^2 + 16\ \(\[Phi]\&^\) - 2\ \(v\&^\)\ \(v\&~\) + \(v\&~\)\^2\))\) + \(v\&^\)\ \((\(-2\)\ \(v\&~\) + \@\(\(v\&^\)\^2 + 16\ \(\[Phi]\&^\) - 2\ \(v\&^\)\ \(v\&~\) + \(v\&~\)\^2\))\))\), s \[Rule] 1\/4\ \(( \(v\&^\) + 3\ \(v\&~\) - \@\(\(v\&^\)\^2 + 16\ \(\[Phi]\&^\) - 2\ \(v\&^\)\ \(v\&~\) + \(v\&~\)\^2\))\)}, { \(\[Phi]\&~\) \[Rule] 1\/4\ \(( \(v\&^\)\^2 + 4\ \(\[Phi]\&^\) + \(v\&~\)\ \((\(v\&~\) + \@\(\(v\&^\)\^2 + 16\ \(\[Phi]\&^\) - 2\ \(v\&^\)\ \(v\&~\) + \(v\&~\)\^2\))\) - \(v\&^\)\ \((2\ \(v\&~\) + \@\(\(v\&^\)\^2 + 16\ \(\[Phi]\&^\) - 2\ \(v\&^\)\ \(v\&~\) + \(v\&~\)\^2\))\))\), s \[Rule] 1\/4\ \(( \(v\&^\) + 3\ \(v\&~\) + \@\(\(v\&^\)\^2 + 16\ \(\[Phi]\&^\) - 2\ \(v\&^\)\ \(v\&~\) + \(v\&~\)\^2\))\)}}\)], "Output"] }, Open ]], Cell["It looks like the second expression is the cleanest... ", "Text"], Cell[CellGroupData[{ Cell[BoxData[{ \(D[{\(v\&^\) + \(\@\(\[Phi]\&^\)\) \[Xi]\/\@\(1 + \[Xi]\/2\), \(\[Phi]\&^\) \((1 + \[Xi])\)}, \[Xi]] /. {\[Xi] -> 0}\), \(D[{\(v\&^\) - \(\@\(\[Phi]\&^\)\) \[Xi]\/\@\(1 + \[Xi]\/2\), \(\[Phi]\&^\) \((1 + \[Xi])\)}, \[Xi]] /. {\[Xi] -> 0}\)}], "Input"], Cell[BoxData[ \({\@\(\[Phi]\&^\), \(\[Phi]\&^\)}\)], "Output"], Cell[BoxData[ \({\(-\@\(\[Phi]\&^\)\), \(\[Phi]\&^\)}\)], "Output"] }, Open ]], Cell[BoxData[ \(Solve[ v\_l - \(\@2\) \((\[Phi]\_m - \[Phi]\_l)\)\/\@\(\[Phi]\_l + \[Phi]\_m\) == v\_r + \(\@2\) \((\[Phi]\_m - \[Phi]\_r)\)\/\@\(\[Phi]\_r + \[Phi]\_m\), \[Phi]\_m] // Simplify\)], "Input"] }, Closed]], Cell[CellGroupData[{ Cell["Page 84", "Section"], Cell[BoxData[ \(For\ isothermal\ equations\ the\ eigenvalues\ and\ eigenvectors\ are \)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(Eigensystem[{{0, 1}, {a\^2 - m\^2\/\[Rho]\^2, 2 m\/\[Rho]}}] // Simplify\)], "Input"], Cell[BoxData[ \({{\(-a\) + m\/\[Rho], a + m\/\[Rho]}, {{\[Rho]\/\(m - a\ \[Rho]\), 1}, { \[Rho]\/\(m + a\ \[Rho]\), 1}}}\)], "Output"] }, Open ]], Cell[BoxData[ \(\(Then\ \[Del]\[Lambda]\_1 = \)\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(grad1 = {\[PartialD]\_\[Rho] \[Lambda]\_1, \[PartialD]\_m \[Lambda]\_1} \)], "Input"], Cell[BoxData[ \({\(-\(m\/\[Rho]\^2\)\), 1\/\[Rho]}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(denom1 = grad1 . {\[Rho]\/\(m - a\ \[Rho]\), 1} // Simplify; \n denom1 = denom1 /. {\[Rho] -> \[Rho][\[Xi]], m -> m[\[Xi]]}\)], "Input"], Cell[BoxData[ \(a\/\(\(-m[\[Xi]]\) + a\ \[Rho][\[Xi]]\)\)], "Output"] }, Open ]], Cell[BoxData[ \(The\ first\ differential\ equation\ is\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(\(\[Rho][\[Xi]]\/\(m[\[Xi]] - a\ \[Rho][\[Xi]]\)\)\/denom1 // Simplify \)], "Input"], Cell[BoxData[ \(\(-\(\[Rho][\[Xi]]\/a\)\)\)], "Output"] }, Open ]], Cell[BoxData[ \(same\ as\ equation\ 8.12\)], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(1\/denom1 // Simplify\)], "Input"], Cell[BoxData[ \(\(-\(m[\[Xi]]\/a\)\) + \[Rho][\[Xi]]\)], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Page 127", "Section"], Cell[CellGroupData[{ Cell[BoxData[ \(\(1\/k\) \((Series[u[x, t + k], {k, 0, 3}] - u[x, t])\) + \(1\/h\) u[x, t] \((u[x, t] - Series[u[x - h, t], {h, 0, 4}])\)\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"(", InterpretationBox[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["u", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "+", RowBox[{\(u[x, t]\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}], "+", RowBox[{\(1\/2\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((0, 2)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", "k"}], "+", RowBox[{\(1\/6\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((0, 3)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", \(k\^2\)}], "+", InterpretationBox[\(O[k]\^3\), SeriesData[ k, 0, {}, 0, 3, 1]]}], SeriesData[ k, 0, { Plus[ Derivative[ 0, 1][ u][ x, t], Times[ u[ x, t], Derivative[ 1, 0][ u][ x, t]]], Times[ Rational[ 1, 2], Derivative[ 0, 2][ u][ x, t]], Times[ Rational[ 1, 6], Derivative[ 0, 3][ u][ x, t]]}, 0, 3, 1]], ")"}], "-", RowBox[{\(1\/2\), " ", RowBox[{"(", RowBox[{\(u[x, t]\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}], ")"}], " ", "h"}], "+", RowBox[{\(1\/6\), " ", \(u[x, t]\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((3, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", \(h\^2\)}], "-", RowBox[{\(1\/24\), " ", RowBox[{"(", RowBox[{\(u[x, t]\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((4, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}], ")"}], " ", \(h\^3\)}], "+", InterpretationBox[\(O[h]\^4\), SeriesData[ h, 0, {}, 0, 4, 1]]}], SeriesData[ h, 0, { SeriesData[ k, 0, { Plus[ Derivative[ 0, 1][ u][ x, t], Times[ u[ x, t], Derivative[ 1, 0][ u][ x, t]]], Times[ Rational[ 1, 2], Derivative[ 0, 2][ u][ x, t]], Times[ Rational[ 1, 6], Derivative[ 0, 3][ u][ x, t]]}, 0, 3, 1], Times[ Rational[ -1, 2], u[ x, t], Derivative[ 2, 0][ u][ x, t]], Times[ Rational[ 1, 6], u[ x, t], Derivative[ 3, 0][ u][ x, t]], Times[ Rational[ -1, 24], u[ x, t], Derivative[ 4, 0][ u][ x, t]]}, 0, 4, 1]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(1\/k\) \((Series[u[x, t + k], {k, 0, 3}] - u[x, t])\) + \(1\/\(2\ h\)\) \((Series[f[u[x + h, t]], {h, 0, 3}] - Series[f[u[x - h, t]], {h, 0, 3}])\)\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"(", InterpretationBox[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["u", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "+", RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}], "+", RowBox[{\(1\/2\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((0, 2)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", "k"}], "+", RowBox[{\(1\/6\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((0, 3)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", \(k\^2\)}], "+", InterpretationBox[\(O[k]\^3\), SeriesData[ k, 0, {}, 0, 3, 1]]}], SeriesData[ k, 0, { Plus[ Derivative[ 0, 1][ u][ x, t], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 1, 0][ u][ x, t]]], Times[ Rational[ 1, 2], Derivative[ 0, 2][ u][ x, t]], Times[ Rational[ 1, 6], Derivative[ 0, 3][ u][ x, t]]}, 0, 3, 1]], ")"}], "+", RowBox[{\(1\/2\), " ", RowBox[{"(", RowBox[{ RowBox[{\(1\/3\), " ", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "3"]}], "+", RowBox[{ RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}], "+", RowBox[{\(1\/3\), " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((3, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}], " ", \(h\^2\)}], "+", InterpretationBox[\(O[h]\^3\), SeriesData[ h, 0, {}, 0, 3, 1]]}], SeriesData[ h, 0, { SeriesData[ k, 0, { Plus[ Derivative[ 0, 1][ u][ x, t], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 1, 0][ u][ x, t]]], Times[ Rational[ 1, 2], Derivative[ 0, 2][ u][ x, t]], Times[ Rational[ 1, 6], Derivative[ 0, 3][ u][ x, t]]}, 0, 3, 1], 0, Times[ Rational[ 1, 2], Plus[ Times[ Rational[ 1, 3], Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 3]], Times[ Derivative[ 2][ f][ u[ x, t]], Derivative[ 1, 0][ u][ x, t], Derivative[ 2, 0][ u][ x, t]], Times[ Rational[ 1, 3], Derivative[ 1][ f][ u[ x, t]], Derivative[ 3, 0][ u][ x, t]]]]}, 0, 3, 1]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Series[ \(1\/k\) \((u[x, t + k] - u[x, t])\) + \(1\/\(2\ h\)\) \((f[u[x + h, t]] - f[u[x - h, t]])\) - \(k\/\(2\ h\^2\)\) \((\(f'\)[\(1\/2\) \((u[x, t] + u[x + h, t])\)]\ \((f[u[x + h, t]] - f[u[x, t]])\) - \(f'\)[\(1\/2\) \((u[x, t] + u[x - h, t])\)]\ \((f[u[x, t]] - f[u[x - h, t]])\))\), \n \t\t{k, 0, 2}, {h, 0, 2}] // Simplify\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"(", InterpretationBox[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["u", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "+", RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}], "+", RowBox[{\(1\/6\), " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "3"]}], "+", RowBox[{"3", " ", RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}], "+", RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((3, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}], " ", \(h\^2\)}], "+", InterpretationBox[\(O[h]\^3\), SeriesData[ h, 0, {}, 0, 3, 1]]}], SeriesData[ h, 0, { Plus[ Derivative[ 0, 1][ u][ x, t], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 1, 0][ u][ x, t]]], 0, Times[ Rational[ 1, 6], Plus[ Times[ Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 3]], Times[ 3, Derivative[ 2][ f][ u[ x, t]], Derivative[ 1, 0][ u][ x, t], Derivative[ 2, 0][ u][ x, t]], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 3, 0][ u][ x, t]]]]}, 0, 3, 1]], ")"}], "+", RowBox[{ RowBox[{"(", InterpretationBox[ RowBox[{ RowBox[{\(1\/2\), " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["u", TagBox[\((0, 2)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "-", RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "2"]}], "+", RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}]}]}], ")"}]}], "+", RowBox[{\(1\/48\), " ", RowBox[{"(", RowBox[{ RowBox[{\(-18\), " ", SuperscriptBox[ RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], "2"], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "2"], " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}], "-", RowBox[{ RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"7", " ", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "4"]}], "+", RowBox[{"4", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"3", " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t \), "]"}], "2"]}], "+", RowBox[{"4", " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t \), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((3, 0)\), Derivative], MultilineFunction->None], "[", \(x, t \), "]"}]}]}], ")"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"3", " ", RowBox[{ SuperscriptBox["f", TagBox[\((4)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "4"]}], "+", RowBox[{"21", " ", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "2"], " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}], "+", RowBox[{"2", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((4, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}]}]}], ")"}], " ", \(h\^2\)}], "+", InterpretationBox[\(O[h]\^3\), SeriesData[ h, 0, {}, 0, 3, 1]]}], SeriesData[ h, 0, { Times[ Rational[ 1, 2], Plus[ Derivative[ 0, 2][ u][ x, t], Times[ -1, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 2, Derivative[ 2][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 2]], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 2, 0][ u][ x, t]]]]]], 0, Times[ Rational[ 1, 48], Plus[ Times[ -18, Power[ Derivative[ 2][ f][ u[ x, t]], 2], Power[ Derivative[ 1, 0][ u][ x, t], 2], Derivative[ 2, 0][ u][ x, t]], Times[ -1, Derivative[ 2][ f][ u[ x, t]], Plus[ Times[ 7, Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 4]], Times[ 4, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 3, Power[ Derivative[ 2, 0][ u][ x, t], 2]], Times[ 4, Derivative[ 1, 0][ u][ x, t], Derivative[ 3, 0][ u][ x, t]]]]]], Times[ -1, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 3, Derivative[ 4][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 4]], Times[ 21, Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 2], Derivative[ 2, 0][ u][ x, t]], Times[ 2, Derivative[ 1][ f][ u[ x, t]], Derivative[ 4, 0][ u][ x, t]]]]]]}, 0, 3, 1]], ")"}], " ", "k"}], "+", RowBox[{\(1\/6\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((0, 3)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", \(k\^2\)}], "+", InterpretationBox[\(O[k]\^3\), SeriesData[ k, 0, {}, 0, 3, 1]]}], SeriesData[ k, 0, { SeriesData[ h, 0, { Plus[ Derivative[ 0, 1][ u][ x, t], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 1, 0][ u][ x, t]]], 0, Times[ Rational[ 1, 6], Plus[ Times[ Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 3]], Times[ 3, Derivative[ 2][ f][ u[ x, t]], Derivative[ 1, 0][ u][ x, t], Derivative[ 2, 0][ u][ x, t]], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 3, 0][ u][ x, t]]]]}, 0, 3, 1], SeriesData[ h, 0, { Times[ Rational[ 1, 2], Plus[ Derivative[ 0, 2][ u][ x, t], Times[ -1, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 2, Derivative[ 2][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 2]], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 2, 0][ u][ x, t]]]]]], 0, Times[ Rational[ 1, 48], Plus[ Times[ -18, Power[ Derivative[ 2][ f][ u[ x, t]], 2], Power[ Derivative[ 1, 0][ u][ x, t], 2], Derivative[ 2, 0][ u][ x, t]], Times[ -1, Derivative[ 2][ f][ u[ x, t]], Plus[ Times[ 7, Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 4]], Times[ 4, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 3, Power[ Derivative[ 2, 0][ u][ x, t], 2]], Times[ 4, Derivative[ 1, 0][ u][ x, t], Derivative[ 3, 0][ u][ x, t]]]]]], Times[ -1, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 3, Derivative[ 4][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 4]], Times[ 21, Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 2], Derivative[ 2, 0][ u][ x, t]], Times[ 2, Derivative[ 1][ f][ u[ x, t]], Derivative[ 4, 0][ u][ x, t]]]]]]}, 0, 3, 1], Times[ Rational[ 1, 6], Derivative[ 0, 3][ u][ x, t]]}, 0, 3, 1]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Coefficient[Out[21], h\^1]\)], "Input"], Cell[BoxData[ \(0\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{\(Out[21]\), "/.", RowBox[{"{", RowBox[{ RowBox[{ SuperscriptBox["u", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "->", RowBox[{ RowBox[{"-", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}]}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], "}"}]}], "//", "Simplify"}]], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"(", InterpretationBox[ RowBox[{ RowBox[{\(1\/6\), " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "3"]}], "+", RowBox[{"3", " ", RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}], "+", RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((3, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}], " ", \(h\^2\)}], "+", InterpretationBox[\(O[h]\^3\), SeriesData[ h, 0, {}, 2, 3, 1]]}], SeriesData[ h, 0, { Times[ Rational[ 1, 6], Plus[ Times[ Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 3]], Times[ 3, Derivative[ 2][ f][ u[ x, t]], Derivative[ 1, 0][ u][ x, t], Derivative[ 2, 0][ u][ x, t]], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 3, 0][ u][ x, t]]]]}, 2, 3, 1]], ")"}], "+", RowBox[{ RowBox[{"(", InterpretationBox[ RowBox[{ RowBox[{\(1\/2\), " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["u", TagBox[\((0, 2)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "-", RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"2", " ", RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "2"]}], "+", RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}]}]}], ")"}]}], "+", RowBox[{\(1\/48\), " ", RowBox[{"(", RowBox[{ RowBox[{\(-18\), " ", SuperscriptBox[ RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], "2"], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "2"], " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}], "-", RowBox[{ RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"7", " ", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "4"]}], "+", RowBox[{"4", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"3", " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t \), "]"}], "2"]}], "+", RowBox[{"4", " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t \), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((3, 0)\), Derivative], MultilineFunction->None], "[", \(x, t \), "]"}]}]}], ")"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"3", " ", RowBox[{ SuperscriptBox["f", TagBox[\((4)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "4"]}], "+", RowBox[{"21", " ", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "2"], " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}], "+", RowBox[{"2", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((4, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}]}]}], ")"}], " ", \(h\^2\)}], "+", InterpretationBox[\(O[h]\^3\), SeriesData[ h, 0, {}, 0, 3, 1]]}], SeriesData[ h, 0, { Times[ Rational[ 1, 2], Plus[ Derivative[ 0, 2][ u][ x, t], Times[ -1, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 2, Derivative[ 2][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 2]], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 2, 0][ u][ x, t]]]]]], 0, Times[ Rational[ 1, 48], Plus[ Times[ -18, Power[ Derivative[ 2][ f][ u[ x, t]], 2], Power[ Derivative[ 1, 0][ u][ x, t], 2], Derivative[ 2, 0][ u][ x, t]], Times[ -1, Derivative[ 2][ f][ u[ x, t]], Plus[ Times[ 7, Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 4]], Times[ 4, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 3, Power[ Derivative[ 2, 0][ u][ x, t], 2]], Times[ 4, Derivative[ 1, 0][ u][ x, t], Derivative[ 3, 0][ u][ x, t]]]]]], Times[ -1, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 3, Derivative[ 4][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 4]], Times[ 21, Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 2], Derivative[ 2, 0][ u][ x, t]], Times[ 2, Derivative[ 1][ f][ u[ x, t]], Derivative[ 4, 0][ u][ x, t]]]]]]}, 0, 3, 1]], ")"}], " ", "k"}], "+", RowBox[{\(1\/6\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((0, 3)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", \(k\^2\)}], "+", InterpretationBox[\(O[k]\^3\), SeriesData[ k, 0, {}, 0, 3, 1]]}], SeriesData[ k, 0, { SeriesData[ h, 0, { Times[ Rational[ 1, 6], Plus[ Times[ Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 3]], Times[ 3, Derivative[ 2][ f][ u[ x, t]], Derivative[ 1, 0][ u][ x, t], Derivative[ 2, 0][ u][ x, t]], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 3, 0][ u][ x, t]]]]}, 2, 3, 1], SeriesData[ h, 0, { Times[ Rational[ 1, 2], Plus[ Derivative[ 0, 2][ u][ x, t], Times[ -1, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 2, Derivative[ 2][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 2]], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 2, 0][ u][ x, t]]]]]], 0, Times[ Rational[ 1, 48], Plus[ Times[ -18, Power[ Derivative[ 2][ f][ u[ x, t]], 2], Power[ Derivative[ 1, 0][ u][ x, t], 2], Derivative[ 2, 0][ u][ x, t]], Times[ -1, Derivative[ 2][ f][ u[ x, t]], Plus[ Times[ 7, Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 4]], Times[ 4, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 3, Power[ Derivative[ 2, 0][ u][ x, t], 2]], Times[ 4, Derivative[ 1, 0][ u][ x, t], Derivative[ 3, 0][ u][ x, t]]]]]], Times[ -1, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 3, Derivative[ 4][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 4]], Times[ 21, Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 2], Derivative[ 2, 0][ u][ x, t]], Times[ 2, Derivative[ 1][ f][ u[ x, t]], Derivative[ 4, 0][ u][ x, t]]]]]]}, 0, 3, 1], Times[ Rational[ 1, 6], Derivative[ 0, 3][ u][ x, t]]}, 0, 3, 1]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{\(Out[23]\), "/.", RowBox[{"{", RowBox[{ RowBox[{ SuperscriptBox["u", TagBox[\((0, 2)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "->", RowBox[{ RowBox[{ "2", " ", \(\(\(f'\)'\)[u[x, t]]\), " ", \(\(f'\)[u[x, t]]\), " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "2"]}], "+", RowBox[{\(\(f'\)[u[x, t]]\^2\), RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}]}], "}"}]}], "//", "Simplify"}]], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"(", InterpretationBox[ RowBox[{ RowBox[{\(1\/6\), " ", RowBox[{"(", RowBox[{ RowBox[{ RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "3"]}], "+", RowBox[{"3", " ", RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}], "+", RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((3, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}], " ", \(h\^2\)}], "+", InterpretationBox[\(O[h]\^3\), SeriesData[ h, 0, {}, 2, 3, 1]]}], SeriesData[ h, 0, { Times[ Rational[ 1, 6], Plus[ Times[ Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 3]], Times[ 3, Derivative[ 2][ f][ u[ x, t]], Derivative[ 1, 0][ u][ x, t], Derivative[ 2, 0][ u][ x, t]], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 3, 0][ u][ x, t]]]]}, 2, 3, 1]], ")"}], "+", RowBox[{ RowBox[{"(", InterpretationBox[ RowBox[{ RowBox[{\(1\/48\), " ", RowBox[{"(", RowBox[{ RowBox[{\(-18\), " ", SuperscriptBox[ RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], "2"], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "2"], " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}], "-", RowBox[{ RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"7", " ", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "4"]}], "+", RowBox[{"4", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"3", " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t \), "]"}], "2"]}], "+", RowBox[{"4", " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t \), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((3, 0)\), Derivative], MultilineFunction->None], "[", \(x, t \), "]"}]}]}], ")"}]}]}], ")"}]}], "-", RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{"(", RowBox[{ RowBox[{"3", " ", RowBox[{ SuperscriptBox["f", TagBox[\((4)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "4"]}], "+", RowBox[{"21", " ", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "2"], " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}], "+", RowBox[{"2", " ", RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((4, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}]}]}], ")"}], " ", \(h\^2\)}], "+", InterpretationBox[\(O[h]\^3\), SeriesData[ h, 0, {}, 2, 3, 1]]}], SeriesData[ h, 0, { Times[ Rational[ 1, 48], Plus[ Times[ -18, Power[ Derivative[ 2][ f][ u[ x, t]], 2], Power[ Derivative[ 1, 0][ u][ x, t], 2], Derivative[ 2, 0][ u][ x, t]], Times[ -1, Derivative[ 2][ f][ u[ x, t]], Plus[ Times[ 7, Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 4]], Times[ 4, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 3, Power[ Derivative[ 2, 0][ u][ x, t], 2]], Times[ 4, Derivative[ 1, 0][ u][ x, t], Derivative[ 3, 0][ u][ x, t]]]]]], Times[ -1, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 3, Derivative[ 4][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 4]], Times[ 21, Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 2], Derivative[ 2, 0][ u][ x, t]], Times[ 2, Derivative[ 1][ f][ u[ x, t]], Derivative[ 4, 0][ u][ x, t]]]]]]}, 2, 3, 1]], ")"}], " ", "k"}], "+", RowBox[{\(1\/6\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((0, 3)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", \(k\^2\)}], "+", InterpretationBox[\(O[k]\^3\), SeriesData[ k, 0, {}, 0, 3, 1]]}], SeriesData[ k, 0, { SeriesData[ h, 0, { Times[ Rational[ 1, 6], Plus[ Times[ Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 3]], Times[ 3, Derivative[ 2][ f][ u[ x, t]], Derivative[ 1, 0][ u][ x, t], Derivative[ 2, 0][ u][ x, t]], Times[ Derivative[ 1][ f][ u[ x, t]], Derivative[ 3, 0][ u][ x, t]]]]}, 2, 3, 1], SeriesData[ h, 0, { Times[ Rational[ 1, 48], Plus[ Times[ -18, Power[ Derivative[ 2][ f][ u[ x, t]], 2], Power[ Derivative[ 1, 0][ u][ x, t], 2], Derivative[ 2, 0][ u][ x, t]], Times[ -1, Derivative[ 2][ f][ u[ x, t]], Plus[ Times[ 7, Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 4]], Times[ 4, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 3, Power[ Derivative[ 2, 0][ u][ x, t], 2]], Times[ 4, Derivative[ 1, 0][ u][ x, t], Derivative[ 3, 0][ u][ x, t]]]]]], Times[ -1, Derivative[ 1][ f][ u[ x, t]], Plus[ Times[ 3, Derivative[ 4][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 4]], Times[ 21, Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 2], Derivative[ 2, 0][ u][ x, t]], Times[ 2, Derivative[ 1][ f][ u[ x, t]], Derivative[ 4, 0][ u][ x, t]]]]]]}, 2, 3, 1], Times[ Rational[ 1, 6], Derivative[ 0, 3][ u][ x, t]]}, 0, 3, 1]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(Series[\(f'\)[\(1\/2\) \((u[x, t] + u[x - h, t])\)], {h, 0, 2}]\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{ SuperscriptBox["f", "\[Prime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], "-", RowBox[{\(1\/2\), " ", RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}], ")"}], " ", "h"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{\(1\/8\), " ", RowBox[{ SuperscriptBox["f", TagBox[\((3)\), Derivative], MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", SuperscriptBox[ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "2"]}], "+", RowBox[{\(1\/4\), " ", RowBox[{ SuperscriptBox["f", "\[DoublePrime]", MultilineFunction->None], "[", \(u[x, t]\), "]"}], " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}], " ", \(h\^2\)}], "+", InterpretationBox[\(O[h]\^3\), SeriesData[ h, 0, {}, 0, 3, 1]]}], SeriesData[ h, 0, { Derivative[ 1][ f][ u[ x, t]], Times[ Rational[ -1, 2], Derivative[ 2][ f][ u[ x, t]], Derivative[ 1, 0][ u][ x, t]], Plus[ Times[ Rational[ 1, 8], Derivative[ 3][ f][ u[ x, t]], Power[ Derivative[ 1, 0][ u][ x, t], 2]], Times[ Rational[ 1, 4], Derivative[ 2][ f][ u[ x, t]], Derivative[ 2, 0][ u][ x, t]]]}, 0, 3, 1]]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(u[x, t] + Series[u[x + h, t], {h, 0, 3}]\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{\(2\ u[x, t]\), "+", RowBox[{ RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", "h"}], "+", RowBox[{\(1\/2\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", \(h\^2\)}], "+", RowBox[{\(1\/6\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((3, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", \(h\^3\)}], "+", InterpretationBox[\(O[h]\^4\), SeriesData[ h, 0, {}, 0, 4, 1]]}], SeriesData[ h, 0, { Times[ 2, u[ x, t]], Derivative[ 1, 0][ u][ x, t], Times[ Rational[ 1, 2], Derivative[ 2, 0][ u][ x, t]], Times[ Rational[ 1, 6], Derivative[ 3, 0][ u][ x, t]]}, 0, 4, 1]]], "Output"] }, Open ]], Cell[BoxData[""], "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(Factor[2 - 3 x + x\^2]\)], "Input"], Cell[BoxData[ \(\((\(-2\) + x)\)\ \((\(-1\) + x)\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(Plot[%, {x, 0, 4}]; \)\)], "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 /Courier findfont 10 scalefont setfont % Scaling calculations 0.0238095 0.238095 0.0382592 0.0941766 [ [.2619 .02576 -3 -9 ] [.2619 .02576 3 0 ] [.5 .02576 -3 -9 ] [.5 .02576 3 0 ] [.7381 .02576 -3 -9 ] [.7381 .02576 3 0 ] [.97619 .02576 -3 -9 ] [.97619 .02576 3 0 ] [.01131 .13244 -6 -4.5 ] [.01131 .13244 0 4.5 ] [.01131 .22661 -6 -4.5 ] [.01131 .22661 0 4.5 ] [.01131 .32079 -6 -4.5 ] [.01131 .32079 0 4.5 ] [.01131 .41497 -6 -4.5 ] [.01131 .41497 0 4.5 ] [.01131 .50914 -6 -4.5 ] [.01131 .50914 0 4.5 ] [.01131 .60332 -6 -4.5 ] [.01131 .60332 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 .03826 m .2619 .04451 L s [(1)] .2619 .02576 0 1 Mshowa .5 .03826 m .5 .04451 L s [(2)] .5 .02576 0 1 Mshowa .7381 .03826 m .7381 .04451 L s [(3)] .7381 .02576 0 1 Mshowa .97619 .03826 m .97619 .04451 L s [(4)] .97619 .02576 0 1 Mshowa .125 Mabswid .07143 .03826 m .07143 .04201 L s .11905 .03826 m .11905 .04201 L s .16667 .03826 m .16667 .04201 L s .21429 .03826 m .21429 .04201 L s .30952 .03826 m .30952 .04201 L s .35714 .03826 m .35714 .04201 L s .40476 .03826 m .40476 .04201 L s .45238 .03826 m .45238 .04201 L s .54762 .03826 m .54762 .04201 L s .59524 .03826 m .59524 .04201 L s .64286 .03826 m .64286 .04201 L s .69048 .03826 m .69048 .04201 L s .78571 .03826 m .78571 .04201 L s .83333 .03826 m .83333 .04201 L s .88095 .03826 m .88095 .04201 L s .92857 .03826 m .92857 .04201 L s .25 Mabswid 0 .03826 m 1 .03826 L s .02381 .13244 m .03006 .13244 L s [(1)] .01131 .13244 1 0 Mshowa .02381 .22661 m .03006 .22661 L s [(2)] .01131 .22661 1 0 Mshowa .02381 .32079 m .03006 .32079 L s [(3)] .01131 .32079 1 0 Mshowa .02381 .41497 m .03006 .41497 L s [(4)] .01131 .41497 1 0 Mshowa .02381 .50914 m .03006 .50914 L s [(5)] .01131 .50914 1 0 Mshowa .02381 .60332 m .03006 .60332 L s [(6)] .01131 .60332 1 0 Mshowa .125 Mabswid .02381 .05709 m .02756 .05709 L s .02381 .07593 m .02756 .07593 L s .02381 .09477 m .02756 .09477 L s .02381 .1136 m .02756 .1136 L s .02381 .15127 m .02756 .15127 L s .02381 .17011 m .02756 .17011 L s .02381 .18894 m .02756 .18894 L s .02381 .20778 m .02756 .20778 L s .02381 .24545 m .02756 .24545 L s .02381 .26428 m .02756 .26428 L s .02381 .28312 m .02756 .28312 L s .02381 .30195 m .02756 .30195 L s .02381 .33962 m .02756 .33962 L s .02381 .35846 m .02756 .35846 L s .02381 .3773 m .02756 .3773 L s .02381 .39613 m .02756 .39613 L s .02381 .4338 m .02756 .4338 L s .02381 .45264 m .02756 .45264 L s .02381 .47147 m .02756 .47147 L s .02381 .49031 m .02756 .49031 L s .02381 .52798 m .02756 .52798 L s .02381 .54681 m .02756 .54681 L s .02381 .56565 m .02756 .56565 L s .02381 .58448 m .02756 .58448 L s .02381 .01942 m .02756 .01942 L s .02381 .00059 m .02756 .00059 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 .5 Mabswid .02381 .22661 m .06244 .18325 L .10458 .14161 L .14415 .10787 L .18221 .08033 L .22272 .05631 L .26171 .03834 L .28158 .03112 L .30316 .02477 L .31323 .02233 L .32392 .02012 L .33393 .01839 L .34309 .0171 L .3525 .01606 L .35734 .01564 L .36262 .01527 L .36748 .01502 L .36985 .01492 L .372 .01485 L .37422 .01479 L .37632 .01475 L .37747 .01474 L .37853 .01472 L .37969 .01472 L .38092 .01472 L .38214 .01472 L .38283 .01472 L .38347 .01473 L .38464 .01474 L .38588 .01476 L .38803 .0148 L .3904 .01486 L .39292 .01495 L .39564 .01507 L .40056 .01535 L .40522 .01569 L .41012 .01613 L .41889 .01711 L .42844 .01846 L .43863 .02024 L .457 .02432 L .47835 .03047 L .4985 .03767 L .53696 .05515 L .57788 .07914 L .61729 .1075 L .65517 .13964 L .69551 .1791 L .73434 .22218 L .77561 .27347 L Mistroke .81537 .32824 L .85362 .38587 L .89432 .45253 L .9335 .52192 L .97513 .60123 L .97619 .60332 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]HGAYHf4PAg9QL6QYHgoooojjb`i>GYMcNiEePJ@00076; YF1fS;o;ecQ5D_ooooogjMSAa/c5^YVDR`0009BDU:f][B4Q8Kfm_MkNgP3o07=cLf=SHd92@P3ooooo 0?l00?l0oclo?oogj@00<`00IP00V@00c00c000c<`0cIP0cV@0cc00co`1V001V<`1VIP1VV@1Vc01V o`2I002I<`2IIP2IV@2Ic02Io`3<003<<`3I03>IIIS>IVC>Ic3>I oc?<03?<0Al000<10000H0030F8000<100403@0000l0 0P5P00030@0001`0104U00@170000`40001P00030@00064000<100403@0000l000<10000G@000`40 0@0J00@1;@040A@000D100000@1S00030@0006400P4=00003`000`40001N00815P040CD0104@0005 0@000040I0000`40001Q00030@0000/0000?00030@0005l000<100003`050Cd01@4<00<1H@050F<0 00<100002`0000l000<10000K@040DL0107R00003`000`40001Y00@1C`040Mh0000?00030@0006H0 0`5G0081g00000l000<10000I0020E`00`7I0000103o0J`10@0000l000<100004@000`40000A0003 0@00014000<100004P000`40000>008100<00@004P000`40000A00030@00018000<100004@000`40 000A00<14@000`40000B00030@00014000<100004@000`40000A00030@00018000<100004@000`40 000A00030@00014000<100004P000`40000900003`000`40000A00030@00014000<100004@000`40 000B00030@0000`00P4300030@00014000<100004@000`40000B00030@00014000<100004@000`40 00020@l000<100004P000`40000A00030@00014000<100004@000`40000B00030@00014000<10000 4@000`40000A00030@00018000<100002@0000l000<10000F`020@D000<10000HP000`4000020081 GP000`40001S00030@0000T0000?00030@0005T00P5c0081cP0000l000<10000E`020GL00P7<0000 3`000`40001E0081N`020LX0000?00030@0005<00P5o0081b00000l00P5B0081P`020LH0000?0003 0@0004l00P670081a00000l000<10000C@020H/00P7200003`000`40001<00030@0008d000<10000 _`0000l000<10000BP020I400P6o00003`000`4000180081U@020Kd0000?00030@0004H00P6I0081 ^`0000l000<10000A@000`40002K00030@000;P0000?0081A0020Il00P6h00003`000`4000110081 X`020KH0000?00030@00040000<10000Y@000`40002c00003`000`40000n0081Z@020K<0000?0003 0@0003d000<10000Z`000`40002`00003`000`40000l00030@000:d000<10000[`0000l000<10000 >P020K400P6_00003`000`40000i00030@000;<000<10000[00000l00P4i00030@000;D00P6/0000 3`000`40000f0081^P000`40002Y00003`000`40000e00030@000;/00P6Y00003`000`40000d0003 0@000;h000<10000YP0000l000<10000T000<10000T@000080 00<100002P000`40000M00030@000>/000<10000T00000l000<1000070000`40003]00030@0008l0 000?00030@0001/000<10000k`000`40002>00003`020A/000<10000l@000`40002=00003`000`40 000I00030@000?<000<10000S00000l000<1000060000`40003e00030@0008/0000?00030@0001L0 00<10000m`020H/0000?00030@0001H000<10000nP000`40002800003`000`40000E00030@000?`0 00<10000Q`0000l000<1000050000`40003n00030@0008H0000?00030@0001<000<10000o`010003 0@0008D0000?00814`000`40003o00<000<10000Q00000l000<100004@000`40003o00D000<10000 P`0000l000<1000040000`40003o00L000<10000PP0000l000<100003`000`40003o00T000<10000 P@0000l000<100003P000`40003o00/000<10000P00000l000<100003@000`40003o00d000<10000 O`0000l000<1000030000`40003o00l000<10000OP0000l000<100002`000`40003o014000<10000 O@0000l00P4;00030@000?l04P000`40001m00003`000`40000:00030@000?l04`000`40001l0000 3`000`40000900030@000?l05@000`40001k00003`000`40000800030@000?l05`000`40001j0000 3`000`40000700030@000?l06@000`40001i00003`000`40000600030@000?l06`000`40001h0000 3`000`40000500030@000?l07@000`40001g00003`000`40000500030@000?l07P000`40001f0000 3`020@D000<10000o`0P00030@0007D0000?00030@0000<000<10000o`0R00030@0007@0000?0003 0@00008000<10000o`0S00030@0007@0000?00050@000040o`0W00030@0007<0000?00040@000Ol0 :@000`40001b00001@4:00040@000Ol0:P000`40001a000000<1000030000`400@3o02`000<10000 L0000003004000`00P7o02h000<10000K`00008000<100002P030Ol0;P000`40001^00000`000`40 000900030@000?l0;P000`40001^000010000`40000800030@000?l0;`000`40001]000000D10000 0@0:00030@000?l0<0000`40001/000000D100000@0:00030@000?l0<@000`40001[000000@00@41 2`000`40003o038000<10000JP0000l000<10000o`0b00030@0006X0000?00030@000?l0<`000`40 001Y00003`020Ol0=@000`40001X00003`000`40003o03D000<10000I`0000l000<10000o`0f0003 0@0006H0000?00030@000?l0=`000`40001U00003`000`40003o03L000<10000I@0000l000<10000 o`0h00030@0006@0000?00030@000?l0>@000`40001S00003`000`40003o03X000<10000HP0000l0 0P7o03`000<10000H@0000l000<10000o`0k00030@000640000?00030@000?l0?0000`40001P0000 3`000`40003o03d000<10000G`0000l000<10000o`0n00030@0005h0000?00030@000?l0?P000`40 001N00003`000`40003o03l000<10000G@0000l000<10000o`1000030@0005`0000?0081o`120003 0@0005/0000?00030@000?l0@@000`40001K00003`000`40003o048000<10000FP0000l000<10000 o`1300030@0005T0000?00030@000?l0@`000`40001I00003`000`40003o04@000<10000F00000l0 00<10000o`1500030@0005L0000?00030@000?l0A@000`40001G00003`020Ol0A`000`40001F0000 3`000`40003o04L000<10000E@0000l000<10000o`1800030@0005@0000?00030@000?l0B0000`40 001D00003`000`40003o04T000<10000D`00000400410@/000<10000o`1:00030@00058000001@40 000100X000<10000o`1:00030@000580000400030@0000P000<10000o`1;00030@00054000040003 0@0000P00`7o04`000<10000D000000400410@/000<10000o`1=00030@0004l0000200030@0000X0 00<10000o`1=00030@0004l0000300030@0000T000<10000o`1>00030@0004h0000400030@0000P0 00<10000o`1?00030@0004d000050@X000<10000o`1?00030@0004d0000?00030@000?l0D0000`40 001<00003`000`40003o054000<10000B`0000l00P7o058000<10000B`0000l000<10000o`1B0003 0@0004X0000?00030@000?l0D`000`40001900003`000`40003o05<000<10000B@0000l000<10000 o`1D00030@0004P0000?00030@000?l0E@000`40001700003`000`40003o05D000<10000A`0000l0 00<10000o`1F00030@0004H0000?0081o`1H00030@0004D0000?00030@000?l0E`000`4000150000 3`000`40003o05P000<10000A00000l000<10000o`1I00030@0004<0000?00030@000?l0F@000`40 001300003`000`40003o05X000<10000@P0000l000<10000o`1K00030@000440000?00030@000?l0 F`000`40001100003`020Ol0G@000`40001000003`000`40003o05d000<10000?`0000l000<10000 o`1M00030@0003l0000?00030@000?l0GP000`40000n00003`000`40003o05l000<10000?@0000l0 00<10000o`1O00030@0003d0000?00030@000?l0H0000`40000l00003`000`40003o060000<10000 ?00000l00P7o068000<10000>`0000l000<10000o`1R00030@0003X0000?00030@000?l0HP000`40 000j00003`000`40003o06<000<10000>@0000l000<10000o`1T00030@0003P0000300030@0000T0 00<10000o`1T00030@0003P0000300030@0000T000<10000o`1U00030@0003L000050@X000<10000 o`1U00030@0003L000001040004;00<1o`1V00030@0003H000001001004;00030@000?l0I`000`40 000e000000@00@012`000`40003o06L000<10000=@0000800P4;00030@000?l0J0000`40000d0000 0`000`40000900030@000?l0J@000`40000c00000`000`40000900030@000?l0J@000`40000c0000 3`000`40003o06X000<10000"], ImageRangeCache->{{{0, 431}, {265.875, 0}} -> {-0.148954, -0.469821, 0.0148139, 0.0374522}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(1\/k\) \((Series[u[x, t + k], {k, 0, 3}] - \(1\/2\) \((Series[u[x - h, t], {h, 0, 3}] + Series[u[x + h, t], {h, 0, 3}])\))\) + \(1\/\(2 h\)\) A \((Series[u[x + h, t], {h, 0, 3}] - Series[u[x - h, t], {h, 0, 3}]) \)\)], "Input"], Cell[BoxData[ InterpretationBox[ RowBox[{ RowBox[{"(", InterpretationBox[ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{ SuperscriptBox["u", TagBox[\((0, 1)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], "+", RowBox[{"A", " ", RowBox[{ SuperscriptBox["u", TagBox[\((1, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}], "+", RowBox[{\(1\/2\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((0, 2)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", "k"}], "+", RowBox[{\(1\/6\), " ", RowBox[{ SuperscriptBox["u", TagBox[\((0, 3)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], " ", \(k\^2\)}], "+", InterpretationBox[\(O[k]\^3\), SeriesData[ k, 0, {}, 0, 3, 1]]}], SeriesData[ k, 0, { Plus[ Derivative[ 0, 1][ u][ x, t], Times[ A, Derivative[ 1, 0][ u][ x, t]]], Times[ Rational[ 1, 2], Derivative[ 0, 2][ u][ x, t]], Times[ Rational[ 1, 6], Derivative[ 0, 3][ u][ x, t]]}, 0, 3, 1]], ")"}], "+", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"-", FractionBox[ RowBox[{ SuperscriptBox["u", TagBox[\((2, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}], \(2\ k\)]}], "+", RowBox[{\(1\/6\), " ", "A", " ", RowBox[{ SuperscriptBox["u", TagBox[\((3, 0)\), Derivative], MultilineFunction->None], "[", \(x, t\), "]"}]}]}], ")"}], " ", \(h\^2\)}], "+", InterpretationBox[\(O[h]\^3\), SeriesData[ h, 0, {}, 0, 3, 1]]}], SeriesData[ h, 0, { SeriesData[ k, 0, { Plus[ Derivative[ 0, 1][ u][ x, t], Times[ A, Derivative[ 1, 0][ u][ x, t]]], Times[ Rational[ 1, 2], Derivative[ 0, 2][ u][ x, t]], Times[ Rational[ 1, 6], Derivative[ 0, 3][ u][ x, t]]}, 0, 3, 1], 0, Plus[ Times[ Rational[ -1, 2], Power[ k, -1], Derivative[ 2, 0][ u][ x, t]], Times[ Rational[ 1, 6], A, Derivative[ 3, 0][ u][ x, t]]]}, 0, 3, 1]]], "Output"] }, Open ]] }, Closed]] }, Open ]] }, FrontEndVersion->"X 3.0", ScreenRectangle->{{0, 1024}, {0, 768}}, WindowSize->{616, 623}, WindowMargins->{{182, Automatic}, {9, Automatic}}, PrintingPageRange->{Automatic, Automatic}, PrintingOptions->{"PaperSize"->{612, 792}, "PaperOrientation"->"Portrait", "Magnification"->1} ] (*********************************************************************** 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[1731, 51, 57, 0, 148, "Title"], Cell[CellGroupData[{ Cell[1813, 55, 26, 0, 54, "Section"], Cell[1842, 57, 242, 5, 62, "Input"], Cell[2087, 64, 89, 1, 44, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[2213, 70, 26, 0, 34, "Section"], Cell[CellGroupData[{ Cell[2264, 74, 276, 4, 98, "Input"], Cell[2543, 80, 485, 9, 104, "Output"] }, Open ]], Cell[3043, 92, 457, 9, 140, "Input"], Cell[CellGroupData[{ Cell[3525, 105, 67, 1, 27, "Input"], Cell[3595, 108, 123, 3, 50, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3755, 116, 323, 6, 113, "Input"], Cell[4081, 124, 593, 15, 60, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4711, 144, 131, 2, 33, "Input"], Cell[4845, 148, 456, 11, 113, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5338, 164, 98, 2, 33, "Input"], Cell[5439, 168, 14864, 464, 186, 4648, 333, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[20340, 637, 148, 3, 74, "Input"], Cell[20491, 642, 16889, 514, 279, 5161, 365, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[37417, 1161, 242, 4, 97, "Input"], Cell[37662, 1167, 139, 3, 53, "Message"], Cell[37804, 1172, 280, 6, 111, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[38121, 1183, 191, 4, 68, "Input"], Cell[38315, 1189, 300, 6, 70, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[38652, 1200, 433, 7, 160, "Input"], Cell[39088, 1209, 971, 15, 201, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[40096, 1229, 113, 2, 46, "Input"], Cell[40212, 1233, 140, 2, 48, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[40389, 1240, 65, 1, 27, "Input"], Cell[40457, 1243, 114, 2, 54, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[40608, 1250, 385, 7, 90, "Input"], Cell[40996, 1259, 105, 2, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[41138, 1266, 247, 4, 69, "Input"], Cell[41388, 1272, 56, 1, 27, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[41493, 1279, 26, 0, 34, "Section"], Cell[CellGroupData[{ Cell[41544, 1283, 67, 1, 27, "Input"], Cell[41614, 1286, 123, 3, 50, "Output"] }, Open ]], Cell[41752, 1292, 273, 5, 48, "Input"], Cell[CellGroupData[{ Cell[42050, 1301, 82, 1, 32, "Input"], Cell[42135, 1304, 835, 18, 144, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[43007, 1327, 70, 1, 28, "Input"], Cell[43080, 1330, 593, 15, 60, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[43710, 1350, 75, 1, 32, "Input"], Cell[43788, 1353, 1400, 32, 155, "Output"] }, Open ]], Cell[45203, 1388, 72, 0, 32, "Text"], Cell[CellGroupData[{ Cell[45300, 1392, 298, 4, 137, "Input"], Cell[45601, 1398, 66, 1, 39, "Output"], Cell[45670, 1401, 71, 1, 39, "Output"] }, Open ]], Cell[45756, 1405, 268, 6, 49, "Input"] }, Closed]], Cell[CellGroupData[{ Cell[46061, 1416, 26, 0, 34, "Section"], Cell[46090, 1418, 108, 2, 27, "Input"], Cell[CellGroupData[{ Cell[46223, 1424, 113, 2, 47, "Input"], Cell[46339, 1428, 157, 3, 43, "Output"] }, Open ]], Cell[46511, 1434, 64, 1, 27, "Input"], Cell[CellGroupData[{ Cell[46600, 1439, 113, 2, 27, "Input"], Cell[46716, 1443, 68, 1, 46, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[46821, 1449, 160, 2, 61, "Input"], Cell[46984, 1453, 73, 1, 43, "Output"] }, Open ]], Cell[47072, 1457, 71, 1, 27, "Input"], Cell[CellGroupData[{ Cell[47168, 1462, 110, 2, 52, "Input"], Cell[47281, 1466, 59, 1, 43, "Output"] }, Open ]], Cell[47355, 1470, 57, 1, 27, "Input"], Cell[CellGroupData[{ Cell[47437, 1475, 54, 1, 44, "Input"], Cell[47494, 1478, 70, 1, 43, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[47613, 1485, 27, 0, 34, "Section"], Cell[CellGroupData[{ Cell[47665, 1489, 164, 3, 79, "Input"], Cell[47832, 1494, 3553, 99, 149, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[51422, 1598, 211, 4, 79, "Input"], Cell[51636, 1604, 4796, 127, 168, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[56469, 1736, 470, 9, 138, "Input"], Cell[56942, 1747, 18845, 434, 406, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[75824, 2186, 59, 1, 31, "Input"], Cell[75886, 2189, 35, 1, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[75958, 2195, 762, 21, 31, "Input"], Cell[76723, 2218, 17688, 404, 385, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[94448, 2627, 1023, 27, 70, "Input"], Cell[95474, 2656, 14493, 331, 304, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[110004, 2992, 99, 2, 43, "Input"], Cell[110106, 2996, 2349, 65, 79, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[112492, 3066, 73, 1, 27, "Input"], Cell[112568, 3069, 1160, 35, 43, "Output"] }, Open ]], Cell[113743, 3107, 26, 0, 27, "Input"], Cell[CellGroupData[{ Cell[113794, 3111, 56, 1, 31, "Input"], Cell[113853, 3114, 68, 1, 27, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[113958, 3120, 57, 1, 27, "Input"], Cell[114018, 3123, 14644, 427, 186, 4096, 292, "GraphicsData", "PostScript", "Graphics"] }, Open ]], Cell[CellGroupData[{ Cell[128699, 3555, 332, 8, 117, "Input"], Cell[129034, 3565, 3158, 86, 86, "Output"] }, Open ]] }, Closed]] }, Open ]] } ] *) (*********************************************************************** End of Mathematica Notebook file. ***********************************************************************)