MathGroup Archive 1998

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

Search the Archive

DSolve&NDSOLve with 3.0.1, w95


I need some advice! 
Both the instructions below return solutions with math 2.2.1. 3.0.1
returns excuses by the page full. Am I missing something, or has
mathematica lost something between 2.2.1 and 3.0.1?  I suspect the
fault is mine since I have not noted discussion of this problem here or
otherwise. 
I would appreciate any help anyone can offer. meanwhile, I will
carefully retain 2.2.1..... 

gv=1;gp=.02;gr=.0001;ar=1;ai=.99;
DSolve[{ac'''[t]+gv ac''[t]+gp ac'[t]+gr ac[t]== (ai-ar) gr, ac[0]==0,
 ac'[0]==(ai-ar) gv, ac''[0]==(ai-ar) gp-gv ac'[0]},  ac[t], t]


Clear[ai,ar,t,v];kv=1;ka=.2;kj=.0016; vt=NDSolve[{v''''[t]+kv v'''[t]+ka
v''[t]+kj v'[t]== ka ar'[t]+kj ar[t], p''''[t]+kv p'''[t]+ka p''[t]==
v'''[t]+kv ar'[t]+ka ar[t],  ar'[t]==.1, v[0]==0,  v'[0]==.9,
v''[0]==.09+kv-.9 kv, v'''[0]==ar'[0] kv +ka ar[0]-v''[0] kv-v'[0] ka, 
ar''[0]==0,  ar'[0]==.1, ar[0]==1, p'''[0]==.29, p''[0]==.9, p'[0] ==0,
p[0]==0}, 
{v[t], p[t], ar[t]}, {t, 0, 1000}]
-- 
ned!


  • Prev by Date: error messages
  • Next by Date: Re: EPS file of Plot with AxesLabel->{"t","f(t)"}]
  • Previous by thread: Re: error messages
  • Next by thread: Is there an online documentation ???