Multiple Executions of .nb file part2

• To: mathgroup at smc.vnet.net
• Subject: [mg90484] Multiple Executions of .nb file part2
• From: "Dave" <david.fourie at ncl.ac.uk>
• Date: Thu, 10 Jul 2008 06:38:43 -0400 (EDT)
• Organization: Newcastle University

```Hi,

I think I maybe over simplified my situation. My project uses a series of
nested differential equations. These run over a given time and then print
either a graph or the plot points or both (depending on my needs).

An example of the code is:

****************************

Off[General::spell]

time = 500;

SmV = -80;

PD2mnax0[PD2u_] := 1/(1 + Exp[-(PD2u + 24.7)/5.29]);

PD2hnax0[PD2u_] := 1/(1 + Exp[(PD2u + 48.9)/5.18]);

PD2mk10x[PD2u_] := 1/(1 + Exp[-(PD2u + 14.2)/11.8]);

PD2taumnax[PD2u_] := 1.32 - 1.26/(1 + Exp[-(PD2u + 120)/25]);

PD2tauhnax[PD2u_] := (0.67/(1 + Exp[-(PD2u + 62.9)/10]))*(1.5 + 1/(1 + Exp[(

PD2u + 34.9)/3.6]));

PD2taumk1[PD2u_] := 7.2 - 6.4/(1 + Exp[-(PD2u + 28.3)/19.2]);

solution = NDSolve[{

PD2mnax'[t] == (PD2mnax0[PD2z[t]] - PD2mnax[t])*(1/

PD2taumnax[PD2z[t]]), PD2mnax[0] == PD2mnax0[-80]*(1/

PD2taumnax[-80]), PD2hnax'[t] == (PD2hnax0[PD2z[t]] - PD2hnax[

t])*(1/PD2tauhnax[PD2z[

t]]), PD2hnax[

0] == PD2hnax0[-80]*(1/PD2tauhnax[-80]),

PD2mkx'[t] == (PD2mk10x[PD2z[t]] - \

PD2mkx[t])*(1/PD2taumk1[PD2z[t]]),

PD2mkx[0] == PD2mk10x[-80]*(1/PD2taumk1[-80]), PD2z'[t] == (-1*(0

+ 1*1100*(PD2mnax[t]^3)*PD2hnax[t]*(PD2z[t] - 50)

+ 1*150*(PD2mkx[t]^4)*(PD2z[t] + 80)

+ 1*0.00081*(PD2z[t] + 50)

)

*(1/6)),

PD2z[0] == SmV

},

{PD2z, PD2mnax,

PD2hnax, PD2mkx}, {t, 0, time}, Method -> Automatic, MaxSteps -> \

Infinity];

ParametricPlot[Evaluate[{t, PD2z[

t]} /. solution], {t, 0, time},

PlotPoints -> time, PlotRange -> {-80, 50}, PlotLabel -> "PD Axon"];

(*

mkw = Table[0, {time-1}];

For[w=1,w<(time),mkw[[w]]=PD2z[w]/.solution;w++];

mkw

*)

*******************************

This produces what looks like a flat line on a graph.

I need to run this file several times changing the variable SmV. Can I
encapsulate the whole thing in a function? I have tried but it doesn't work.

fun[]:= (SmV;

***code***

)

fun[-80]

out- fun[-80]

Thank you

David Fourie

```

• Prev by Date: Re: distribution fitting
• Next by Date: Glitch in Exporting PDFs from Mathematica to Illustrator, one solution
• Previous by thread: Re: How to do that with version 6?
• Next by thread: Re: Multiple Executions of .nb file part2