PHYS 115 - Rec 5 - Solution Key - Travis Hoppe
restart:
R := 1: V:= 12: C[1] := 1: C[2] := 2: C[3] := 3:
deq[1] := V - diff(q[0](t),t)*R - diff(q[1](t),t)*R - q[1](t)/C[1]=0;
deq[2] := V - diff(q[0](t),t)*R - diff(q[2](t),t)*R - q[2](t)/C[2]=0;
deq[3] := V - diff(q[0](t),t)*R - diff(q[3](t),t)*R - q[3](t)/C[3]=0;
deq[4] := diff(q[0](t),t) = diff(q[1](t),t) + diff(q[2](t),t) + diff(q[3](t),t);
ICS := q[1](0)=0, q[2](0)=0, q[3](0)=0;
ODE := deq[i]$i=1..4;
LUklaGVscEc2IjYjUSdkc29sdmVGJA==
sol := dsolve({ODE,ICS}):
assign(sol);
Q := 'Re(evalf(q[i](t)))'$i=1..4:
plot( [Q[1],Q[2],Q[3]], t=0..30,title="Charge across each Cap",legend=["Q1","Q2","Q3"]);
VC := 'Re(evalf(q[i](t)))/C[i]' $i=1..3:
plot( [VC[1],VC[2],VC[3]], t=0..30,title="Potential across each Cap",legend=["C1","C2","C3"]);
IC := 'Re(evalf(diff(q[i](t),t)))' $i=1..3:
I0 := IC[1] + IC[2] + IC[3]:
plot( [I0, IC], t=0..15,title="Current across each Res.",legend=["I0","I1","I2","I3"]);
JSFH