MathGroup Archive 2009

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

Search the Archive

Re: Piecewise with ODE and constant

  • To: mathgroup at smc.vnet.net
  • Subject: [mg104770] Re: Piecewise with ODE and constant
  • From: sean <sean_incali at yahoo.com>
  • Date: Mon, 9 Nov 2009 05:45:33 -0500 (EST)
  • References: <hd6c4q$p43$1@smc.vnet.net>

>  eq=v'[t]*1/c *(i-v[t]/r)

This is not an equation. DSolve requires an equation.
v[t]==v[t] for instance.


>  r=100 x 10^6;i=0.3 x 10^-9;c=100 x 10^-12;
>

x is not multiplication. * is used instead. So r= 100*10^6


>  ans=DSolve[{eq,v[0]==0},v[t],t]
>
>  p1=Plot[v[t]/.ans,{t,0,140 x 10^-3}]

Again x is not multiplication. 140*10^-3


>
>  Plot[Piecewise[{{-70 x 10^-3,0*t*10 x 10^-3},{ans,10 x 10^-3*t*100 x 10^-3},{-70 x 10^-3,100 x 10^-3*t*140 x 10^-3}}],{t,0,140 x 10^-3}]


I'm not sure what you want above. If you want to define an ode
piecewise then you have to define is like this. But what yousaid you
want is i0[t_]:= i0[t]/; 10 <= t <= 100;

i0[t_] := 0/; t < 0;
i0[t_] := 0/; 0 <= t <= 10;
i0[t_] := 1/; 10 <= t <= 100;
deqn = i0[t] == v[t]/r + c v'[t];

NDSolve will solve that.



On Nov 8, 4:04 am, "Becky" <noslow... at comcast.net> wrote:
> I have tried so many different ways to do this problem, that I am stuck.
> Its like being back I school.
>
> I would like to plot my ODE which is a curve that goes starts at 10 msec and
> ends at 100 msec).  In the beginning there is no current delivered to t he
> cell until time equal 10 msec.  Once the current probe in inserted to the
> cell, I have a ODE, which is modeled as an RC circuit.  When time equal 100
> msec, the current it turned off, and the curve should go back to the
> original state.  I have tried using If statements, Piecewise plots, and so
> on, but I cannot get the curve to work with the initial and final conditions
> of -70 mV.
>
> I would like to stay with solving the ODE using Dsolve, vs., NDSolve, if
> that is ok.  I cannot figure out my problem, I need help, please.
>
> Sincerely Yours
> Prof. Jake
>
>  eq=v'[t]*1/c *(i-v[t]/r)
>
>  r=100 x 10^6;i=0.3 x 10^-9;c=100 x 10^-12;
>
>  ans=DSolve[{eq,v[0]==0},v[t],t]
>
>  p1=Plot[v[t]/.ans,{t,0,140 x 10^-3}]
>
>  Plot[Piecewise[{{-70 x 10^-3,0*t*10 x 10^-3},{ans,10 x 10^-3*t*100 x 10^-3},{-70 x 10^-3,100 x 10^-3*t*140 x 10^-3}}],{t,0,140 x 10^-3}]



  • Prev by Date: Give arrow a certai color
  • Next by Date: Re: Thikness of edges in a 3D graph
  • Previous by thread: Re: Piecewise with ODE and constant
  • Next by thread: Re: How can I reduce the frequency of a sine wave tia