t=var('t')
Y=function('Y',t)
DE = 10*diff(diff(Y,t),t) + 10*Y-10 
d(t)=desolve (DE,[Y,t])
d1(t)=diff(d(t),t)
k1,k2=var('k1,k2')
solve([d(0)==0,
       d1(0)==0],k1,k2)
d=d(t).subs(k1 = 0, k2 = -1)
print d

SageDays/Kreyszig-2.12-13 (last edited 2010-08-12 08:34:46 by hardikgajera)