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=8,L=2,C=0.1,E0=160,omega=5)
print 'steady-state current is'
show(I(t))

SageDays/Kreyszig-2.12-5 (last edited 2010-08-10 12:10:47 by lokeshpimpale)