Re: capacitor equation solution not working
- To: mathgroup at smc.vnet.net
- Subject: [mg119716] Re: capacitor equation solution not working
- From: "J and B" <noslowski at comcast.net>
- Date: Sat, 18 Jun 2011 06:22:14 -0400 (EDT)
- References: <201106170409.AAA28826@smc.vnet.net> <201106170953.FAA02174@smc.vnet.net>
Sorry for not typing my code in correctly. I have the same code as you,
that is mentioned in this email
i[t_] := If[t < 10 || t > 60, 0, 1]
eq = v'[t] == i[t]/c;
c = 1;
soleq = v[t] /. First@NDSolve[{eq, v[0] == 0}, v[t], {t, 0, 100}];
Plot[soleq, {t, 0, 100}]
But I thought it had to be wrong, since the capacitor never came back to
zero. I did not realize that I had to have a resistor in the circuit, as
several people pointed out.
Thanks for your help
-----Original Message-----
From: DrMajorBob [mailto:btreat1 at austin.rr.com]
Sent: Friday, June 17, 2011 5:53 AM
To: mathgroup at smc.vnet.net
Subject: [mg119716] Re: capacitor equation solution not working
Surely that is NOT your code; it has at least four misspellings and two
fatal syntax errors -- == where you need = (as in Soleq = ...) and the comma
before [[1]].
If you ran that code, you'd see error messages, not just a constant plot.
This works, however:
i[t_] := If[t < 10 || t > 60, 0, 1]
eq = v'[t] == i[t]/c;
c = 1;
soleq = v[t] /. First@NDSolve[{eq, v[0] == 0}, v[t], {t, 0, 100}];
Plot[soleq, {t, 0, 100}]
I'd also suggest THIS for the first line:
i[t_] := Boole[10 <= t <= 60]
Bobby
On Thu, 16 Jun 2011 23:09:13 -0500, J and B <noslowski at comcast.net> wrote:
> I am trying to solve a very simple circuit, one that has only a
> capacitor in it. What am I doing wrong? The final graph should show
> that the capacitor has discharged, but its staying constant. I have
> current starting at t=10 and ending at t=60.
> I hope I have used the correct way to put in my code, I am using a new
> computer and I having trouble with my email program. Any advice would
> be appreciated.
>
> Thanks
>
> My code.
> I[t_]:= If[t<10||t>60,0,1]
> Plot[i[t],{t,0,100}]
> Eq=v'[t]==i[t]/c
> c=1;
> Soleq==v[t]/.NDSolve[{eq,v[0]==0},v[t],{t,0,100}],[[1]]
> Plot[soleq,{t,0,100}]
--
DrMajorBob at yahoo.com
- References:
- capacitor equation solution not working
- From: "J and B" <noslowski@comcast.net>
- Re: capacitor equation solution not working
- From: DrMajorBob <btreat1@austin.rr.com>
- capacitor equation solution not working