t=var('t')
Q=function('Q',t)
DE = 2*diff(diff(Q,t),t) + 8*diff(Q,t) + 10*Q 
de(t)=desolve (DE,[Q,t])
I1(t)=diff(de,t)          #  I(0)=0, Q(0)=0, but E(0)=10 , hence I'=E/L
k1,k2=var('k1,k2')
solve([de(0)==0,
       I1(0)==10/2],k1,k2)
de(t)=de(t).subs(k1 =5, k2 =0)
print 'I is' 
print de

SageDays/Kreyszig-2.12-11 (last edited 2010-08-11 12:33:45 by hardikgajera)