MathGroup Archive 2008

[Date Index] [Thread Index] [Author Index]

Search the Archive

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