(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 4.2' 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[ 19765, 352]*) (*NotebookOutlinePosition[ 20409, 374]*) (* CellTagsIndexPosition[ 20365, 370]*) (*WindowFrame->Normal*) Notebook[{ Cell[BoxData[ \(\(\( \ (*\(--\(--\(--\(--\(--\(--\(--\(-\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(\ --\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(----\)\ \)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)*) \)\ \(\[IndentingNewLine]\)\( (*\ \ MATHEMATICA\ CODE\ to\ run\ the\ model\ \ described\ in\ Goodreau\ et\ al . \ JID\ 2003. \ *) \)\(\[IndentingNewLine]\)\( (*\ \ Last\ \(updated : \ \ \ 8\ May\ 2003\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ *) \)\(\[IndentingNewLine]\)\( (*\(--\(\ --\(--\(--\(--\(--\(--\(--\(--\(--\(-\(--\(--\(--\(--\(--\(--\(--\(--\(--\(\ --\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(----\)\)\)\)\)\ \)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)*) \)\(\ \[IndentingNewLine]\)\(\[IndentingNewLine]\)\( (*\ \ INTRO : \ \[IndentingNewLine]\[IndentingNewLine]\ \ \ The\ basic\ structure\ is\ a\ \ compartmental\ model\ defined\ by\ a\ system\ of\ six\ ODE' s, \ one\ per\ compartment . \ \[IndentingNewLine]\[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ in[ t]\ = \ \(#\ of\ insertive\ negatives\ at\ time\ t\ \[IndentingNewLine] \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ ip[ t]\ = \ \(#\ of\ insertive\ positives\ at\ time\ t\ \[IndentingNewLine] \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ rn[ t]\ = \ \(#\ of\ receptive\ negatives\ at\ time\ t\ \[IndentingNewLine] \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ rp[ t]\ = \ \(#\ of\ receptive\ positives\ at\ time\ t\ \[IndentingNewLine] \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ vn[ t]\ = \ \(#\ of\ versatile\ negatives\ at\ time\ t\ \[IndentingNewLine] \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ vp[ t]\ = \ #\ of\ versatile\ positives\ at\ time\ t\)\)\)\)\ \)\[IndentingNewLine]*) \)\(\[IndentingNewLine]\)\(\[IndentingNewLine]\)\( \ (*\(--\(--\(--\(--\(--\(--\(--\(-\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(\ --\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(----\)\ \)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)*) \)\ \(\[IndentingNewLine]\)\( \ (*\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(-PARAMETERS\)\)\)\)\)\)\)\ \)\)\)\)\)\)\ TO\ \ \(\(\(\(\(\(\(\(\(\(\(\(\(\(\(CHANGE--\)--\)--\)--\)--\)--\)--\)--\)--\)--\)\ --\)--\)--\)--\)--\)*) \)\(\[IndentingNewLine]\)\( (*\(--\(--\(--\(--\(--\(\ --\(--\(-\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(\ --\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(----\)\)\)\)\)\)\)\)\)\)\)\)\)\ \)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)*) \)\(\[IndentingNewLine]\)\(\ \[IndentingNewLine]\)\(\(Off[ General::spell1];\)\[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \[IndentingNewLine] \(singlerun\ = \ True;\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ If\ True, \ a\ single\ run\ is\ conducted\ using\ the\ user - provided\ values\ "\", \ "\", \[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ and\ "\" . \ Returned\ values\ are\ the\ size\ of\ each\ compartment\ every\ \ "\"\[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ number\ of\ days\ for\ the\ time\ for\ the\ time\ period\ of\ length\ \ "\" . \[IndentingNewLine]\[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ If\ False, \ a\ full\ set\ of\ 605\ runs\ are\ conducted\ using\ each\ combination\ \ of\ the\ \(values : \[IndentingNewLine]\[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ vvalue\)\ = \ 0, \ 0.1, \ 0.2, \ 0.3, \ 0.4, \ 0.5, \ 0.6, \ 0.7, \ 0.8, \ 0.9, \ 1.0\ \[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ hvalue\ = \ 0, \ 0.1, \ 0.2, \ 0.3, \ 0.4, \ 0.5, \ 0.6, \ 0.7, \ 0.8, \ 0.9, \ 1.0\ \[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ cvalue\ = \ 18, \ 21, \ 24, \ 27, \ 30. \[IndentingNewLine]\[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ Returned\ values\ are\ the\ size\ of\ each\ compartment\ at\ \ the\ end\ of\ each\ simulation\ \[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \((time = \ "\")\)\ for\ each\ run . \ \ \[IndentingNewLine]\[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \(REMINDER : \ BE\ AWARE\ OF\ CASE - SENSITIVITY\ IN\ CHOOSING\ True\ OR\ False\)\ \ *) \[IndentingNewLine]\[IndentingNewLine] \(vvalue\ = \ 1.0;\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ fraction\ of\ population\ that\ is\ versatile; \ ranges\ from\ 0\ to\ 1; \ matters\ only\ if\ singlerun\ = \ True\ *) \[IndentingNewLine] \(hvalue\ = \ 0.5;\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ role\ mixing\ level; \ ranges\ from\ 0\ to\ 1; \ matters\ only\ if\ singlerun\ = \ True\ *) \[IndentingNewLine] \(cvalue\ = \ 24;\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ average\ number\ of\ partners\ per\ person\ per\ year\ \((actually\ per\ \ 360\ days)\); \ can\ take\ any\ non - negative\ value; \ matters\ only\ if\ singlerun\ = \ True\ \ *) \[IndentingNewLine]\ \[IndentingNewLine] \(duration\ = 20*365;\)\ \ \ \ \ \ \ \ \ \ \ (*\ duration\ in\ days; \ note, \ if\ less\ than\ 20*365\ then\ an\ error\ will\ occur\ for\ \ multiple\ runs\ *) \ \ \ \ \ \ \ \ \ \ \ \ \[IndentingNewLine] \(samplewidth = 1000;\)\ \ \ \ \ \ \ \ \ \ \ \ \ (*\ number\ of\ days\ between\ samples\ for\ a\ single\ run\ *) \[IndentingNewLine] \(initn = 10000;\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ initial\ total\ pop\ size\ *) \[IndentingNewLine] \(initp = 0.01;\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ initial\ HIV\ prevalence\ *) \[IndentingNewLine] \(\[Mu] = 2;\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ number\ of\ new\ susceptibles\ per\ day\ *) \[IndentingNewLine] \(\[Delta] = \(1/365\)/ 30;\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ fraction\ removed\ per\ day\ from\ general\ causes\ \((all\ 6\ \ compartments)\)\ \ *) \n \(z = \(1/365\)/10;\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ fraction\ removed\ per\ day\ from\ AIDS\ \((3\ postive\ compartments\ \ only)\)\ \ \ *) \n \(infi = 0.005;\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ prob\ of\ infection\ per\ day\ for\ a\ positive - receptive\ negative - insertive\ contact\ *) \n \(infr = 0.025;\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ prob\ of\ infection\ per\ day\ for\ a\ positive - insertive\ negative - receptive\ contact\ *) \[IndentingNewLine]\[IndentingNewLine] \ (*\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(\ --\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(----\)\)\)\)\)\)\)\)\)\)\)\ \)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)*) \[IndentingNewLine] (*\ CHANGES\ BELOW\ THIS\ POINT\ WILL\ AFFECT\ THE\ STRUCTURE\ OF\ THE\ \ MODEL\ \ *) \[IndentingNewLine] (*\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(\ --\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(\ --\(--\(----\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\)\ *) \[IndentingNewLine]\[IndentingNewLine]\[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ defining\ group\ relationships\ *) \[IndentingNewLine]\ \[IndentingNewLine] \(i[t] = in[t] + ip[t];\)\n \(r[t] = rn[t] + rp[t];\)\n \(v[t] = vn[t] + vp[t];\)\n \(n[t] = i[t] + r[t] + v[t];\)\[IndentingNewLine]\[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ select\ run\ values\ *) \[IndentingNewLine]\[IndentingNewLine] \(If\ [ singlerun == \ True, \ {lowvindex = vvalue, highvindex = vvalue, lowhindex = hvalue, highhindex = hvalue, lowcindex = cvalue, \ highcindex = cvalue}, {lowvindex = 1, highvindex = 11, lowhindex = 1, highhindex = 11, lowcindex = 1, \ highcindex = 5}];\)\[IndentingNewLine]\[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ \(--\(--\(\ --\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\(-main\)\)\)\)\)\)\)\)\)\)\)\)\ \)\)\)\ \(\(\(\(\(\(\(\(\(\(\(\(\(\(\(code\ --\)--\)--\)--\)--\)--\)--\)--\)\ --\)--\)--\)--\)--\)--\)--\)*) \[IndentingNewLine]\[IndentingNewLine] \(Do[\[IndentingNewLine]\ \[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ for\ a\ single\ run, \ looping\ indices\ do\ not\ have\ to\ be\ integral; \ \[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ cindex, \ vindex, \ and\ hindex\ can\ be\ given\ whatever\ values\ we\ want\ for\ c, v, h\[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ for\ \ multiple\ runs, \ looping\ indices\ must\ be\ integral\ because\ they\ are\ used\ as\ \ array\ positions; \[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ cindex, \ vindex\ and\ hindex\ must\ therefore\ be\ transformed\ before\ \ being\ assigned\ to\ c, v, h\ *) \[IndentingNewLine]\[IndentingNewLine]\ If[singlerun \[Equal] True, \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \[IndentingNewLine]{\[IndentingNewLine]\ \ c = cindex/360, \[IndentingNewLine]\ \ initv = Which[vindex \[Equal] 0, 0.00001, vindex \[Equal] 1, .99999, True, vindex], \[IndentingNewLine]\ \ h = \((hindex)\)\ \[IndentingNewLine]}, \n\ \ \ \ \ \ \ \ \ \ {\ \[IndentingNewLine]\ c = \((cindex + 5)\)/120, \[IndentingNewLine]\ \ initv = Which[vindex \[Equal] 1, 0.00001, vindex \[Equal] 11, .99999, True, \((vindex - 1)\)/ 10], \[IndentingNewLine]\ \ h = \((hindex - 1)\)/ 10\n\ \ \ \ \ \ \ \ \ \ \ }, \n\ \ \ \ \ \ \ \ ]; \ \[IndentingNewLine]\[IndentingNewLine]initi = \((1 - initv)\)/ 2; \[IndentingNewLine]initr = \((1 - initv)\)/ 2; \n\ \ \ \[IndentingNewLine]\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (*\ \ \ \ \ \ \ \ \ ODE\ core\ \ \ \ \ \ \ \ \ *) \[IndentingNewLine]\ \[IndentingNewLine]Solution = NDSolve[{\(in'\)[ t] \[Equal] \[Mu]* initi - \((infi*in[t]* rp[t]*\(c/i[t]\)/ r[t])\) \((\((1 - h)\) \((n[t]/2 + \(v[t]^2/2\)/n[t] - v[t])\) + \((h)\) \((i[t] + r[t])\)/ 2)\) - \((infi*in[t]* vp[t]*\(c/i[t]\)/v[t])\) \((1 - h)\) \((i[t]* v[t]/n[t])\) - \[Delta]*in[t], \(rn'\)[ t] \[Equal] \[Mu]* initr - \((infr*rn[t]* ip[t]*\(c/i[t]\)/ r[t])\) \((\((1 - h)\) \((n[t]/2 + \(v[t]^2/2\)/n[t] - v[t])\) + \((h)\) \((i[t] + r[t])\)/ 2)\) - \((infr*rn[t]* vp[t]*\(c/r[t]\)/v[t])\) \((1 - h)\) \((r[t]* v[t]/n[t])\) - \[Delta]*rn[t], \(vn'\)[ t] \[Equal] \[Mu]* initv - \((infr*vn[t]*ip[t]*\(c/v[t]\)/i[t])\) \((1 - h)\) \((v[t]*i[t]/n[t])\) - \((infi*vn[t]* rp[t]*\(c/v[t]\)/r[t])\) \((1 - h)\) \((v[t]* r[t]/n[t])\) - \((infi + infr)\)*\((vn[t]*vp[t]* c/v[t]^2)\)*\((\((\(v[t]^2/2\)/n[t])\) \((1 - h)\) + v[t]*h/2)\) - \[Delta]*vn[t], \(ip'\)[ t] \[Equal] \((infi*in[t]* rp[t]*\(c/i[t]\)/ r[t])\) \((\((1 - h)\) \((n[t]/2 + \(v[t]^2/2\)/n[t] - v[t])\) + \((h)\) \((i[t] + r[t])\)/ 2)\) + \((infi*in[t]* vp[t]*\(c/i[t]\)/v[t])\) \((1 - h)\) \((i[t]* v[t]/n[t])\) - \[Delta]*ip[t] - z*ip[t], \(rp'\)[ t] \[Equal] \((infr*rn[t]* ip[t]*\(c/i[t]\)/ r[t])\) \((\((1 - h)\) \((n[t]/2 + \(v[t]^2/2\)/n[t] - v[t])\) + \((h)\) \((i[t] + r[t])\)/ 2)\) + \((infr*rn[t]* vp[t]*\(c/r[t]\)/v[t])\) \((1 - h)\) \((r[t]* v[t]/n[t])\) - \[Delta]*rp[t] - z*rp[t], \(vp'\)[ t] \[Equal] \((infr*vn[t]*ip[t]*\(c/v[t]\)/i[t])\) \((1 - h)\) \((v[t]*i[t]/n[t])\) + \((infi*vn[t]* rp[t]*\(c/v[t]\)/r[t])\) \((1 - h)\) \((v[t]* r[t]/n[t])\) + \((infi + infr)\)*\((vn[t]*vp[t]* c/v[t]^2)\)*\((\((\(v[t]^2/2\)/n[t])\) \((1 - h)\) + v[t]*h/2)\) - \[Delta]*vp[t] - z*vp[t], \[IndentingNewLine]in[0] \[Equal] initn*initi*\((1 - initp)\), rn[0] \[Equal] initn*initr*\((1 - initp)\), vn[0] \[Equal] initn*initv*\((1 - initp)\), ip[0] \[Equal] initn*initi*initp, rp[0] \[Equal] initn*initr*initp, vp[0] \[Equal] initn*initv*initp}, {in, rn, vn, ip, rp, vp}, {t, 0, duration}]; \[IndentingNewLine]\[IndentingNewLine]ineg\ = \ Solution[\([1, 1, 2]\)]; \[IndentingNewLine]rneg\ = \ Solution[\([1, 2, 2]\)]; \[IndentingNewLine]vneg\ = \ Solution[\([1, 3, 2]\)]; \[IndentingNewLine]ipos\ = \ Solution[\([1, 4, 2]\)]; \[IndentingNewLine]rpos\ = \ Solution[\([1, 5, 2]\)]; \[IndentingNewLine]vpos\ = \ Solution[\([1, 6, 2]\)]; \[IndentingNewLine]\[IndentingNewLine]If[ singlerun == False, \[IndentingNewLine]{ineg20[cindex, vindex, hindex]\ = \ ineg[7300], \[IndentingNewLine]rneg20[cindex, vindex, hindex]\ = \ rneg[7300], \[IndentingNewLine]vneg20[cindex, vindex, hindex]\ = \ vneg[7300], \[IndentingNewLine]ipos20[cindex, vindex, hindex]\ = \ ipos[7300], \[IndentingNewLine]rpos20[cindex, vindex, hindex]\ = \ rpos[7300], \[IndentingNewLine]vpos20[cindex, vindex, hindex]\ = \ vpos[ 7300], \[IndentingNewLine]}], \[IndentingNewLine]\ \[IndentingNewLine]{cindex, lowcindex, highcindex}, {vindex, lowvindex, highvindex}, {hindex, lowhindex, highhindex}];\)\[IndentingNewLine]\[IndentingNewLine]\ \[IndentingNewLine] (*\(--\(--\(--\(--\(--\(--\(--\(--\(--\(--\ \ \ \ Final\)\ \)\)\)\)\)\)\)\)\)\ processing\ into\ \(tables\ \ \ \ \ \ \ --\) \(--\(--\(\ --\(--\(--\(--\(--\(--\(--\(--\(--\(---\)\)\)\)\)\)\)\)\)\)\)\)\ *) \[IndentingNewLine]\[IndentingNewLine] If\ [singlerun \[Equal] False, \[IndentingNewLine]{inegAllRunsY20Table\ = \ Table[ineg20[i, j, k], {i, 5}, {j, 11}, {k, 11}], \[IndentingNewLine]rnegAllRunsY20Table\ = \ Table[rneg20[i, j, k], {i, 5}, {j, 11}, {k, 11}], \[IndentingNewLine]vnegAllRunsY20Table\ = \ Table[vneg20[i, j, k], {i, 5}, {j, 11}, {k, 11}], \[IndentingNewLine]iposAllRunsY20Table\ = \ Table[ipos20[i, j, k], {i, 5}, {j, 11}, {k, 11}], \[IndentingNewLine]rposAllRunsY20Table\ = \ Table[rpos20[i, j, k], {i, 5}, {j, 11}, {k, 11}], \[IndentingNewLine]vposAllRunsY20Table\ = \ Table[vpos20[i, j, k], {i, 5}, {j, 11}, {k, 11}], \[IndentingNewLine]}, \[IndentingNewLine]{\ \[IndentingNewLine]inegOneRunTable\ = \ Table[ineg[i], {i, 0, duration, samplewidth}], \[IndentingNewLine]rnegOneRunTable\ = \ Table[rneg[i], {i, 0, duration, samplewidth}], \[IndentingNewLine]vnegOneRunTable\ = \ Table[vneg[i], {i, 0, duration, samplewidth}], \[IndentingNewLine]iposOneRunTable\ = \ Table[ipos[i], {i, 0, duration, samplewidth}], \[IndentingNewLine]rposOneRunTable\ = \ Table[rpos[i], {i, 0, duration, samplewidth}], \[IndentingNewLine]vposOneRunTable\ = \ Table[vpos[i], {i, 0, duration, samplewidth}]\[IndentingNewLine]}\ \[IndentingNewLine]\ ]\ \[IndentingNewLine]\[IndentingNewLine]\[IndentingNewLine]\[IndentingNewLine] \ (*\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ OK, \ now\ it' s\ up\ to\ you\ to\ export\ them\ however\ you\ \(want!\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ *) \ \[IndentingNewLine]\[IndentingNewLine] \)\)\)], "Input"] }, FrontEndVersion->"4.2 for Microsoft Windows", ScreenRectangle->{{0, 1280}, {0, 908}}, WindowSize->{1272, 876}, 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[1754, 51, 18007, 299, 3150, "Input"] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)