(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.0' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. 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[ 47097, 1568]*) (*NotebookOutlinePosition[ 47762, 1591]*) (* CellTagsIndexPosition[ 47718, 1587]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[TextData[{ StyleBox["Prediction of Conversion \nwith an Axially-Dispersed \nPlug Flow \ Reactor ", FontSize->24], StyleBox["\n", FontWeight->"Plain", FontVariations->{"CompatibilityType"->0}], StyleBox["Author's Data", FontSize->14, FontWeight->"Bold"], StyleBox[": ", FontSize->14], StyleBox["Housam BINOUS\nDepartment of Chemical Engineering\nNational \ Institute of Applied Sciences and Technology\nTunis, TUNISIA\nEmail: \ binoushousam@yahoo.com ", FontSize->14, FontWeight->"Plain"] }], "Title", TextAlignment->Center, FontSize->18, Background->RGBColor[0.68751, 0.996109, 0.796887]], Cell[CellGroupData[{ Cell["\<\ Analytical solution for the conversion of an isothermal axially dispersed PFR \ accomplishing a first order reaction. \ \>", "Subsubtitle", Background->RGBColor[0.717647, 1, 1]], Cell["Clear[u,Di,k,L];", "Input"], Cell[CellGroupData[{ Cell[BoxData[ \(solution = First[DSolve[{0 \[Equal] 1/Pe\ D[c[z], {z, 2}] - D[c[z], {z, 1}] - Da\ c[z]}, c[z], z]] // Simplify\)], "Input"], Cell[BoxData[ \({c[z] \[Rule] \[ExponentialE]\^\(1\/2\ \((Pe\ z - \@Pe\ \@\(4\ Da + \ Pe\)\ z)\)\)\ \((C[ 1] + \[ExponentialE]\^\(\@Pe\ \@\(4\ Da + Pe\)\ z\)\ C[ 2])\)}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["Conc=c[z]/.solution/.C[1]->C1/.C[2]->C2//Factor", "Input"], Cell[BoxData[ \(\[ExponentialE]\^\(\(Pe\ z\)\/2 - 1\/2\ \@Pe\ \@\(4\ Da + Pe\)\ z\)\ \ \((C1 + C2\ \[ExponentialE]\^\(\@Pe\ \@\(4\ Da + Pe\)\ z\))\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(CONCENTRATION = Conc /. \@\(4\ Da + Pe\) -> q\ \@Pe\)], "Input"], Cell[BoxData[ \(\[ExponentialE]\^\(\(Pe\ z\)\/2 - \(Pe\ q\ z\)\/2\)\ \((C1 + C2\ \[ExponentialE]\^\(Pe\ q\ z\))\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Danckwerts boundary conditions for the so-called \"closed\" configuration\ \>", "Subsubtitle", Background->RGBColor[1, 1, 0.65098]], Cell[CellGroupData[{ Cell["\<\ sol=Solve[{(D[%,{z,1}]/.z->1)\[Equal]0,(%/.z->0)-1/Pe (D[%,{z,1}]/.z->0)\ \[Equal]1},{C1,C2}]//Simplify\ \>", "Input"], Cell[BoxData[ \({{C1 \[Rule] \(2\ \[ExponentialE]\^\(Pe\ q\)\ \((1 + \ q)\)\)\/\(\(-\((\(-1\) + q)\)\^2\) + \[ExponentialE]\^\(Pe\ q\)\ \((1 + \ q)\)\^2\), C2 \[Rule] \(2\ \((\(-1\) + q)\)\)\/\(\(-\((\(-1\) + q)\)\^2\) + \ \[ExponentialE]\^\(Pe\ q\)\ \((1 + q)\)\^2\)}}\)], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Conversion", "Subsubtitle", Background->RGBColor[1, 1, 0.65098]], Cell[CellGroupData[{ Cell["CONCENTRATIONEXIT=CONCENTRATION/.z->1/.sol//Simplify", "Input"], Cell[BoxData[ \({\(4\ \[ExponentialE]\^\(1\/2\ Pe\ \((1 + q)\)\)\ q\)\/\(\(-\((\(-1\) + \ q)\)\^2\) + \[ExponentialE]\^\(Pe\ q\)\ \((1 + q)\)\^2\)}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["X=First[1-CONCENTRATIONEXIT]", "Input"], Cell[BoxData[ \(1 - \(4\ \[ExponentialE]\^\(1\/2\ Pe\ \((1 + q)\)\)\ \ q\)\/\(\(-\((\(-1\) + q)\)\^2\) + \[ExponentialE]\^\(Pe\ q\)\ \((1 + \ q)\)\^2\)\)], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Peclet Number", "Subsubtitle", Background->RGBColor[1, 1, 0.65098]], Cell[CellGroupData[{ Cell["Pe=(u L)/Di", "Input"], Cell[BoxData[ \(\(L\ u\)\/Di\)], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Damkohler number", "Subsubtitle", Background->RGBColor[1, 1, 0.65098]], Cell[CellGroupData[{ Cell["Da=k L/u", "Input"], Cell[BoxData[ \(\(k\ L\)\/u\)], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Numerical evaluation", "Subsubtitle", Background->RGBColor[1, 1, 0.65098]], Cell["Di=0.005;k=0.01;L=1;u=0.01;", "Input"], Cell[CellGroupData[{ Cell["Da", "Input"], Cell[BoxData[ \(1.`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["Pe", "Input"], Cell[BoxData[ \(2.`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(q = \@\(1 + 4\ Da/Pe\)\)], "Input"], Cell[BoxData[ \(1.7320508075688772`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["X", "Input"], Cell[BoxData[ \(0.5526014772421974`\)], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Plot of concentration profile in the PFR", "Subsubtitle", Background->RGBColor[1, 1, 0.65098]], Cell[CellGroupData[{ Cell["\<\ plt1=Plot[CONCENTRATION/.sol,{z,0,1},PlotStyle->RGBColor[1,0,0]]\ \>", "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.952381 -0.885734 2.01263 [ [.21429 .10808 -9 -9 ] [.21429 .10808 9 0 ] [.40476 .10808 -9 -9 ] [.40476 .10808 9 0 ] [.59524 .10808 -9 -9 ] [.59524 .10808 9 0 ] [.78571 .10808 -9 -9 ] [.78571 .10808 9 0 ] [.97619 .10808 -3 -9 ] [.97619 .10808 3 0 ] [.01131 .01995 -24 -4.5 ] [.01131 .01995 0 4.5 ] [.01131 .22121 -24 -4.5 ] [.01131 .22121 0 4.5 ] [.01131 .32185 -18 -4.5 ] [.01131 .32185 0 4.5 ] [.01131 .42248 -24 -4.5 ] [.01131 .42248 0 4.5 ] [.01131 .52311 -18 -4.5 ] [.01131 .52311 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 .21429 .12058 m .21429 .12683 L s [(0.2)] .21429 .10808 0 1 Mshowa .40476 .12058 m .40476 .12683 L s [(0.4)] .40476 .10808 0 1 Mshowa .59524 .12058 m .59524 .12683 L s [(0.6)] .59524 .10808 0 1 Mshowa .78571 .12058 m .78571 .12683 L s [(0.8)] .78571 .10808 0 1 Mshowa .97619 .12058 m .97619 .12683 L s [(1)] .97619 .10808 0 1 Mshowa .125 Mabswid .07143 .12058 m .07143 .12433 L s .11905 .12058 m .11905 .12433 L s .16667 .12058 m .16667 .12433 L s .2619 .12058 m .2619 .12433 L s .30952 .12058 m .30952 .12433 L s .35714 .12058 m .35714 .12433 L s .45238 .12058 m .45238 .12433 L s .5 .12058 m .5 .12433 L s .54762 .12058 m .54762 .12433 L s .64286 .12058 m .64286 .12433 L s .69048 .12058 m .69048 .12433 L s .7381 .12058 m .7381 .12433 L s .83333 .12058 m .83333 .12433 L s .88095 .12058 m .88095 .12433 L s .92857 .12058 m .92857 .12433 L s .25 Mabswid 0 .12058 m 1 .12058 L s .02381 .01995 m .03006 .01995 L s [(0.45)] .01131 .01995 1 0 Mshowa .02381 .22121 m .03006 .22121 L s [(0.55)] .01131 .22121 1 0 Mshowa .02381 .32185 m .03006 .32185 L s [(0.6)] .01131 .32185 1 0 Mshowa .02381 .42248 m .03006 .42248 L s [(0.65)] .01131 .42248 1 0 Mshowa .02381 .52311 m .03006 .52311 L s [(0.7)] .01131 .52311 1 0 Mshowa .125 Mabswid .02381 .04008 m .02756 .04008 L s .02381 .0602 m .02756 .0602 L s .02381 .08033 m .02756 .08033 L s .02381 .10046 m .02756 .10046 L s .02381 .14071 m .02756 .14071 L s .02381 .16084 m .02756 .16084 L s .02381 .18096 m .02756 .18096 L s .02381 .20109 m .02756 .20109 L s .02381 .24134 m .02756 .24134 L s .02381 .26147 m .02756 .26147 L s .02381 .28159 m .02756 .28159 L s .02381 .30172 m .02756 .30172 L s .02381 .34197 m .02756 .34197 L s .02381 .3621 m .02756 .3621 L s .02381 .38222 m .02756 .38222 L s .02381 .40235 m .02756 .40235 L s .02381 .4426 m .02756 .4426 L s .02381 .46273 m .02756 .46273 L s .02381 .48286 m .02756 .48286 L s .02381 .50298 m .02756 .50298 L s .02381 .54324 m .02756 .54324 L s .02381 .56336 m .02756 .56336 L s .02381 .58349 m .02756 .58349 L s .02381 .60361 m .02756 .60361 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 1 0 0 r .5 Mabswid .02381 .60332 m .06244 .56156 L .10458 .51766 L .14415 .47796 L .18221 .44117 L .22272 .40349 L .26171 .36866 L .30316 .33319 L .34309 .30053 L .3815 .27053 L .42237 .24014 L .46172 .21239 L .49955 .18715 L .53984 .16185 L .57861 .13909 L .61984 .11666 L .65954 .09685 L .69774 .07952 L .73838 .06305 L .77751 .04922 L .81909 .03683 L .85916 .02729 L .89771 .02051 L .91765 .01798 L .92782 .01696 L .93871 .01607 L .94885 .01544 L .9581 .01504 L .96255 .0149 L .96489 .01484 L .96739 .01479 L .96956 .01476 L .97072 .01474 L .97196 .01473 L .9731 .01472 L .97414 .01472 L .97514 .01472 L .97619 .01472 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]HGAYHf4PAg9QL6QYHgeoo`il00Moo`0000Eoo`00Oomoo`0000Qoo`D000Aoo`03001oogoo009oo`800>Ao o`Yl01Eoo`0000Eoo`00Oomoo`0000Qoo`04001oogoo0003Ool20005Ool00`00Oomoo`3NOol5O00O Ool00005Ool007ooOol00009Ool00`00Ool00003Ool00`00Oomoo`04Ool00`00Oomoo`3JOol4O00T Ool0009oo`8000]oo`8000=oo`<000Aoo`03001oogoo0=Ioo`Al02Qoo`006goo00<007ooOol0dWoo 17`0;7oo000KOol2003@Ool3O00`Ool001]oo`03001oogoo0<]oo`Al03=oo`006goo00<007ooOol0 agoo17`0=goo000KOol00`00Oomoo`34Ool3O00kOol001]oo`03001oogoo0<9oo`9l03ioo`006goo 00<007ooOol0_goo0g`0@7oo000KOol2002nOol2O013Ool001]oo`03001oogoo0;Yoo`=l04Eoo`00 6goo00<007ooOol0]goo0g`0B7oo000KOol00`00Oomoo`2eOol2O01;Ool001]oo`03001oogoo0;9o o`=l04eoo`006goo0P00:Goo0P0017oo0P000goo100097oo0P0017oo0P001Goo0`008goo0P0017oo 0P000goo0`0057oo0W`03Woo0P0017oo0P0017oo0P00:Woo1@001Goo000KOol00`00Oomoo`0WOol0 1000Oomoo`0027oo00<007ooOol097oo00@007ooOol000]oo`03001oogoo025oo`04001oogoo0008 Ool01000Oomoo`004Goo0W`03goo00@007ooOol000Qoo`04001oogoo000[Ool00`00Oomoo`05Ool0 01]oo`03001oogoo02Moo`04001oogoo0009Ool00`00Oomoo`0SOol01000Oomoo`0027oo1@008Woo 00@007ooOol000Qoo`04001oogoo000?Ool2O00AOol01000Oomoo`0027oo00@007ooOol002]oo`03 001oogoo00Eoo`006goo00<007ooOol09goo00@007ooOol000Yoo`03001oogoo029oo`04001oogoo 0008Ool01000Oomoo`008goo00@007ooOol000Qoo`<000ioo`9l01=oo`04001oogoo0009Ool2000/ Ool00`00Oomoo`05Ool001]oo`03001oogoo02Moo`04001oogoo0008Ool01000Oomoo`008goo00@0 07ooOol000Uoo`03001oo`0002=oo`04001oogoo0009Ool00`00Oomoo`0;Ool2O00EOol01000Oomo o`0027oo00@007ooOol002]oo`03001oogoo00Eoo`006goo00<007ooOol0:7oo0P002Woo0P009Goo 0P002goo0P0097oo0P002Woo0`0027oo0g`067oo0P002Woo0P00:goo0P001goo000KOol2002TOol2 O01MOol001]oo`03001oogoo0:5oo`9l05moo`006goo00<007ooOol0Wgoo0W`0HGoo000KOol00`00 Oomoo`2MOol2O01SOol001]oo`03001oogoo09]oo`9l06Eoo`0057ooX`000W`0I`00000KOol00`00 Oomoo`09Ool00`00Oomoo`0:Ool00`00Oomoo`0:Ool00`00Oomoo`09Ool00`00Oomoo`0:Ool00`00 Oomoo`0:Ool00`00Oomoo`09Ool00`00Oomoo`0:Ool00`00Oomoo`0:Ool00`00Oomoo`09Ool00`00 Oomoo`0:Ool00`00Oomoo`0:Ool01000Ooml07`027oo00<007ooOol02Woo00<007ooOol02Woo00<0 07ooOol02Goo00<007ooOol02Woo00<007ooOol02Woo00<007ooOol02Goo00<007ooOol02Woo00<0 07ooOol01Goo000KOol00`00Oomoo`0_Ool00`00Oomoo`0`Ool00`00Oomoo`0`Ool2O00`Ool00`00 Oomoo`0`Ool00`00Oomoo`05Ool001]oo`03001oogoo09=oo`9l06eoo`006goo00<007ooOol0TGoo 0W`0Kgoo000KOol2002AOol00g`0Oomoo`1_Ool001]oo`03001oogoo08ioo`9l079oo`006goo00<0 07ooOol0S7oo0W`0M7oo000KOol00`00Oomoo`2;Ool00g`0Oomoo`1dOol001]oo`03001oogoo08Uo o`9l07Moo`006goo00<007ooOol0Qgoo0W`0NGoo000KOol20027Ool00g`0Oomoo`1iOol001]oo`03 001oogoo08Aoo`9l07aoo`006goo00<007ooOol0Pgoo00=l07ooOol0O7oo000KOol00`00Oomoo`21 Ool2O01oOol001]oo`03001oogoo07moo`9l085oo`006goo0P00Ogoo00=l07ooOol0PGoo000KOol0 0`00Oomoo`1lOol2O024Ool001]oo`03001oogoo07]oo`03O01oogoo08Aoo`006goo00<007ooOol0 NGoo0W`0Qgoo000KOol00`00Oomoo`1gOol2O029Ool001]oo`8007Moo`03O01oogoo08Uoo`006goo 00<007ooOol0M7oo0W`0S7oo000KOol00`00Oomoo`1bOol2O02>Ool001]oo`03001oogoo075oo`03 O01oogoo08ioo`000Woo0P0017oo0P0017oo0P0017oo0P001Goo00<007ooOol0L7oo00=l07ooOol0 Sgoo00001Goo001oogoo000027oo00@007ooOol0009oo`04001oogoo0004Ool00`00Oomoo`1^Ool2 O02BOol00005Ool007ooOol0000;Ool00`00Oomoo`03Ool00`00Oomoo`02Ool2001^Ool00g`0Oomo o`2BOol00005Ool007ooOol00009Ool20004Ool20005Ool00`00Oomoo`1/Ool00g`0Oomoo`2COol0 0005Ool007ooOol00009Ool00`00Oomoo`03Ool00`00Oomoo`04Ool00`00Oomoo`1[Ool00g`0Oomo o`2DOol0009oo`8000Yoo`<000=oo`<000Aoo`03001oogoo06Uoo`9l09Moo`006goo00<007ooOol0 J7oo00=l07ooOol0Ugoo000KOol2001XOol00g`0Oomoo`2HOol001]oo`03001oogoo06Eoo`9l09]o o`006goo00<007ooOol0I7oo00=l07ooOol0Vgoo000KOol00`00Oomoo`1SOol00g`0Oomoo`2LOol0 01]oo`03001oogoo065oo`9l09moo`006goo0P00HGoo00=l07ooOol0Wgoo000KOol00`00Oomoo`1O Ool00g`0Oomoo`2POol001]oo`03001oogoo05eoo`9l0:=oo`006goo00<007ooOol0G7oo00=l07oo Ool0Xgoo000KOol00`00Oomoo`1KOol00g`0Oomoo`2TOol001]oo`03001oogoo05Uoo`9l0:Moo`00 6goo0P00FGoo00=l07ooOol0Ygoo000KOol00`00Oomoo`1GOol00g`0Oomoo`2XOol001]oo`03001o ogoo05Ioo`03O01oogoo0:Uoo`006goo00<007ooOol0E7oo0W`0[7oo000KOol00`00Oomoo`1COol0 0g`0Oomoo`2/Ool001]oo`8005=oo`03O01oogoo0:eoo`006goo00<007ooOol0D7oo0W`0/7oo000K Ool00`00Oomoo`1?Ool00g`0Oomoo`2`Ool000Qoo`8000Aoo`8000=oo`<000Eoo`03001oogoo04io o`03O01oogoo0;5oo`001goo00@007ooOol000Qoo`04001oogoo0004Ool00`00Oomoo`1goo00=l 07ooOol0a7oo000KOol00`00Oomoo`0jOol00g`0Oomoo`35Ool001]oo`03001oogoo03Uoo`03O01o ogoo0Goo00=l07ooOol0agoo000KOol00`00Oomoo`0fOol2O03:Ool001]oo`03 001oogoo03Eoo`03O01oogoo05oo`006goo00<007ooOol07Goo00=l07ooOol0hWoo 000KOol2000MOol00g`0Oomoo`3SOol001]oo`03001oogoo01]oo`03O01oogoo0>Aoo`006goo00<0 07ooOol06Woo00=l07ooOol0iGoo000KOol00`00Oomoo`0IOol00g`0Oomoo`3VOol001]oo`03001o ogoo01Qoo`03O01oogoo0>Moo`006goo00<007ooOol05goo00=l07ooOol0j7oo000KOol2000GOol0 0g`0Oomoo`3YOol001]oo`03001oogoo01Eoo`03O01oogoo0>Yoo`006goo00<007ooOol057oo00=l 07ooOol0jgoo0008Ool20004Ool20005Ool00`00Oomoo`03Ool00`00Oomoo`0COol00g`0Oomoo`3/ Ool000Moo`04001oogoo000:Ool00`00Oomoo`03Ool00`00Oomoo`0BOol00g`0Oomoo`3]Ool000Mo o`04001oogoo000:Ool00`00Oomoo`03Ool2000BOol00g`0Oomoo`3^Ool000Moo`04001oogoo000; Ool00`00Oomoo`02Ool00`00Oomoo`0@Ool00g`0Oomoo`3_Ool000Moo`04001oogoo0008Ool01000 Oomoo`0017oo00<007ooOol03goo00=l07ooOol0l7oo0008Ool20009Ool40004Ool00`00Oomoo`0> Ool00g`0Oomoo`3aOol001]oo`03001oogoo00eoo`03O01oogoo0?9oo`006goo0P003Goo00=l07oo Ool0lgoo000KOol00`00Oomoo`0;Ool00g`0Oomoo`3dOol001]oo`03001oogoo00Yoo`03O01oogoo 0?Eoo`006goo00<007ooOol02Goo00=l07ooOol0mWoo000KOol00`00Oomoo`08Ool00g`0Oomoo`3g Ool001]oo`03001oogoo00Moo`03O01oogoo0?Qoo`006goo0P001goo00=l07ooOol0nGoo000KOol0 0`00Oomoo`05Ool00g`0Oomoo`3jOol001]oo`03001oogoo00Aoo`03O01oogoo0?]oo`006goo00<0 07ooOol00goo00=l07ooOol0o7oo000KOol00`00Oomoo`02Ool00g`0Oomoo`3mOol001]oo`8000=o o`03O01oogoo0?eoo`006goo00D007ooOomoog`00?moo`5oo`006goo00@007ooOoml0?moo`9oo`00 6goo00<007ooO000ogoo0goo000KOol00`00O01oo`3oOol3Ool001]oo`03O00007oo0?moo`=oo`00 6goo00<007ooOol0ogoo0goo000KOol00`00Oomoo`3oOol3Ool001]oo`03001oogoo0?moo`=oo`00 6goo00<007ooOol0ogoo0goo003oOolQOol00?moob5oo`00ogoo8Goo003oOolQOol00?moob5oo`00 ogoo8Goo0000\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.107864, 0.428378, \ 0.00394728, 0.00186786}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["\<\ Numerical solution for the conversion of an isothermal axially dispersed PFR \ accomplishing a second order reaction. The Danckwerts boundary conditions for \ the so-called \"closed\" configuration are used. We use the shooting method \ since it is a split-boundary-value problem. Damkohler number (Da) and Inverse \ Peclet number (1/Pe) must have same order of magnitude unless one gets a \ stiff ODE which cannot be solved easily.\ \>", "Subsubtitle", Background->RGBColor[0.717647, 1, 1]], Cell["Di=0.005;k=0.01;L=1;u=0.01;", "Input"], Cell[CellGroupData[{ Cell["Da", "Input"], Cell[BoxData[ \(1.`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["Pe", "Input"], Cell[BoxData[ \(2.`\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{\(SOL[\[Lambda]_]\), ":=", RowBox[{"NDSolve", "[", RowBox[{\({0 \[Equal] 1/Pe\ D[c[z], {z, 2}] - D[c[z], {z, 1}] - Da\ c[z]^2, \(c'\)[ 0] \[Equal] \[Lambda], \(-1\)/Pe\ \[Lambda] + c[0] \[Equal] 1}\), ",", \(c[z]\), ",", \({z, 0, 1}\), ",", RowBox[{"Method", "->", FormBox[ FormBox["StiffnessSwitching", "TraditionalForm"], "TraditionalForm"]}]}], "]"}]}]], "Input"], Cell[BoxData[ RowBox[{\(General::"spell1"\), \(\(:\)\(\ \)\), "\<\"Possible spelling \ error: new symbol name \\\"\\!\\(SOL\\)\\\" is similar to existing symbol \ \\\"\\!\\(sol\\)\\\". \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \ ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::spell1\\\"]\\)\"\>"}]], "Message"] }, Open ]], Cell["\<\ BCend[\[Lambda]_?NumericQ]:=Flatten[(D[c[z]/.SOL[\[Lambda]],{z,1}])/.z\[Rule]\ 1]\ \>", "Input"], Cell["\<\ \[Xi]=FindRoot[BCend[\[Lambda]]\[Equal]0,{\[Lambda],1,2},MaxIterations\[Rule] \ 1000][[1,2]];\ \>", "Input"], Cell[CellGroupData[{ Cell["\[Xi]", "Input"], Cell[BoxData[ \(\(-0.4127058035716387`\)\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["BCend[\[Xi]]", "Input"], Cell[BoxData[ \({3.0923909266622474`*^-12}\)], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["Conversion", "Subsubtitle", Background->RGBColor[1, 1, 0.65098]], Cell[CellGroupData[{ Cell["1-c[z]/.SOL[\[Xi]]/.z->1", "Input"], Cell[BoxData[ \({0.42774428455308755`}\)], "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Plot of concentration profile in the PFR", "Subsubtitle", Background->RGBColor[1, 1, 0.65098]], Cell[CellGroupData[{ Cell["CONC=First[c[z]/.SOL[\[Xi]]]", "Input"], Cell[BoxData[ RowBox[{ TagBox[\(InterpolatingFunction[{{0.`, 1.`}}, "<>"]\), False, Editable->False], "[", "z", "]"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell["\<\ plt2=Plot[CONC,{z,0,1},PlotStyle->RGBColor[0,0,1]]\ \>", "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.952381 -1.50672 2.65866 [ [.21429 .07598 -9 -9 ] [.21429 .07598 9 0 ] [.40476 .07598 -9 -9 ] [.40476 .07598 9 0 ] [.59524 .07598 -9 -9 ] [.59524 .07598 9 0 ] [.78571 .07598 -9 -9 ] [.78571 .07598 9 0 ] [.97619 .07598 -3 -9 ] [.97619 .07598 3 0 ] [.01131 .22141 -24 -4.5 ] [.01131 .22141 0 4.5 ] [.01131 .35434 -18 -4.5 ] [.01131 .35434 0 4.5 ] [.01131 .48728 -24 -4.5 ] [.01131 .48728 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 .21429 .08848 m .21429 .09473 L s [(0.2)] .21429 .07598 0 1 Mshowa .40476 .08848 m .40476 .09473 L s [(0.4)] .40476 .07598 0 1 Mshowa .59524 .08848 m .59524 .09473 L s [(0.6)] .59524 .07598 0 1 Mshowa .78571 .08848 m .78571 .09473 L s [(0.8)] .78571 .07598 0 1 Mshowa .97619 .08848 m .97619 .09473 L s [(1)] .97619 .07598 0 1 Mshowa .125 Mabswid .07143 .08848 m .07143 .09223 L s .11905 .08848 m .11905 .09223 L s .16667 .08848 m .16667 .09223 L s .2619 .08848 m .2619 .09223 L s .30952 .08848 m .30952 .09223 L s .35714 .08848 m .35714 .09223 L s .45238 .08848 m .45238 .09223 L s .5 .08848 m .5 .09223 L s .54762 .08848 m .54762 .09223 L s .64286 .08848 m .64286 .09223 L s .69048 .08848 m .69048 .09223 L s .7381 .08848 m .7381 .09223 L s .83333 .08848 m .83333 .09223 L s .88095 .08848 m .88095 .09223 L s .92857 .08848 m .92857 .09223 L s .25 Mabswid 0 .08848 m 1 .08848 L s .02381 .22141 m .03006 .22141 L s [(0.65)] .01131 .22141 1 0 Mshowa .02381 .35434 m .03006 .35434 L s [(0.7)] .01131 .35434 1 0 Mshowa .02381 .48728 m .03006 .48728 L s [(0.75)] .01131 .48728 1 0 Mshowa .125 Mabswid .02381 .11506 m .02756 .11506 L s .02381 .14165 m .02756 .14165 L s .02381 .16824 m .02756 .16824 L s .02381 .19482 m .02756 .19482 L s .02381 .248 m .02756 .248 L s .02381 .27458 m .02756 .27458 L s .02381 .30117 m .02756 .30117 L s .02381 .32776 m .02756 .32776 L s .02381 .38093 m .02756 .38093 L s .02381 .40752 m .02756 .40752 L s .02381 .4341 m .02756 .4341 L s .02381 .46069 m .02756 .46069 L s .02381 .51386 m .02756 .51386 L s .02381 .54045 m .02756 .54045 L s .02381 .56704 m .02756 .56704 L s .02381 .59362 m .02756 .59362 L s .02381 .06189 m .02756 .06189 L s .02381 .0353 m .02756 .0353 L s .02381 .00872 m .02756 .00872 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 0 0 1 r .5 Mabswid .02381 .60332 m .06244 .55974 L .10458 .5143 L .14415 .47353 L .18221 .436 L .22272 .3978 L .26171 .3627 L .30316 .32714 L .34309 .29457 L .3815 .26476 L .42237 .2347 L .46172 .20735 L .49955 .18255 L .53984 .15776 L .57861 .13552 L .61984 .11364 L .65954 .09435 L .69774 .07751 L .73838 .06152 L .77751 .04811 L .81909 .03611 L .85916 .02688 L .89771 .02032 L .91765 .01788 L .92782 .01689 L .93871 .01603 L .94885 .01542 L .9581 .01503 L .96255 .01489 L .96489 .01484 L .96739 .01479 L .96956 .01476 L .97072 .01474 L .97196 .01473 L .9731 .01472 L .97414 .01472 L .97514 .01472 L .97619 .01472 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]HGAYHf4PAg9QL6QYHg01l7Ool001]oo`03001oogoo0>=oo`X07aEoo`006goo00<007ooOol0gWoo1@0O7goo000K Ool00`00Oomoo`3IOol501lTOol001]oo`03001oogoo0==oo`H07bUoo`006goo0P00d7oo100O;goo 000KOol00`00Oomoo`3Ool00`0OOomoo`2aOol001]oo`03001oogoo04eoo`0301moogoo0;9oo`006goo 00<007ooOol0C7oo00<07gooOol0/goo000KOol00`00Oomoo`1;Ool00`0OOomoo`2dOol001]oo`03 001oogoo04Uoo`807kMoo`006goo00<007ooOol0B7oo00<07gooOol0]goo000KOol20018Ool00`0O Oomoo`2hOol001]oo`03001oogoo04Ioo`0301moogoo0;Uoo`006goo00<007ooOol0A7oo0P0O_7oo 000KOol00`00Oomoo`13Ool00`0OOomoo`2lOol001]oo`03001oogoo049oo`0301moogoo0;eoo`00 27oo0P0017oo0P001Goo00<007ooOol00goo00<007ooOol0@Goo00<07gooOol0_Woo0007Ool01000 Oomoo`002Woo00<007ooOol00goo00<007ooOol0@7oo00<07gooOol0_goo0007Ool01000Oomoo`00 2Woo00<007ooOol00goo0P00?goo0P0O`Woo0007Ool01000Oomoo`002goo00<007ooOol00Woo00<0 07ooOol0?Goo00<07gooOol0`Woo0007Ool01000Oomoo`0027oo00@007ooOol000Aoo`03001oogoo 03aoo`0301moogoo0<=oo`0027oo0P002Goo100017oo00<007ooOol0>goo00<07gooOol0a7oo000K Ool00`00Oomoo`0jOol00`0OOomoo`35Ool001]oo`03001oogoo03Uoo`0301moogoo07oo00<07gooOol0agoo000KOol2000hOol00`0OOomoo`38Ool001]oo`03001oogoo 03Ioo`0301moogoo0Ool001]oo`80035oo`0301moogoo 09oo`000Woo0P0017oo0P001Goo00<007ooOol00Woo0P00 1Goo00<007ooOol077oo00<07gooOol0hgoo00001Goo001oogoo00002Woo00D007ooOomoo`00009o o`03001oogoo009oo`03001oogoo01]oo`0301moogoo0>Aoo`0000Eoo`00Oomoo`0000Yoo`03001o ogoo00Aoo`03001oogoo009oo`8001]oo`0301moogoo0>Eoo`0000Eoo`00Oomoo`0000]oo`04001o ogooOol20005Ool00`00Oomoo`0IOol00`0OOomoo`3VOol00005Ool007ooOol00008Ool01000Oomo o`000goo00<007ooOol017oo00<007ooOol067oo00<07gooOol0igoo0002Ool20009Ool40003Ool3 0004Ool00`00Oomoo`0GOol00`0OOomoo`3XOol001]oo`03001oogoo01Ioo`0301moogoo0>Uoo`00 6goo00<007ooOol05Goo00<07gooOol0jWoo000KOol00`00Oomoo`0DOol00`0OOomoo`3[Ool001]o o`03001oogoo01=oo`0301moogoo0>aoo`006goo0P004goo00<07gooOol0kGoo000KOol00`00Oomo o`0AOol00`0OOomoo`3^Ool001]oo`03001oogoo011oo`0301moogoo0>moo`006goo00<007ooOol0 3goo00<07gooOol0l7oo000KOol00`00Oomoo`0>Ool00`0OOomoo`3aOol001]oo`03001oogoo00eo o`0301moogoo0?9oo`006goo00<007ooOol037oo00<07gooOol0lgoo000KOol2000=Ool00`0OOomo o`3cOol001]oo`03001oogoo00]oo`0301moogoo0?Aoo`006goo00<007ooOol02Woo00<07gooOol0 mGoo000KOol00`00Oomoo`09Ool00`0OOomoo`3fOol001]oo`03001oogoo00Qoo`0301moogoo0?Mo o`006goo00<007ooOol01goo00<07gooOol0n7oo000KOol00`00Oomoo`06Ool00`0OOomoo`3iOol0 01]oo`8000Ioo`0301moogoo0?Yoo`006goo00<007ooOol017oo00<07gooOol0ngoo000KOol00`00 Oomoo`03Ool00`0OOomoo`3lOol001]oo`03001oogoo009oo`0301moogoo0?eoo`006goo00<007oo Ool00Woo00<07gooOol0oGoo000KOol01@00Oomoogoo01l0ogoo0Goo000KOol01000Oomoo`0Oogoo 0Woo000KOol200000`0OOomoo`3oOol1Ool001]oo`0300007goo0?moo`=oo`006goo00<07gooOol0 ogoo0goo000KOol00`00Oomoo`3oOol3Ool001]oo`03001oogoo0?moo`=oo`006goo00<007ooOol0 ogoo0goo000KOol00`00Oomoo`3oOol3Ool00?moob5oo`00ogoo8Goo003oOolQOol00?moob5oo`00 ogoo8Goo003oOolQOol00001\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.107864, 0.557857, \ 0.00394728, 0.00141399}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Comparing first and second order reactions", "Subsubtitle", Background->RGBColor[1, 1, 0.65098]], Cell[CellGroupData[{ Cell["Show[plt1,plt2]", "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.952381 -0.745838 1.69995 [ [.21429 .09163 -9 -9 ] [.21429 .09163 9 0 ] [.40476 .09163 -9 -9 ] [.40476 .09163 9 0 ] [.59524 .09163 -9 -9 ] [.59524 .09163 9 0 ] [.78571 .09163 -9 -9 ] [.78571 .09163 9 0 ] [.97619 .09163 -3 -9 ] [.97619 .09163 3 0 ] [.01131 .01914 -24 -4.5 ] [.01131 .01914 0 4.5 ] [.01131 .18913 -24 -4.5 ] [.01131 .18913 0 4.5 ] [.01131 .27413 -18 -4.5 ] [.01131 .27413 0 4.5 ] [.01131 .35913 -24 -4.5 ] [.01131 .35913 0 4.5 ] [.01131 .44412 -18 -4.5 ] [.01131 .44412 0 4.5 ] [.01131 .52912 -24 -4.5 ] [.01131 .52912 0 4.5 ] [.01131 .61412 -18 -4.5 ] [.01131 .61412 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 .21429 .10413 m .21429 .11038 L s [(0.2)] .21429 .09163 0 1 Mshowa .40476 .10413 m .40476 .11038 L s [(0.4)] .40476 .09163 0 1 Mshowa .59524 .10413 m .59524 .11038 L s [(0.6)] .59524 .09163 0 1 Mshowa .78571 .10413 m .78571 .11038 L s [(0.8)] .78571 .09163 0 1 Mshowa .97619 .10413 m .97619 .11038 L s [(1)] .97619 .09163 0 1 Mshowa .125 Mabswid .07143 .10413 m .07143 .10788 L s .11905 .10413 m .11905 .10788 L s .16667 .10413 m .16667 .10788 L s .2619 .10413 m .2619 .10788 L s .30952 .10413 m .30952 .10788 L s .35714 .10413 m .35714 .10788 L s .45238 .10413 m .45238 .10788 L s .5 .10413 m .5 .10788 L s .54762 .10413 m .54762 .10788 L s .64286 .10413 m .64286 .10788 L s .69048 .10413 m .69048 .10788 L s .7381 .10413 m .7381 .10788 L s .83333 .10413 m .83333 .10788 L s .88095 .10413 m .88095 .10788 L s .92857 .10413 m .92857 .10788 L s .25 Mabswid 0 .10413 m 1 .10413 L s .02381 .01914 m .03006 .01914 L s [(0.45)] .01131 .01914 1 0 Mshowa .02381 .18913 m .03006 .18913 L s [(0.55)] .01131 .18913 1 0 Mshowa .02381 .27413 m .03006 .27413 L s [(0.6)] .01131 .27413 1 0 Mshowa .02381 .35913 m .03006 .35913 L s [(0.65)] .01131 .35913 1 0 Mshowa .02381 .44412 m .03006 .44412 L s [(0.7)] .01131 .44412 1 0 Mshowa .02381 .52912 m .03006 .52912 L s [(0.75)] .01131 .52912 1 0 Mshowa .02381 .61412 m .03006 .61412 L s [(0.8)] .01131 .61412 1 0 Mshowa .125 Mabswid .02381 .03614 m .02756 .03614 L s .02381 .05314 m .02756 .05314 L s .02381 .07014 m .02756 .07014 L s .02381 .08714 m .02756 .08714 L s .02381 .12113 m .02756 .12113 L s .02381 .13813 m .02756 .13813 L s .02381 .15513 m .02756 .15513 L s .02381 .17213 m .02756 .17213 L s .02381 .20613 m .02756 .20613 L s .02381 .22313 m .02756 .22313 L s .02381 .24013 m .02756 .24013 L s .02381 .25713 m .02756 .25713 L s .02381 .29113 m .02756 .29113 L s .02381 .30813 m .02756 .30813 L s .02381 .32513 m .02756 .32513 L s .02381 .34213 m .02756 .34213 L s .02381 .37613 m .02756 .37613 L s .02381 .39313 m .02756 .39313 L s .02381 .41013 m .02756 .41013 L s .02381 .42712 m .02756 .42712 L s .02381 .46112 m .02756 .46112 L s .02381 .47812 m .02756 .47812 L s .02381 .49512 m .02756 .49512 L s .02381 .51212 m .02756 .51212 L s .02381 .54612 m .02756 .54612 L s .02381 .56312 m .02756 .56312 L s .02381 .58012 m .02756 .58012 L s .02381 .59712 m .02756 .59712 L s .02381 .00214 m .02756 .00214 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 1 0 0 r .5 Mabswid .02381 .51187 m .06244 .4766 L .10458 .43952 L .14415 .40599 L .18221 .37491 L .22272 .34309 L .26171 .31367 L .30316 .28371 L .34309 .25613 L .3815 .23078 L .42237 .20511 L .46172 .18168 L .49955 .16036 L .53984 .13899 L .57861 .11977 L .61984 .10082 L .65954 .08409 L .69774 .06946 L .73838 .05554 L .77751 .04386 L .81909 .03339 L .85916 .02534 L .89771 .01961 L .91765 .01748 L .92782 .01661 L .93871 .01586 L .94885 .01533 L .9581 .01499 L .96255 .01487 L .96489 .01482 L .96739 .01478 L .96956 .01475 L .97072 .01474 L .97196 .01473 L .9731 .01472 L .97414 .01472 L .97514 .01472 L .97619 .01472 L s 0 0 1 r .02381 .60332 m .06244 .57546 L .10458 .5464 L .14415 .52033 L .18221 .49633 L .22272 .47191 L .26171 .44947 L .30316 .42673 L .34309 .4059 L .3815 .38685 L .42237 .36762 L .46172 .35014 L .49955 .33428 L .53984 .31843 L .57861 .30421 L .61984 .29022 L .65954 .27788 L .69774 .26712 L .73838 .2569 L .77751 .24832 L .81909 .24065 L .85916 .23475 L .89771 .23055 L .91765 .22899 L .92782 .22836 L .93871 .22781 L .94885 .22742 L .9581 .22716 L .96255 .22708 L .96489 .22704 L .96739 .22701 L .96956 .22699 L .97072 .22698 L .97196 .22698 L .9731 .22697 L .97414 .22697 L .97514 .22697 L .97619 .22697 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]HGAYHf4PAg9QL6QYHgOol00g`0Oomoo`1N Ool301m@Ool000Moo`04001oogoo0008Ool01000Oomoo`0017oo00<007ooOol0C7oo0W`0GWoo0`0O Dgoo0007Ool01000Oomoo`0027oo00@007ooOol000Aoo`8004aoo`03O01oogoo05Yoo`@07eIoo`00 1goo00@007ooOol000Qoo`<000Eoo`03001oogoo04Yoo`03O01oogoo05Qoo`<07eYoo`001goo00@0 07ooOol000Uoo`03001oogoo00Aoo`03001oogoo04Qoo`9l05Moo`@07eeoo`0027oo0P002Woo0`00 17oo00<007ooOol0Agoo00=l07ooOol0Dgoo100OHGoo000KOol00`00Oomoo`15Ool2O01COol301mU Ool001]oo`8004Eoo`03O01oogoo051oo`<07fQoo`006goo00<007ooOol0@goo00=l07ooOol0Cgoo 0P0OJgoo000KOol00`00Oomoo`11Ool2O01?Ool301m]Ool001]oo`03001oogoo041oo`03O01oogoo 04aoo`<07g1oo`006goo0P00@7oo00=l07ooOol0Bgoo0P0OLgoo000KOol00`00Oomoo`0mOol2O01; Ool301meOol001]oo`03001oogoo03aoo`03O01oogoo04Uoo`807gQoo`006goo00<007ooOol0>Woo 0W`0BGoo0`0ONWoo000KOol00`00Oomoo`0iOol00g`0Oomoo`16Ool301mmOol001]oo`8003Uoo`03 O01oogoo04Eoo`807h1oo`006goo00<007ooOol0=Woo0W`0AGoo0`0OPWoo000KOol00`00Oomoo`0e Ool00g`0Oomoo`12Ool301n5Ool001]oo`03001oogoo03=oo`9l04=oo`807hQoo`006goo0P00goo0P0OUgoo00001Goo001oogoo00002Goo 00<007ooOol00goo00<007ooOol017oo00<007ooOol0:Goo0W`0>goo0`0OVGoo0002Ool2000:Ool3 0003Ool30004Ool00`00Oomoo`0XOol00g`0Oomoo`0iOol201nLOol001]oo`8002Qoo`03O01oogoo 03Qoo`807iioo`006goo00<007ooOol09Woo00=l07ooOol0=goo0P0OX7oo000KOol00`00Oomoo`0T Ool2O00hOol201nROol001]oo`03001oogoo02=oo`03O01oogoo03Ioo`807jAoo`006goo00<007oo Ool08Woo00=l07ooOol0=Goo0P0OYWoo000KOol2000ROol00g`0Oomoo`0dOol201nXOol001]oo`03 001oogoo01moo`9l03Eoo`807jYoo`006goo00<007ooOol07Woo00=l07ooOol09oo`0000Eo o`00Oomoo`0000Yoo`05001oogooOol00002Ool00`00Oomoo`02Ool00`00Oomoo`0KOol201oUOol0 0005Ool007ooOol0000:Ool00`00Oomoo`04Ool00`00Oomoo`02Ool2000JOol201oWOol00005Ool0 07ooOol0000;Ool01000Oomoogoo0P001Goo00<007ooOol05goo0P0OjGoo00001Goo001oogoo0000 27oo00@007ooOol000=oo`03001oogoo00Aoo`03001oogoo01Eoo`807n]oo`000Woo0P002Goo1000 0goo0`0017oo00<007ooOol04goo0P0OkGoo000KOol2000COol00`0OOomoo`3]Ool001]oo`03001o ogoo011oo`807o1oo`006goo00<007ooOol03goo00<07gooOol0l7oo000KOol00`00Oomoo`0>Ool0 0`0OOomoo`3aOol001]oo`03001oogoo00aoo`807oAoo`006goo0P0037oo00<07gooOol0m7oo000K Ool00`00Oomoo`0:Ool00`0OOomoo`3eOol001]oo`03001oogoo00Qoo`807oQoo`006goo00<007oo Ool01goo00<07gooOol0n7oo000KOol20007Ool00`0OOomoo`3iOol001]oo`03001oogoo00Aoo`80 7oaoo`006goo00<007ooOol00goo00<07gooOol0o7oo000KOol00`00Oomoo`02Ool00`0OOomoo`3m Ool001]oo`05001oogooOol07`3oOol1Ool001]oo`8000807omoo`9oo`006goo00<0000OOol0ogoo 0goo0008Ool20004Ool20004Ool20005Ool00`0OOomoo`3oOol3Ool000Moo`04001oogoo0008Ool0 1000Oomoo`0017oo00<007ooOol0ogoo0goo0007Ool01000Oomoo`0027oo00@007ooOol000Aoo`80 0?moo`Aoo`001goo00@007ooOol000Uoo`8000Eoo`03001oogoo0?moo`=oo`001goo00@007ooOol0 00Qoo`04001oogoo003oOol:Ool000Qoo`8000Yoo`800?moo`]oo`00ogoo8Goo003oOolQOol00?mo ob5oo`00ogoo8Goo003oOolQOol00?moob5oo`00\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-0.107864, 0.428704, \ 0.00394728, 0.00221143}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]] }, Closed]] }, Open ]] }, FrontEndVersion->"5.0 for Microsoft Windows", ScreenRectangle->{{0, 1024}, {0, 695}}, CellGrouping->Manual, WindowSize->{495, 537}, WindowMargins->{{0, Automatic}, {Automatic, 0}} ] (******************************************************************* 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[1776, 53, 644, 20, 319, "Title"], Cell[CellGroupData[{ Cell[2445, 77, 189, 4, 66, "Subsubtitle"], Cell[2637, 83, 33, 0, 30, "Input"], Cell[CellGroupData[{ Cell[2695, 87, 180, 4, 110, "Input"], Cell[2878, 93, 220, 4, 45, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3135, 102, 64, 0, 30, "Input"], Cell[3202, 104, 166, 2, 45, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3405, 111, 84, 1, 32, "Input"], Cell[3492, 114, 142, 2, 36, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3671, 121, 143, 3, 66, "Subsubtitle"], Cell[CellGroupData[{ Cell[3839, 128, 128, 3, 48, "Input"], Cell[3970, 133, 296, 5, 70, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[4315, 144, 72, 1, 48, "Subsubtitle"], Cell[CellGroupData[{ Cell[4412, 149, 69, 0, 30, "Input"], Cell[4484, 151, 165, 2, 70, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4686, 158, 45, 0, 30, "Input"], Cell[4734, 160, 169, 3, 70, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[4952, 169, 75, 1, 48, "Subsubtitle"], Cell[CellGroupData[{ Cell[5052, 174, 28, 0, 30, "Input"], Cell[5083, 176, 46, 1, 70, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[5178, 183, 78, 1, 48, "Subsubtitle"], Cell[CellGroupData[{ Cell[5281, 188, 25, 0, 30, "Input"], Cell[5309, 190, 45, 1, 70, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[5403, 197, 82, 1, 48, "Subsubtitle"], Cell[5488, 200, 44, 0, 30, "Input"], Cell[CellGroupData[{ Cell[5557, 204, 19, 0, 30, "Input"], Cell[5579, 206, 37, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5653, 212, 19, 0, 30, "Input"], Cell[5675, 214, 37, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5749, 220, 55, 1, 32, "Input"], Cell[5807, 223, 53, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[5897, 229, 18, 0, 30, "Input"], Cell[5918, 231, 53, 1, 29, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[6020, 238, 102, 1, 48, "Subsubtitle"], Cell[CellGroupData[{ Cell[6147, 243, 89, 2, 48, "Input"], Cell[6239, 247, 12075, 374, 186, 3799, 267, "GraphicsData", "PostScript", \ "Graphics"], Cell[18317, 623, 130, 3, 29, "Output"] }, Open ]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[18508, 633, 503, 8, 188, "Subsubtitle"], Cell[19014, 643, 44, 0, 30, "Input"], Cell[CellGroupData[{ Cell[19083, 647, 19, 0, 30, "Input"], Cell[19105, 649, 37, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[19179, 655, 19, 0, 30, "Input"], Cell[19201, 657, 37, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[19275, 663, 522, 11, 110, "Input"], Cell[19800, 676, 354, 5, 52, "Message"] }, Open ]], Cell[20169, 684, 106, 3, 48, "Input"], Cell[20278, 689, 118, 3, 48, "Input"], Cell[CellGroupData[{ Cell[20421, 696, 22, 0, 30, "Input"], Cell[20446, 698, 58, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[20541, 704, 29, 0, 30, "Input"], Cell[20573, 706, 60, 1, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[20670, 712, 72, 1, 46, "Subsubtitle"], Cell[CellGroupData[{ Cell[20767, 717, 41, 0, 30, "Input"], Cell[20811, 719, 56, 1, 29, "Output"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[20916, 726, 102, 1, 48, "Subsubtitle"], Cell[CellGroupData[{ Cell[21043, 731, 45, 0, 30, "Input"], Cell[21091, 733, 155, 4, 29, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[21283, 742, 75, 2, 30, "Input"], Cell[21361, 746, 11352, 343, 186, 3388, 240, "GraphicsData", "PostScript", \ "Graphics"], Cell[32716, 1091, 130, 3, 29, "Output"] }, Open ]] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[32907, 1101, 104, 1, 48, "Subsubtitle"], Cell[CellGroupData[{ Cell[33036, 1106, 32, 0, 30, "Input"], Cell[33071, 1108, 13853, 450, 186, 4820, 334, "GraphicsData", "PostScript", \ "Graphics"], Cell[46927, 1560, 130, 3, 29, "Output"] }, Open ]] }, Closed]] }, Open ]] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)