Iterated NDSolve. HELP.
- To: mathgroup at smc.vnet.net
- Subject: [mg57047] Iterated NDSolve. HELP.
- From: Chrisantha <ctf20 at sussex.ac.uk>
- Date: Thu, 12 May 2005 22:44:33 -0400 (EDT)
- References: <200505120632.CAA08891@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
> Dear All, I need to detect an event in my system of dif. equations, and then change one of the variables (in a manner which is not specified by the dif. equations), and then re-run the dif equations with the new values. How can I most easily do this? Here is what I have so far. (************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.1' 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 at 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[ 51623, 1349]*) (*NotebookOutlinePosition[ 52259, 1371]*) (* CellTagsIndexPosition[ 52215, 1367]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell[BoxData[{ \(\(data\ = \ Table[0, {i, 100}, {j, 25}];\)\), "\[IndentingNewLine]", \(\(iter\ = \ 0;\)\ \), "\[IndentingNewLine]", \(\({meg, \ temp = 280, melta = \ 280, \ meltb = \ \ 300, \ meltc = \ \ \ 320, \ meltg = \ 340, \[Beta]a = \ 10 + \(80000\ \[ExponentialE]\^temp\)\/\(\[ExponentialE]\^melta + \ \[ExponentialE]\^temp\), \ \[Beta]b = 10 + \(80000\ \[ExponentialE]\^temp\)\/\(\[ExponentialE]\^meltb + \ \[ExponentialE]\^temp\), \ \[Beta]c = 10 + \(80000\ \[ExponentialE]\^temp\)\/\(\[ExponentialE]\^meltc + \ \[ExponentialE]\^temp\), \[Beta]g = 10 + \ \(80000\ \[ExponentialE]\^temp\)\/\(\[ExponentialE]\^meltg + \ \[ExponentialE]\^temp\), lig = 10^2\ 0.6*Exp[\(-25\)/\((8.31*10^\(-3\)*\ temp)\)], \ deg = 0.2*10^17\ * Exp[\(-122\)/\((8.31*10^\(-3\)*\ temp)\)], \ \ \[Alpha] = 200000, \ \ \ d = 0.1\ deg, \ dd = 0.5\ deg\ , \ dg = \ \ deg, \ tmax = 1000000, \ kb = \ lig, \ kc = lig, \ kg = \ lig, k6a\ = \ lig, \ k6b\ = \ lig/2\ , \ k6c\ = \ lig/4\ , \ k6g\ = \ lig/8\ , \ ke\ = \ 0.0001, \ Tsize\ = \ 1, \ k1\ = \ 0.005, \ k2\ = \ 0.001, \ k3\ = \ 0.001, \ k4\ = \ 0.001, \ k5\ = \ 0.001, \ k8\ = \ 0.001, \ k9\ = \ 0.001, \ k10\ = \ 0.001, \ k1d\ = \ 0.00001, \ k2d\ = \ 0.00001, \ k3d\ = \ 0.00001, \ k4d\ = \ 0.00001, \ k5d\ = \ 0.00001, \ \ k9d\ = \ 0.00001\ };\)\[IndentingNewLine]\), "\[IndentingNewLine]", \(Chemoton[data, \ iter, \ meg, \ t0_\ \ , ima_, imb_, imc_, imd_, ime_, iab_, iaab_, ibb_, ibc_, ibbc_, icc_, icg_, iccg_, \ igg_, ia_, ib_, ic_, ig_, ie_, iVdash_, iR_, iTdash_, \ iTstar_, iT_, iS_] := \ Module[{}, \ \[IndentingNewLine]iter\ = \ iter + \ 1; \ \[IndentingNewLine]meg = NDSolve[{\ \[IndentingNewLine]\ \(ma'\)[ t] \[Equal] \ \((2 \((k5\ me[t]\ - \ k5d\ ma[t]\ ma[t])\)\ - \ k1\ ma[t]\ x[t]\ + \ k1d\ mb[t])\), \[IndentingNewLine]\(mb'\)[ t] \[Equal] \ \((k1\ ma[t]\ x[t]\ - \ k1d\ mb[t]\ - \ k2\ mb[t]\ + \ k2d\ mc[t]\ y[t])\)\ , \[IndentingNewLine]\(mc'\)[ t] \[Equal] \ \((k2\ mb[t]\ - \ k2d\ mc[t]\ y[t]\ - \ k3\ mc[t]\ + \ k3d\ md[t]\ Vdash[t])\), \[IndentingNewLine]\(md'\)[ t] \[Equal] \ \((k3\ mc[t]\ - \ k3d\ md[t]\ Vdash[t]\ - \ k4\ md[t]\ + \ k4d\ me[t]\ Tdash[t])\), \[IndentingNewLine]\(me'\)[ t] \[Equal] \ \((k4\ md[t]\ - \ k4d\ me[t]\ Tdash[t]\ - \ k5\ me[t]\ + \ k5d\ ma[t]\ ma[t])\)\ , \[IndentingNewLine]\(a'\)[ t] \[Equal] \ \((k6a\ Vdash[t]\ a[t]\ - \ \[Alpha]\ a[t] b[t] + \ \[Beta]a\ \ ab[t]\ - \ \[Alpha]\ a[t] ab[t] + \ \[Beta]a\ aab[t] + \ d\ \ b[t]/2 - ke\ a[t] e[t])\), \[IndentingNewLine]\(ab'\)[ t] \[Equal] \ \((\[Alpha]\ a[t] b[t] - \[Beta]a\ \ ab[t] - \[Alpha]\ a[t] ab[t] + \ \[Beta]a\ aab[t]\ - \ d\ ab[t])\)\ , \[IndentingNewLine]\(bb'\)[ t] \[Equal] \ \((k6b\ Vdash[t]\ b[t]\ + \ kb\ aab[t] + \ \[Alpha]\ b[t]^2 - \ \[Beta]b\ bb[t] - d\ bb[t]/10\ + \ dd\ cc[t]/20)\)\ , \[IndentingNewLine]\(aab'\)[ t] \[Equal] \ \((\ \[Alpha]\ a[t] ab[t]\ - \ kb\ aab[t] - \ \[Beta]a\ aab[t] - \ d\ aab[t]/10)\)\ , \ \(b'\)[ t] \[Equal] \ \ \ \ \((\ \(-\ \[Alpha]\)\ b[t] a[t] + \ \[Beta]a\ \ ab[t]\ \ - \ \[Alpha]\ b[t] c[t] + \ \[Beta]b\ \ bc[t]\ - \ \[Alpha]\ b[t] bc[t] + \ \[Beta]b\ bbc[t] - 2 \[Alpha]\ \ b[t]^2\ + \ 2\ \[Beta]b\ \ bb[t] - d\ b[t] + \ dd\ c[t]/2 - ke\ \ b[t] e[t])\)\ , \[IndentingNewLine]\(bc'\)[ t] \[Equal] \ \((\[Alpha]\ b[t] c[t] - \[Beta]b\ \ bc[t] - \[Alpha]\ b[t] bc[t] + \ \[Beta]b\ bbc[t]\ - \ dd\ bc[t])\), \[IndentingNewLine]\(bbc'\)[ t] \[Equal] \ \((\ \[Alpha]\ b[t] bc[t]\ - \ kc\ \ bbc[t] - \ \[Beta]b\ bbc[t] - \ dd\ bbc[t]/10)\), \[IndentingNewLine]\(cc'\)[ t] \[Equal] \ \((k6c\ Vdash[t]\ c[t]\ + \ kc\ bbc[t] + \ \[Alpha]\ c[t]^2 - \ \[Beta]c\ cc[t]\ - \ dd\ cc[t]/10\ + \ dg\ gg[t]/20)\), \[IndentingNewLine]\(c'\)[ t] \[Equal] \ \((\ \(-\ \[Alpha]\)\ b[t] c[t] + \ \[Beta]b\ \ bc[t] - 2 \[Alpha]\ \ c[t]^2\ + \ 2\ \[Beta]c\ \ cc[t] - \ dd\ c[t]\ \ - \ \[Alpha]\ c[t] g[t] + \ \[Beta]c\ \ cg[t] - \[Alpha]\ c[t] cg[t] + \ \[Beta]c\ ccg[t]\ + \ dg\ \ g[t]/2 - ke\ c[t] e[t])\)\ , \[IndentingNewLine]\(cg'\)[ t] \[Equal] \((\ \[Alpha]\ c[t] g[t] - \[Beta]c\ \ cg[t] - \[Alpha]\ c[t] cg[t] + \ \[Beta]c\ ccg[t]\ - \ dg\ cg[t])\), \[IndentingNewLine]\(ccg'\)[ t] \[Equal] \ \((\ \[Alpha]\ c[t] cg[t]\ - \ kg\ \ ccg[t] - \ \[Beta]c\ ccg[t]\ - \ dg\ ccg[t]/10)\), \[IndentingNewLine]\(gg'\)[ t] \[Equal] \ \((k6g\ Vdash[t]\ g[t]\ + \ kg\ ccg[t] + \ \[Alpha]\ g[t]^2 - \ \[Beta]g\ gg[t]\ - \ dg\ gg[t]/10)\), \[IndentingNewLine]\(g'\)[ t] \[Equal] \ \((\ \(-\ \[Alpha]\)\ c[t] g[t] + \ \[Beta]c\ \ cg[t] - 2 \[Alpha]\ \ g[t]^2\ + \ 2\ \[Beta]g\ \ gg[t] - \ dg\ g[t] - ke\ g[t] e[t])\)\ , \[IndentingNewLine]\(e'\)[ t]\ \ \ \ \ \[Equal] \((ke\ \ a[t] e[t]\ + \ ke\ \ b[t] e[t] + \ ke\ \ c[t] e[t] + \ ke\ \ g[t]\ e[t] - dg*0\ e[t])\), \[IndentingNewLine]\(Vdash'\)[ t]\ \ \ \ \ \[Equal] \ \((\ k3\ mc[t]\ - \ k3d\ md[t]\ Vdash[t]\ - \ \((k6a\ Vdash[t]\ a[t]\ + \ k6b\ Vdash[t]\ b[t] + \ \ k6c\ Vdash[t]\ c[ t] + \ \ k6g\ Vdash[t]\ g[ t])\))\)\ , \ \[IndentingNewLine]\ \(R'\)[ t]\ \ \ \ \ \[Equal] \ \ \((\((k6a\ Vdash[t]\ a[t]\ + \ k6b\ Vdash[t]\ b[t] + \ \ k6c\ Vdash[t]\ c[ t] + \ \ k6g\ Vdash[t]\ g[t])\)\ + \ k9d\ T[t]\ - \ k9\ Tstar[t]\ R[t])\)\ , \[IndentingNewLine]\ \(Tdash'\)[ t]\ \ \ \ \ \[Equal] \((\ \ k4\ md[t]\ - \ k4d\ me[t]\ Tdash[t]\ - \ k8\ Tdash[t]\ )\)\ , \[IndentingNewLine]\ \(Tstar'\)[ t]\ \ \ \ \ \[Equal] \((\ \ k8\ Tdash[t]\ - \ k9\ Tstar[t]\ R[t]\ + \ k9d\ T[t])\)\ , \[IndentingNewLine]\ \(T'\)[ t]\ \ \ \ \ \[Equal] \ \ \((k9\ Tstar[t]\ R[t]\ - \ k9d\ T[t]\ - \ k10\ T[t]\ S[t])\), \[IndentingNewLine]\ \(S'\)[ t]\ \ \ \[Equal] \ k10\ T[t]\ Tsize\ S[t]\ , \[IndentingNewLine]\(x'\)[ t]\ \ \ \ \ \[Equal] \ 0, \[IndentingNewLine]\ \(y'\)[t]\ \ \ \ \ \[Equal] \ 0, \[IndentingNewLine]\ ma[t0]\ \ \ \ \ \[Equal] \ ima, \ mb[t0]\ \ \ \ \ \[Equal] \ imb, mc[t0]\ \ \ \ \ \[Equal] \ imc, \ md[t0]\ \ \ \ \ \[Equal] \ imd, me[t0]\ \ \ \ \ \[Equal] \ ime, ab[t0] \[Equal] \ iab, \ aab[t0]\ \[Equal] \ iaab, bb[t0] \[Equal] ibb, bc[t0] \[Equal] ibc, bbc[t0] \[Equal] ibbc, \ cc[t0] \[Equal] \ icc, \ cg[t0] \[Equal] icg, \ ccg[t0] \[Equal] iccg, \ gg[t0] \[Equal] \ igg, \ a[t0] \[Equal] \ ia, \ b[t0] \[Equal] \ ib, \ c[t0] \[Equal] \ ic, \ g[t0] \[Equal] \ ig, \ e[t0] \[Equal] \ ie, \ \ \ Vdash[t0]\ \[Equal] \ iVdash, \ R[t0]\ \[Equal] \ iR, \[IndentingNewLine]Tdash[t0]\ \[Equal] \ iTdash, Tstar[t0]\ \[Equal] \ iTstar, T[t0]\ \[Equal] \ iT, S[t0]\ \[Equal] \ iS, x[t0]\ \[Equal] \ 0.01, \[IndentingNewLine]y[t0]\ \[Equal] \ 0.000001}, \ {\ ma, mb, mc, md, me, a, ab, bb, \ aab, \ b, \ bc, \ bbc, \ cc, \ c, \ cg, ccg, \ gg, \ g, e, Vdash, R, Tdash, Tstar, T, S, x, y}, \ {t, \ t0, \ t0 + tmax}, \ MaxSteps \[Rule] Infinity, \ Method\ \[Rule] \ {EventLocator, \ "\<Event\>"\ \[Rule] \ \((S[ t] - 2)\), "\<EventAction\>" \[RuleDelayed] \ {data[\([iter, 1]\)]\ = \ ma[t]; data[\([iter, 2]\)]\ = \ mb[t]; data[\([iter, 3]\)] = \ mc[t]; data[\([iter, 4]\)] = \ md[t]\ ; \ data[\([iter, 5]\)]\ = \ me[t]; \[IndentingNewLine]data[\([iter, 5]\)] = ab[t]; \[IndentingNewLine]data[\([iter, 6]\)] = \ aab[t]; \ \[IndentingNewLine]data[\([iter, 7]\)] = bb[t]; \ \[IndentingNewLine]data[\([iter, 8]\)] = bc[t]; \ \[IndentingNewLine]data[\([iter, 9]\)]\ = \ bbc[t]; \ \[IndentingNewLine]data[\([iter, 10]\)]\ = \ cc[t]; \ \[IndentingNewLine]data[\([iter, 11]\)] cg[t]; \ \[IndentingNewLine]data[\([iter, 12]\)] ccg[t]; \ \[IndentingNewLine]data[\([iter, 13]\)] gg[t]; \ \[IndentingNewLine]data[\([iter, 14]\)] a[t]; \ \[IndentingNewLine]data[\([iter, 15]\)] b[t]; \ \[IndentingNewLine]data[\([iter, 16]\)]\ c[ t]; \ \[IndentingNewLine]data[\([iter, 17]\)] g[t]; \ \[IndentingNewLine]data[\([iter, 18]\)]\ e[ t]; \ \ \ \[IndentingNewLine]data[\([iter, 19]\)]\ = \ Vdash[t]; \ \[IndentingNewLine]data[\([iter, 20]\)] = \ R[t]; \[IndentingNewLine]data[\([iter, 21]\)]\ = Tdash[t]; \[IndentingNewLine]data[\([iter, 22]\)]\ = \ Tstar[t]; \[IndentingNewLine]data[\([iter, 23]\)] = T[t]; \[IndentingNewLine]data[\([iter, 24]\)] = \ S\ [t]; \[IndentingNewLine]data[\([iter, 25]\)]\ = \ t; \ \[IndentingNewLine]Throw[ end\ \ = \ t, \ "\<StopIntegration\>"]}}]\[IndentingNewLine]\ \[IndentingNewLine]]\), "\[IndentingNewLine]", \(Chemoton[data, \ iter, \ meg, \ 0\ \ , 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0.001, 0, 0, 0, 0, 0, 1]\)}], "Input"], Cell[BoxData[ RowBox[{\(General::"spell1"\), \(\(:\)\(\ \)\), "\<\"Possible spelling \ error: new symbol name \\\"\\!\\(meltb\\)\\\" is similar to existing symbol \ \\\"\\!\\(melta\\)\\\". \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \ ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::spell1\\\"]\\)\"\>"}]], "Message"], Cell[BoxData[ RowBox[{\(General::"spell"\), \(\(:\)\(\ \)\), "\<\"Possible spelling \ error: new symbol name \\\"\\!\\(meltc\\)\\\" is similar to existing symbols \ \\!\\({melta, meltb}\\). \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \ ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::spell\\\"]\\)\"\>"}]], "Message"], Cell[BoxData[ RowBox[{\(General::"spell"\), \(\(:\)\(\ \)\), "\<\"Possible spelling \ error: new symbol name \\\"\\!\\(meltg\\)\\\" is similar to existing symbols \ \\!\\({melta, meltb, meltc}\\). \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \ ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::spell\\\"]\\)\"\>"}]], "Message"], Cell[BoxData[ RowBox[{\(General::"spell1"\), \(\(:\)\(\ \)\), "\<\"Possible spelling \ error: new symbol name \\\"\\!\\(\[Beta]b\\)\\\" is similar to existing \ symbol \\\"\\!\\(\[Beta]a\\)\\\". \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \ ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::spell1\\\"]\\)\"\>"}]], "Message"], Cell[BoxData[ RowBox[{\(General::"spell"\), \(\(:\)\(\ \)\), "\<\"Possible spelling \ error: new symbol name \\\"\\!\\(\[Beta]c\\)\\\" is similar to existing \ symbols \\!\\({\[Beta]a, \[Beta]b}\\). \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\ \\\", ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::spell\\\"]\\)\"\>"}]], "Message"], Cell[BoxData[ RowBox[{\(General::"stop"\), \(\(:\)\(\ \)\), "\<\"Further output of \ \\!\\(General :: \\\"spell\\\"\\) will be suppressed during this calculation. \ \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", ButtonStyle->\\\"RefGuideLinkText\ \\\", ButtonFrame->None, ButtonData:>\\\"General::stop\\\"]\\)\"\>"}]], \ "Message"], Cell[BoxData[ RowBox[{\(General::"spell1"\), \(\(:\)\(\ \)\), "\<\"Possible spelling \ error: new symbol name \\\"\\!\\(Vdash\\)\\\" is similar to existing symbol \ \\\"\\!\\(iVdash\\)\\\". \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \ ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::spell1\\\"]\\)\"\>"}]], "Message"], Cell[BoxData[ RowBox[{\(General::"spell"\), \(\(:\)\(\ \)\), "\<\"Possible spelling \ error: new symbol name \\\"\\!\\(Tdash\\)\\\" is similar to existing symbols \ \\!\\({iTdash, Vdash}\\). \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \ ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::spell\\\"]\\)\"\>"}]], "Message"], Cell[BoxData[ RowBox[{\(General::"spell1"\), \(\(:\)\(\ \)\), "\<\"Possible spelling \ error: new symbol name \\\"\\!\\(Tstar\\)\\\" is similar to existing symbol \ \\\"\\!\\(iTstar\\)\\\". \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \ ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::spell1\\\"]\\)\"\>"}]], "Message"], Cell[BoxData[ RowBox[{\(General::"spell1"\), \(\(:\)\(\ \)\), "\<\"Possible spelling \ error: new symbol name \\\"\\!\\(end\\)\\\" is similar to existing symbol \ \\\"\\!\\(End\\)\\\". \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \ ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::spell1\\\"]\\)\"\>"}]], "Message"], Cell[BoxData[ RowBox[{\(General::"stop"\), \(\(:\)\(\ \)\), "\<\"Further output of \ \\!\\(General :: \\\"spell1\\\"\\) will be suppressed during this \ calculation. \\!\\(\\*ButtonBox[\\\"More\[Ellipsis]\\\", \ ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \ ButtonData:>\\\"General::stop\\\"]\\)\"\>"}]], "Message"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{"ma", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"mb", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"mc", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"md", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"me", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"a", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"ab", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"bb", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"aab", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"b", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"bc", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"bbc", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"cc", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"c", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"cg", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"ccg", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"gg", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"g", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"e", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"Vdash", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"R", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"Tdash", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"Tstar", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"T", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"S", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"x", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}], ",", RowBox[{"y", "\[Rule]", TagBox[\(InterpolatingFunction[{{0.`, 140554.44684630845`}}, "<>"]\), False, Editable->False]}]}], "}"}], "}"}]], "Output"] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ \(\(\(Plot[{e[t] /. meg, a[t]\ /. meg, b[t] + 2\ bb[t] /. meg, c[t] + 2\ cc[t] /. meg, \ g[t] + 2\ gg[t] /. meg}\ , \ {t, 0, 140554}, \ PlotStyle \[Rule] {RGBColor[0.5, 0.5, 0.5], RGBColor[1, 0, 1], RGBColor[1, 0, 0], RGBColor[0, 1, 0], RGBColor[0, 0, 1]}, \ PlotRange\ \[Rule] All]\)\(\[IndentingNewLine]\) \)\)], "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 6.77591e-06 0.0147143 0.615682 [ [.15933 .00221 -15 -9 ] [.15933 .00221 15 0 ] [.29485 .00221 -15 -9 ] [.29485 .00221 15 0 ] [.43036 .00221 -15 -9 ] [.43036 .00221 15 0 ] [.56588 .00221 -15 -9 ] [.56588 .00221 15 0 ] [.7014 .00221 -18 -9 ] [.7014 .00221 18 0 ] [.83692 .00221 -18 -9 ] [.83692 .00221 18 0 ] [.97244 .00221 -18 -9 ] [.97244 .00221 18 0 ] [.01131 .13785 -18 -4.5 ] [.01131 .13785 0 4.5 ] [.01131 .26099 -18 -4.5 ] [.01131 .26099 0 4.5 ] [.01131 .38412 -18 -4.5 ] [.01131 .38412 0 4.5 ] [.01131 .50726 -18 -4.5 ] [.01131 .50726 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 .15933 .01471 m .15933 .02096 L s [(20000)] .15933 .00221 0 1 Mshowa .29485 .01471 m .29485 .02096 L s [(40000)] .29485 .00221 0 1 Mshowa .43036 .01471 m .43036 .02096 L s [(60000)] .43036 .00221 0 1 Mshowa .56588 .01471 m .56588 .02096 L s [(80000)] .56588 .00221 0 1 Mshowa .7014 .01471 m .7014 .02096 L s [(100000)] .7014 .00221 0 1 Mshowa .83692 .01471 m .83692 .02096 L s [(120000)] .83692 .00221 0 1 Mshowa .97244 .01471 m .97244 .02096 L s [(140000)] .97244 .00221 0 1 Mshowa .125 Mabswid .05769 .01471 m .05769 .01846 L s .09157 .01471 m .09157 .01846 L s .12545 .01471 m .12545 .01846 L s .19321 .01471 m .19321 .01846 L s .22709 .01471 m .22709 .01846 L s .26097 .01471 m .26097 .01846 L s .32873 .01471 m .32873 .01846 L s .3626 .01471 m .3626 .01846 L s .39648 .01471 m .39648 .01846 L s .46424 .01471 m .46424 .01846 L s .49812 .01471 m .49812 .01846 L s .532 .01471 m .532 .01846 L s .59976 .01471 m .59976 .01846 L s .63364 .01471 m .63364 .01846 L s .66752 .01471 m .66752 .01846 L s .73528 .01471 m .73528 .01846 L s .76916 .01471 m .76916 .01846 L s .80304 .01471 m .80304 .01846 L s .8708 .01471 m .8708 .01846 L s .90468 .01471 m .90468 .01846 L s .93856 .01471 m .93856 .01846 L s .25 Mabswid 0 .01471 m 1 .01471 L s .02381 .13785 m .03006 .13785 L s [(0.2)] .01131 .13785 1 0 Mshowa .02381 .26099 m .03006 .26099 L s [(0.4)] .01131 .26099 1 0 Mshowa .02381 .38412 m .03006 .38412 L s [(0.6)] .01131 .38412 1 0 Mshowa .02381 .50726 m .03006 .50726 L s [(0.8)] .01131 .50726 1 0 Mshowa .125 Mabswid .02381 .0455 m .02756 .0455 L s .02381 .07628 m .02756 .07628 L s .02381 .10707 m .02756 .10707 L s .02381 .16863 m .02756 .16863 L s .02381 .19942 m .02756 .19942 L s .02381 .2302 m .02756 .2302 L s .02381 .29177 m .02756 .29177 L s .02381 .32256 m .02756 .32256 L s .02381 .35334 m .02756 .35334 L s .02381 .41491 m .02756 .41491 L s .02381 .44569 m .02756 .44569 L s .02381 .47648 m .02756 .47648 L s .02381 .53804 m .02756 .53804 L s .02381 .56883 m .02756 .56883 L s .02381 .59961 m .02756 .59961 L s .25 Mabswid .02381 0 m .02381 .61803 L s .5 .5 .5 r .5 Mabswid .02381 .01533 m .06244 .01533 L .10458 .01533 L .14415 .01534 L .18221 .01534 L .22272 .01534 L .26171 .01534 L .30316 .01534 L .34309 .01535 L .3815 .01535 L .40095 .01535 L .42237 .01536 L .46172 .01536 L .48113 .01537 L .502 .01538 L .52314 .01538 L .54323 .0154 L .56303 .01541 L .58142 .01543 L .60194 .01545 L .62055 .01548 L .64142 .01551 L .66062 .01556 L .69917 .01567 L .71799 .01575 L .73866 .01585 L .7598 .01597 L .77908 .01612 L .81799 .01652 L .83707 .0168 L .85783 .0172 L .87761 .01769 L .89862 .01839 L .91809 .01926 L .93637 .02036 L .95655 .02201 L .97505 .0241 L .97619 .02426 L s 1 0 1 r .02381 .01717 m .02846 .01718 L .03279 .01718 L .03754 .01718 L .04262 .01718 L .06244 .01721 L .07379 .01722 L .08451 .01724 L .10458 .01728 L .12455 .01734 L .14299 .0174 L .16331 .01749 L .18234 .01759 L .22018 .01787 L .24043 .01806 L .25895 .01828 L .27954 .01858 L .29866 .01892 L .33685 .01986 L .35738 .02057 L .37598 .02138 L .39523 .02244 L .41605 .02394 L .43599 .02582 L .45705 .02847 L .47669 .03175 L .49502 .03575 L .51538 .04155 L .53393 .04836 L .57378 .06877 L .61211 .09554 L .65289 .12894 L .69216 .16415 L .72991 .20188 L .77011 .2487 L .8088 .30143 L .84994 .36592 L .88956 .43569 L .92767 .50817 L .96823 .58775 L .97619 .60332 L s 1 0 0 r .02381 .0172 m .02499 .0159 L .02605 .01577 L .02846 .01555 L .03068 .01541 L .03279 .01531 L .03507 .01522 L .0375 .01514 L .04016 .01508 L .04262 .01503 L .04496 .01499 L .04749 .01496 L .05205 .01492 L .05483 .0149 L .05748 .01488 L .06244 .01485 L .0675 .01484 L .07289 .01482 L .07798 .01481 L .08259 .0148 L .08803 .01479 L .09391 .01479 L .10458 .01478 L .11408 .01477 L .12452 .01477 L .13541 .01476 L .14558 .01476 L .16594 .01476 L .18449 .01475 L .20358 .01475 L .22434 .01475 L .24411 .01475 L .2541 .01475 L .26512 .01475 L .27454 .01475 L .28447 .01475 L .28966 .01475 L .29535 .01475 L .30051 .01475 L .30532 .01475 L .30795 .01475 L .31037 .01475 L .31164 .01475 L .31303 .01475 L .31373 .01475 L .31448 .01475 L .3152 .01475 L .31585 .01475 L .31711 .01475 L .31828 .01475 L Mistroke .31957 .01475 L .32029 .01475 L .32094 .01475 L .32214 .01475 L .32343 .01475 L .32578 .01475 L .32803 .01475 L .33049 .01475 L .33547 .01475 L .34007 .01475 L .34437 .01475 L .35485 .01475 L .36476 .01475 L .38333 .01475 L .40328 .01475 L .42171 .01475 L .442 .01475 L .46103 .01475 L .49883 .01476 L .51905 .01476 L .53757 .01477 L .5474 .01477 L .55653 .01477 L .56719 .01478 L .57725 .01478 L .58709 .01479 L .59774 .0148 L .60743 .01482 L .61634 .01483 L .62606 .01485 L .63515 .01488 L .64511 .01491 L .6558 .01496 L .67701 .01509 L .69713 .01524 L .70677 .0153 L .71205 .01533 L .71701 .01535 L .72576 .01539 L .73041 .0154 L .73543 .01541 L .74491 .01543 L .75505 .01545 L .76617 .01549 L .77653 .01553 L .79727 .01566 L .81613 .01581 L .85608 .01619 L .89451 .01667 L .9354 .01736 L Mistroke .97477 .01824 L .97619 .01828 L Mfstroke 0 1 0 r .02381 .01613 m .02499 .01521 L .02605 .01519 L .02729 .01516 L .02846 .01513 L .03279 .015 L .03507 .01495 L .0375 .0149 L .04016 .01486 L .04262 .01483 L .04496 .01481 L .04749 .01479 L .05205 .01477 L .05483 .01476 L .05748 .01475 L .06244 .01474 L .0651 .01474 L .0675 .01473 L .07009 .01473 L .07289 .01473 L .07789 .01473 L .08259 .01472 L .08816 .01472 L .09404 .01472 L .09907 .01472 L .10458 .01472 L .10982 .01472 L .11478 .01472 L .12409 .01472 L .13411 .01472 L .145 .01472 L .1544 .01472 L .16478 .01472 L .17574 .01472 L .18579 .01472 L .20513 .01472 L .22599 .01472 L .24491 .01472 L .26561 .01472 L .28522 .01472 L .30617 .01472 L .32752 .01472 L .34767 .01472 L .36769 .01472 L .38614 .01472 L .39583 .01472 L .4049 .01472 L .41485 .01472 L .42044 .01472 L .42554 .01472 L Mistroke .428 .01472 L .43032 .01472 L .43161 .01472 L .43233 .01472 L .43299 .01472 L .43541 .01472 L .43772 .01472 L .43987 .01472 L .44104 .01472 L .44211 .01472 L .44329 .01472 L .44456 .01472 L .44585 .01472 L .4465 .01472 L .44722 .01472 L .44849 .01472 L .44967 .01472 L .45035 .01472 L .45107 .01472 L .45234 .01472 L .45522 .01472 L .45772 .01472 L .46004 .01472 L .4653 .01472 L .47037 .01472 L .47585 .01472 L .4858 .01472 L .49548 .01472 L .5044 .01472 L .5232 .01472 L .53388 .01472 L .54385 .01472 L .55429 .01472 L .56418 .01472 L .58273 .01472 L .59237 .01472 L .60262 .01472 L .61136 .01472 L .62103 .01472 L .63111 .01472 L .64035 .01472 L .6505 .01472 L .65561 .01472 L .6612 .01472 L .66647 .01472 L .67146 .01472 L .68079 .01472 L .68601 .01472 L .69083 .01472 L .70173 .01472 L Mistroke .71141 .01472 L .72047 .01473 L .7411 .01473 L .77989 .01474 L .79962 .01475 L .82114 .01475 L .84174 .01476 L .86086 .01478 L .90001 .0148 L .94162 .01484 L .97619 .01487 L Mfstroke 0 0 1 r .02381 .01647 m .03279 .01827 L .03793 .01928 L .04262 .01994 L .04517 .0202 L .04758 .02041 L .05034 .0206 L .05288 .02074 L .05747 .02092 L .05988 .021 L .06244 .02106 L .06743 .02116 L .07023 .0212 L .0728 .02123 L .07777 .02127 L .0825 .02131 L .08797 .02134 L .09305 .02136 L .09863 .02137 L .10458 .02139 L .11487 .02141 L .12066 .02141 L .12599 .02142 L .13642 .02143 L .14616 .02144 L .16499 .02145 L .18565 .02146 L .22456 .02147 L .26592 .02148 L .30576 .0215 L .34409 .02152 L .38488 .02154 L .42414 .02156 L .46189 .02159 L .5021 .02163 L .54078 .02167 L .58192 .0217 L .60245 .02172 L .62155 .02173 L .64115 .02175 L .65966 .02175 L .66947 .02176 L .67501 .02176 L .68011 .02176 L .68914 .02176 L .69414 .02176 L .6987 .02176 L .70102 .02176 L .70354 .02176 L .70593 .02176 L Mistroke .70814 .02176 L .71071 .02177 L .71305 .02177 L .71574 .02177 L .71828 .02177 L .72286 .02178 L .72708 .02179 L .73171 .0218 L .73659 .02182 L .74609 .02187 L .75627 .02194 L .77461 .02212 L .79589 .02239 L .81603 .02272 L .83598 .02315 L .85441 .02369 L .87311 .02442 L .89373 .02548 L .91311 .0268 L .93398 .02865 L .95406 .03096 L .97517 .03409 L .97619 .03426 L Mfstroke 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 177.938}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg<PAVmbKF5d0`40004P0000/B000`400?l00000o`00003ooooooolQ ooooo`00oooooolQooooo`00oooooolQooooo`00oooooolQooooo`00:Oooool5o`0000?ooooo0_l0 0004ooooo`;o00001?ooool2o`0000Cooooo0_l0000:ooooo`?o00000oooool2o`0000Cooooo0_l0 0004ooooo`;o00001?ooool2o`0000Sooooo0ol00004ooooo`;o00001?ooool2o`0000Cooooo0_l0 0004ooooo`;o00002?ooool3o`0000Cooooo0_l00004ooooo`;o00001?ooool2o`0000Cooooo0_l0 0006ooooo`?o00001?ooool2o`0000Cooooo0_l00004ooooo`;o00001?ooool2o`0000Cooooo0_l0 0002ooooo`?o00000_ooool5o`0000?ooooo0_l00004ooooo`;o00001?ooool2o`0000Cooooo0_l0 0003ooooo`?o00001?ooool3o`0000?ooooo0_l00004ooooo`;o00001?ooool2o`0000Cooooo0_l0 0003ooooo`00:_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo 00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00002ooooo`04o`000?oooooooooo o`0000[ooooo00Go0000ooooooooooooooooo`000002ooooo`04o`000?ooooooooooo`0000;ooooo 00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00002ooooo`03o`000?oooooooooo 00Cooooo00Go0000ooooooooooooooooo`000002ooooo`04o`000?ooooooooooo`0000;ooooo00Co 0000oooooooooooo00000_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`00 00Kooooo00Go0000ooooooooooooooooo`000002ooooo`04o`000?ooooooooooo`0000;ooooo00Co 0000oooooooooooo00000_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`00 00Kooooo00?o0000ooooooooool00_ooool01?l0003oooooooooool00002ooooo`04o`000?oooooo ooooo`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00002ooooo`04 o`000?ooooooooooo`0000;ooooo00?o0000ooooooooool00_ooool01?l0003oooooooooool00002 ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003ooooo ooooool00002ooooo`04o`000?ooooooooooo`0000?ooooo00?o0000ooooooooool01?ooool01Ol0 003ooooooooooooooooo000000;ooooo00Co0000oooooooooooo00000_ooool01?l0003ooooooooo ool00002ooooo`04o`000?ooooooooooo`0000;ooooo0Ol00001ooooo`7ooooo000[ooooo`03o`00 0?oooooooooo00;ooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00002oooo o`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00001oooool5o`0000;ooooo00Co 0000oooooooooooo00000_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`00 00;ooooo00Co0000oooooooooooo00001_ooool01Ol0003ooooooooooooooooo000000;ooooo00Co 0000oooooooooooo00000_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`00 00;ooooo00Co0000oooooooooooo00001_ooool01Ol0003ooooooooooooooooo000000;ooooo00Co 0000oooooooooooo00000_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`00 00;ooooo00Co0000oooooooooooo00001_ooool00ol0003oooooooooo`02ooooo`04o`000?oooooo ooooo`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00002ooooo`04 o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool00ol0003oooooooooo`03 ooooo`03o`000?oooooooooo00;ooooo00Co0000oooooooooooo00000_ooool01?l0003ooooooooo ool00002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000oooool01?l0 003oooooooooooooool5o`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003ooooooooo ool00002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool002co oooo00Go0000ooooooooooooooooo`000002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000 oooooooooooo00000_ooool01?l0003oooooooooool00002ooooo`03o`000?oooooooooo00Gooooo 00Co0000oooooooooooo00000oooool01?l0003oooooooooool00002ooooo`04o`000?oooooooooo o`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00006ooooo`Co0000 0oooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooo oooooooo00000_ooool01?l0003oooooooooool00007ooooo`?o00000oooool01?l0003ooooooooo ool00002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0 003oooooooooool00006ooooo`03o`000?oooooooooo00;ooooo00Co0000oooooooooooo00000_oo ool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooo oooo00000_ooool01?l0003oooooooooool00002ooooo`03o`000?oooooooooo00Cooooo00Go0000 ooooooooooooooooo`000002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo 00000_ooool01?l0003oooooooooool00002ooooo`05o`000?ooooooooooooooool000000oooool0 1?l0003oooooooooool00003ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo 00000_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo000Yoooo o`05o`000?ooooooooooooooool000000_ooool01?l0003oooooooooool00002ooooo`04o`000?oo ooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00008oooo o`03o`000?oooooo000000?ooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool0 0002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00001_ooool00ol0003o ooooooooo`04ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool0 1?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000Kooooo00Go0000oooooooooooo ooooo`000002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool0 1?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000Kooooo00?o0000ooooooooool0 0_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooo oooooooo00000_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo 00Go0000ooooooooooooooooo`000003ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooo oooooooo00000_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo 00Go0000ooooooooooooooooo`000004ooooo`03o`000?oooooo000000?ooooo00Co0000oooooooo oooo00000_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo00Co 0000oooooooooooo00000_ooool002Wooooo00Go0000ooooooooooooooooo`000002ooooo`04o`00 0?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00002 ooooo`04o`000?ooooooooooo`0000Wooooo0_l00003ooooo`04o`000?ooooooooooo`0000;ooooo 00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00002ooooo`04o`000?oooooooooo o`0000Oooooo00?o0000ooooooooool00oooool01?l0003oooooooooool00002ooooo`04o`000?oo ooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00006oooo o`05o`000?ooooooooooooooool000000_ooool01?l0003oooooooooool00002ooooo`04o`000?oo ooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00005oooo o`;o00001?ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo00Co 0000oooooooooooo00000_ooool01?l0003oooooooooool00002ooooo`05o`000?ooooooooooo`00 0?ooool00_l00003ooooo`05o`000?ooooooooooooooool000000_ooool01?l0003oooooooooool0 0002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003o ooooooooool00002ooooo`;o00001Oooool2o`0000?ooooo00Co0000oooooooooooo00000_ooool0 1?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo 00000_ooool002[ooooo0ol00004ooooo`;o00001?ooool2o`0000Cooooo0_l00004ooooo`;o0000 2oooool00ol0003oooooooooo`02ooooo`;o00001?ooool2o`0000Cooooo0_l00004ooooo`;o0000 2Oooool2o`0000Cooooo0_l00004ooooo`;o00001?ooool2o`0000Cooooo0_l00008ooooo`?o0000 1?ooool2o`0000Cooooo0_l00004ooooo`;o00001?ooool2o`0000Oooooo00?o0000ooooooooool0 0oooool2o`0000Cooooo0_l00004ooooo`;o00001?ooool2o`0000Cooooo0_l00003ooooo`03o`00 0?oooooooooo00;ooooo0ol00004ooooo`;o00001?ooool2o`0000Cooooo0_l00004ooooo`;o0000 1?ooool00ol0003oooooooooo`04ooooo`03o`000?oooooooooo00;ooooo0_l00004ooooo`;o0000 1?ooool2o`0000Cooooo0_l00003ooooo`00oooooolQooooo`005Oooool00ol0003oooooooooo`3o ooooo`Wooooo000Eooooo`03o`000?oooooooooo0?oooooo2Oooool001Gooooo00?o0000oooooooo ool0ooooool9ooooo`003oooool7o`000?So0?l01ol0000;ooooo`005Oooool3o`00o`03ool00?oo 0?ooo`3o04Soo`3oH?moOgm;ool001;ooooo000Eooooo`03o`000?oooooooooo0=3o003o00?ooooo o`000?ooool04_ooool?ogmoO`03ooooool0003ooooo017ooooo000Eooooo`03o`000?oooooooooo 06;ooooo2_oo0?mToooooaCo003o4?ooool2ogmoOa;ooooo000Eooooo`03o`000?oooooooooo06co oooo1_oo0?mbooooo`[o003o6_ooool001Gooooo00?o0000ooooooooool0L_ooool5ool0ogOooooo 2?l00?lBooooo`005Oooool00ol0003oooooooooo`1gooooo`Goo`3oS?ooool001Gooooo00?o0000 ooooooooool0O?ooool3ool0ohWooooo000Eooooo`;o0000P?ooool3ool0ohKooooo000Eooooo`03 o`000?oooooooooo08;ooooo0_oo0?n4ooooo`005Oooool00ol0003oooooooooo`24ooooo`;oo`3o P_ooool001Gooooo00?o0000ooooooooool0Q_ooool00ooo0?ooooooooooo`1oooooo`005Oooool0 0ol0003oooooooooo`27ooooo`;oo`3oOoooool001Gooooo00?o0000ooooooooool0ROooool2ool0 oggooooo000Eooooo`03o`000?oooooooooo08_ooooo00?oo`3oooooooooool0N_ooool001Gooooo 00?o0000ooooooooool0S?ooool2ool0og[ooooo000Eooooo`;o0000Soooool2ool0ogSooooo000E ooooo`03o`000?oooooooooo093ooooo00?oo`3oooooooooool0MOooool001Gooooo00?o0000oooo ooooool0TOooool2ool0ogGooooo000Eooooo`03o`000?oooooooooo09?ooooo0_oo0?mcooooo`00 5Oooool00ol0003oooooooooo`2Eooooo`03ool0oooooooooooo073ooooo000Eooooo`03o`000?oo oooooooo09Kooooo0_oo0?m`ooooo`005Oooool00ol0003oooooooooo`2Hooooo`03ool0oooooooo oooo06gooooo000Eooooo`03o`000?oooooooooo09Wooooo00?oo`3oooooooooool0K?ooool001Go oooo00?o0000ooooooooool0V_ooool00ooo0?ooooooooooo`1[ooooo`005Oooool2o`0009cooooo 0_oo0?m[ooooo`005Oooool00ol0003oooooooooo`2Mooooo`03ool0oooooooooooo06Sooooo000E ooooo`03o`000?oooooooooo09kooooo00?oo`3oooooooooool0Ioooool001Gooooo00?o0000oooo ooooool0Woooool00ooo0?ooooooooooo`1Vooooo`005Oooool00ol0003oooooooooo`2Pooooo`03 ool0oooooooooooo06Gooooo0002ooooo`;o00001?ooool01?l0003oooooooooooooool5o`0000Co oooo00?o0000ooooooooool0XOooool2ool0ofGooooo00001Ooooooo0000oooooooooooo000000So oooo00Co0000oooooooooooo00001?ooool00ol0003oooooooooo`2Sooooo`03ool0oooooooooooo 06;ooooo00001Ooooooo0000oooooooooooo000000Wooooo00?o0000ooooooooool01?ooool00ol0 003oooooooooo`2Tooooo`03ool0oooooooooooo067ooooo00001Ooooooo0000oooooooooooo0000 00[ooooo00?o0000ooooooooool00oooool2o`000:Kooooo00?oo`3oooooooooool0H?ooool00005 ooooool0003oooooooooool000001oooool01Ol0003ooooooooooooooooo000000Cooooo00?o0000 ooooooooool0Y_ooool00ooo0?ooooooooooo`1Oooooo`0000Goooooo`000?ooooooooooo`000007 ooooo`05o`000?ooooooooooooooool000001?ooool00ol0003oooooooooo`2Wooooo`03ool0oooo oooooooo05kooooo0002ooooo`;o00002Oooool3o`0000Gooooo00?o0000ooooooooool0Z?ooool0 0ooo0?ooooooooooo`1Mooooo`005Oooool00ol0003oooooooooo`2Yooooo`03ool0oooooooooooo 05cooooo000Eooooo`03o`000?oooooooooo0:[ooooo00?oo`3oooooooooool0Foooool001Gooooo 00?o0000ooooooooool0Zoooool00ooo0?ooooooooooo`1Jooooo`005Oooool00ol0003ooooooooo o`2/ooooo`03ool0oooooooooooo05Wooooo000Eooooo`;o0000[_ooool00ooo0?ooooooooooo`1H ooooo`005Oooool00ol0003oooooooooo`2^ooooo`03ool0oooooooooooo05Oooooo000Eooooo`03 o`000?oooooooooo0:oooooo00?oo`3oooooooooool0E_ooool001Gooooo00?o0000ooooooooool0 /?ooool00ooo0?ooooooooooo`1Eooooo`005Oooool00ol0003oooooooooo`2aooooo`03ool0oooo oooooooo05Cooooo000Eooooo`03o`000?oooooooooo0;;ooooo00?oo`3oooooooooool0Doooool0 01Gooooo00?o0000ooooooooool0/oooool00ooo0?ooooooooooo`1Booooo`005Oooool00ol0003o ooooooooo`2dooooo`03ool0oooooooooooo057ooooo000Eooooo`;o0000]_ooool00ooo0?oooooo ooooo`1@ooooo`005Oooool00ol0003oooooooooo`2eooooo`03ool0oooooooooooo053ooooo000E ooooo`03o`000?oooooooooo0;Kooooo00?oo`3oooooooooool0Coooool001Gooooo00?o0000oooo ooooool0]oooool00ooo0?ooooooooooo`1>ooooo`005Oooool00ol0003oooooooooo`2hooooo`03 ool0oooooooooooo04gooooo000Eooooo`03o`000?oooooooooo0;Wooooo00?oo`3oooooooooool0 C?ooool001Gooooo00?o0000ooooooooool0^_ooool00ooo0?ooooooooooo`1;ooooo`005Oooool0 0ol0003oooooooooo`2kooooo`03ool0oooooooooooo04[ooooo000Eooooo`;o0000_Oooool00ooo 0?ooooooooooo`19ooooo`005Oooool00ol0003oooooooooo`2mooooo`03ool0oooooooooooo04So oooo000Eooooo`03o`000?oooooooooo0;kooooo00?oo`3oooooooooool0Aoooool001Gooooo00?o 0000ooooooooool0_oooool00ooo0?ooooooooooo`16ooooo`005Oooool00ol0003oooooooooo`30 ooooo`03ool0oooooooooooo04Gooooo0002ooooo`;o00001?ooool00ol0003oooooooooo`03oooo o`?o00001?ooool00ol0003oooooooooo`30ooooo`03ool0oooooooooooo04Gooooo00001Ooooooo 0000oooooooooooo000000[ooooo00?o0000ooooooooool00oooool00ol0003oooooooooo`31oooo o`03ool0oooooooooooo04Cooooo00001Ooooooo0000oooooooooooo000000Oooooo1Ol00004oooo o`03o`000?oooooooooo0<;ooooo00?oo`3oooooooooool0@oooool00005ooooool0003ooooooooo ool000001oooool01?l0003oooooooooool00005ooooo`;o0000a?ooool00ooo0?ooooooooooo`12 ooooo`0000Goooooo`000?ooooooooooo`000008ooooo`03o`000?oooooo000000Gooooo00?o0000 ooooooooool0`oooool00ooo0?ooooooooooo`12ooooo`0000Goooooo`000?ooooooooooo`000009 ooooo`;o00001Oooool00ol0003oooooooooo`34ooooo`03ool0oooooooooooo047ooooo0002oooo o`;o00002oooool00ol0003oooooooooo`03ooooo`03o`000?oooooooooo0<Gooooo00?oo`3ooooo ooooool0@?ooool001Gooooo00?o0000ooooooooool0aOooool00ooo0?ooooooooooo`10ooooo`00 5Oooool00ol0003oooooooooo`36ooooo`03ool0oooooooooooo03oooooo000Eooooo`03o`000?oo oooooooo0<Oooooo00?oo`3oooooooooool0?_ooool001Gooooo00?o0000ooooooooool0b?ooool0 0ooo0?ooooooooooo`0mooooo`005Oooool2o`000<Wooooo00?oo`3oooooooooool0?Oooool001Go oooo00?o0000ooooooooool0bOooool00ooo0?ooooooooooo`0looooo`005Oooool00ol0003ooooo ooooo`3:ooooo`03ool0oooooooooooo03_ooooo000Eooooo`03o`000?oooooooooo0<[ooooo00?o o`3oooooooooool0>oooool001Gooooo00?o0000ooooooooool0boooool00ooo0?ooooooooooo`0j ooooo`005Oooool00ol0003oooooooooo`3<ooooo`03ool0oooooooooooo03Wooooo000Eooooo`03 o`000?oooooooooo0<cooooo00?oo`3oooooooooool0>Oooool001Gooooo00?o0000ooooooooool0 cOooool00ooo0?ooooooooooo`0hooooo`005Oooool2o`000<oooooo00?oo`3oooooooooool0=ooo ool001Gooooo00?o0000ooooooooool0c_ooool00ooo0?ooooooooooo`0gooooo`005Oooool00ol0 003oooooooooo`3?ooooo`03ool0oooooooooooo03Kooooo000Eooooo`03o`000?oooooooooo0=3o oooo00?oo`3oooooooooool0=Oooool001Gooooo00?o0000ooooooooool0d?ooool00ooo0?oooooo ooooo`0eooooo`005Oooool00ol0003oooooooooo`3Aooooo`03ool0oooooooooooo03Cooooo000E ooooo`03o`000?oooooooooo0=7ooooo00?oo`3oooooooooool0=?ooool001Gooooo00?o0000oooo ooooool0d_ooool00ooo0?ooooooooooo`0cooooo`005Oooool2o`000=Cooooo00?oo`3ooooooooo ool0<_ooool001Gooooo00?o0000ooooooooool0doooool00ooo0?ooooooooooo`0booooo`005Ooo ool00ol0003oooooooooo`3Dooooo`03ool0oooooooooooo037ooooo000Eooooo`03o`000?oooooo oooo0=Gooooo00?oo`3oooooooooool0<?ooool001Gooooo00?o0000ooooooooool0eOooool00ooo 0?ooooooooooo`0`ooooo`000_ooool2o`0000Cooooo00?o0000ooooooooool00_ooool3o`0000Go oooo00?o0000ooooooooool0e_ooool00ooo0?ooooooooooo`0_ooooo`0000Goooooo`000?oooooo ooooo`000007ooooo`05o`000?ooooooooooooooool000001?ooool00ol0003oooooooooo`3Foooo o`03ool0oooooooooooo02oooooo00001Ooooooo0000oooooooooooo000000Oooooo00Go0000oooo ooooooooooooo`000004ooooo`03o`000?oooooooooo0=Oooooo00?oo`3oooooooooool0;_ooool0 0005ooooool0003oooooooooool000001oooool4o`0000Gooooo0_l0003Iooooo`03ool0oooooooo oooo02gooooo00001Ooooooo0000oooooooooooo000000Oooooo00?o0000ooooooooool01_ooool0 0ol0003oooooooooo`3Hooooo`03ool0oooooooooooo02gooooo00001Ooooooo0000oooooooooooo 000000Sooooo00?o0000ooooooooool01Oooool00ol0003oooooooooo`3Iooooo`03ool0oooooooo oooo02cooooo0002ooooo`;o00002_ooool2o`0000Gooooo00?o0000ooooooooool0fOooool00ooo 0?ooooooooooo`0/ooooo`005Oooool00ol0003oooooooooo`3Jooooo`03ool0oooooooooooo02_o oooo000Eooooo`03o`000?oooooooooo0=[ooooo00?oo`3oooooooooool0:oooool001Gooooo00?o 0000ooooooooool0foooool00ooo0?ooooooooooo`0Zooooo`005Oooool00ol0003oooooooooo`3K ooooo`03ool0oooooooooooo02[ooooo000Eooooo`;o0000gOooool00ooo0?ooooooooooo`0Yoooo o`005Oooool00ol0003oooooooooo`3Mooooo`03ool0oooooooooooo02Sooooo000Eooooo`03o`00 0?oooooooooo0=gooooo00?oo`3oooooooooool0:?ooool001Gooooo00?o0000ooooooooool0g_oo ool00ooo0?ooooooooooo`0Wooooo`005Oooool00ol0003oooooooooo`3Nooooo`03ool0oooooooo oooo02Oooooo000Eooooo`03o`000?oooooooooo0=oooooo00?oo`3oooooooooool09_ooool001Go oooo00?o0000ooooooooool0goooool00ooo0?ooooooooooo`0Vooooo`005Oooool00ol0003ooooo ooooo`3Pooooo`03ool0oooooooooooo02Gooooo000Eooooo`;o0000hOooool00ooo0?oooooooooo o`0Uooooo`005Oooool00ol0003oooooooooo`3Qooooo`03ool0oooooooooooo02Cooooo000Eoooo o`03o`000?oooooooooo0>7ooooo00?oo`3oooooooooool09?ooool001Gooooo00?o0000oooooooo ool0h_ooool00ooo0?ooooooooooo`0Sooooo`005Oooool00ol0003oooooooooo`3Rooooo`03ool0 oooooooooooo02?ooooo000Eooooo`03o`000?oooooooooo0>?ooooo00?oo`3oooooooooool08_oo ool001Gooooo00?o0000ooooooooool0hoooool00ooo0?ooooooooooo`0Rooooo`005Oooool00ol0 003oooooooooo`3Tooooo`03ool0oooooooooooo027ooooo000Eooooo`;o0000i_ooool00ooo0?oo ooooooooo`0Pooooo`005Oooool00ol0003oooooooooo`3Uooooo`03ool0oooooooooooo023ooooo 000Eooooo`03o`000?oooooooooo0>Kooooo00?oo`3oooooooooool07oooool001Gooooo00?o0000 ooooooooool0i_ooool00ooo0?ooooooooooo`0Oooooo`005Oooool00ol0003oooooooooo`3Woooo o`03ool0oooooooooooo01kooooo000Eooooo`03o`000?oooooooooo0>Oooooo00?oo`3ooooooooo ool07_ooool000;ooooo0_l00004ooooo`03o`000?oooooooooo00;ooooo0ol00005ooooo`03o`00 0?oooooooooo0>Sooooo00?oo`3oooooooooool07Oooool00005ooooool0003oooooooooool00000 1oooool01Ol0003ooooooooooooooooo000000Cooooo00?o0000ooooooooool0j?ooool00ooo0?oo ooooooooo`0Mooooo`0000Goooooo`000?ooooooooooo`000007ooooo`05o`000?oooooooooooooo ool000001?ooool2o`000>[ooooo00?oo`3oooooooooool07?ooool00005ooooool0003ooooooooo ool000002?ooool3o`0000Gooooo00?o0000ooooooooool0jOooool00ooo0?ooooooooooo`0Loooo o`0000Goooooo`000?ooooooooooo`000007ooooo`05o`000?ooooooooooooooool000001?ooool0 0ol0003oooooooooo`3Zooooo`03ool0oooooooooooo01_ooooo00001Ooooooo0000oooooooooooo 000000Oooooo00Go0000ooooooooooooooooo`000004ooooo`03o`000?oooooooooo0>[ooooo00?o o`3oooooooooool06oooool000;ooooo0_l00009ooooo`?o00001Oooool00ol0003oooooooooo`3[ ooooo`03ool0oooooooooooo01[ooooo000Eooooo`03o`000?oooooooooo0>_ooooo00?oo`3ooooo ooooool06_ooool001Gooooo00?o0000ooooooooool0k?ooool00ooo0?ooooooooooo`0Iooooo`00 5Oooool00ol0003oooooooooo`3/ooooo`03ool0oooooooooooo01Wooooo000Eooooo`;o0000k_oo ool00ooo0?ooooooooooo`0Hooooo`005Oooool00ol0003oooooooooo`3]ooooo`03ool0oooooooo oooo01Sooooo000Eooooo`03o`000?oooooooooo0>kooooo00?oo`3oooooooooool05oooool001Go oooo00?o0000ooooooooool0k_ooool00ooo0?ooooooooooo`0Gooooo`005Oooool00ol0003ooooo ooooo`3^ooooo`03ool0oooooooooooo01Oooooo000Eooooo`03o`000?oooooooooo0>oooooo00?o o`3oooooooooool05_ooool001Gooooo00?o0000ooooooooool0koooool00ooo0?ooooooooooo`0F ooooo`005Oooool00ol0003oooooooooo`3`ooooo`03ool0oooooooooooo01Gooooo000Eooooo`;o 0000lOooool00ooo0?ooooooooooo`0Eooooo`005Oooool00ol0003oooooooooo`3aooooo`03ool0 oooooooooooo01Cooooo000Eooooo`03o`000?oooooooooo0?7ooooo00?oo`3oooooooooool05?oo ool001Gooooo00?o0000ooooooooool0l_ooool00ooo0?ooooooooooo`0Cooooo`005Oooool00ol0 003oooooooooo`3booooo`03ool0oooooooooooo01?ooooo000Eooooo`03o`000?oooooooooo0??o oooo00?oo`3oooooooooool04_ooool001Gooooo00?o0000ooooooooool0loooool00ooo0?oooooo ooooo`0Booooo`005Oooool00ol0003oooooooooo`3dooooo`03ool0oooooooooooo017ooooo000E ooooo`;o0000mOooool00ooo0?ooooooooooo`0Aooooo`005Oooool00ol0003oooooooooo`3eoooo o`03ool0oooooooooooo013ooooo000Eooooo`03o`000?oooooooooo0?oooooo2Oooool001Gooooo 00?o0000ooooooooool0ooooool9ooooo`005Oooool00ol0003oooooooooo`3oooooo`Wooooo000E ooooo`03o`000?oooooooooo0?oooooo2Oooool00?oooooo8Oooool00?oooooo8Oooool00?oooooo 8Oooool00001\ \>"], ImageRangeCache->{{{0, 287}, {176.938, 0}} -> {-12017.7, -0.0984178, \ 565.126, 0.00621951}}], Cell[BoxData[ TagBox[\(\[SkeletonIndicator] Graphics \[SkeletonIndicator]\), False, Editable->False]], "Output"] }, Open ]], Cell[BoxData[ \(But\ now\ I\ want\ to\ be\ able\ to\ take\ the\ values\ at\ the\ end\ \ condition, \ and\ feed\ them\ back\ into\ the\ chemoton\ function\ slightly\ changed . \ \ Please\ can\ someone\ show\ me\ how\ to\ do\ this \(\(??\)\(\ \)\)\)], \ "Input"] }, FrontEndVersion->"5.1 for Macintosh", ScreenRectangle->{{0, 1024}, {0, 706}}, WindowSize->{1020, 560}, WindowMargins->{{Automatic, 0}, {Automatic, 4}} ] (******************************************************************* 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, 11635, 185, 1368, "Input"], Cell[13414, 240, 358, 5, 21, "Message"], Cell[13775, 247, 358, 5, 21, "Message"], Cell[14136, 254, 365, 5, 21, "Message"], Cell[14504, 261, 364, 5, 21, "Message"], Cell[14871, 268, 367, 5, 21, "Message"], Cell[15241, 275, 332, 5, 21, "Message"], Cell[15576, 282, 359, 5, 21, "Message"], Cell[15938, 289, 359, 5, 21, "Message"], Cell[16300, 296, 359, 5, 21, "Message"], Cell[16662, 303, 354, 5, 21, "Message"], Cell[17019, 310, 333, 5, 21, "Message"], Cell[17355, 317, 5734, 165, 235, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[23126, 487, 390, 7, 75, "Input"], Cell[23519, 496, 27683, 838, 186, 9122, 604, "GraphicsData", "PostScript", \ "Graphics"], Cell[51205, 1336, 130, 3, 27, "Output"] }, Open ]], Cell[51350, 1342, 269, 5, 43, "Input"] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************) But now I want to be able to take the values at the end condition, and feed \ them back into the chemoton function slightly changed. Please can someone \ show me how to do this?? Yours ever, Chrisantha www.chrisantha.com
- References:
- Re: Representation and Simulation of Dynamic Systems
- From: "Caffa Vittorio Dr." <Caffa@iabg.de>
- Re: Representation and Simulation of Dynamic Systems