Re: Writing to an External file, ie .txt
- To: mathgroup at smc.vnet.net
- Subject: [mg90532] Re: Writing to an External file, ie .txt
- From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
- Date: Sat, 12 Jul 2008 05:31:01 -0400 (EDT)
- Organization: The Open University, Milton Keynes, UK
- References: <g56ti4$41a$1@smc.vnet.net>
Dave wrote:
> I have a matica simulation that generates a Table. Currently it prints in at
> the bottom of the notebook and then I copy and paste it into a text file.
>
> below is the code that generates the Table.
>
> **************************************************
> 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] == -80
>
> },
>
>
>
> {PD2z, PD2mnax, PD2hnax, PD2mkx}, {t, 0, time}, Method -> Automatic,
> MaxSteps -> Infinity];
>
>
> mkw = Table[0, {time - 1}];
>
>
> For[w = 1, w < (time), mkw[[w]] = PD2z[w] /. solution; w++];
>
>
> mkw
>
> ******************************************************
What about *exporting* the table?
[...]
In[11]:= mkw
Out[11]= {{-79.996}, {-79.9919}, {-79.9879}, {-79.9838}, {-79.9798},
{-79.9757}, \
{-79.9717}, {-79.9676}, {-79.9636}}
In[12]:= Export["le.txt", mkw, "Table"]
Out[12]= "le.txt"
In[13]:= FilePrint["le.txt"]
During evaluation of In[13]:=
-79.99595026144705
-79.99190108087485
-79.9878524469452
-79.98380435954662
-79.9797568186066
-79.97570982405419
-79.97166337581847
-79.96761747382855
-79.96357211801356
In[14]:= Import["le.txt", "Table"]
Out[14]= {{-79.996}, {-79.9919}, {-79.9879}, {-79.9838}, {-79.9798},
{-79.9757}, \
{-79.9717}, {-79.9676}, {-79.9636}}
In[15]:= Directory[]
Out[15]= "/Users/jean-marcgulliet"
Regards,
-- Jean-Marc