R,L,C,E0,a,b,omega=var('R,L,C,E0,a,b,omega')
S=L*omega - 1/(omega*C)
I(t)=a*cos(omega*t)+b*sin(omega*t)
I(t)=I(t).subs(a=-(E0*S)/(R^2+S^2),b=(E0*R)/(R^2+S^2))
I(t)=I(t).subs(R=2,L=1,C=0.5,E0=50,omega=1)
print 'steady-state current is'
show(I(t))

SageDays/Kreyszig-2.12-4 (last edited 2010-08-10 12:15:09 by lokeshpimpale)