t = var('t')
x = function('x',t)
DE = diff(x, t, 2) - diff(x,t)
m = desolve(DE, [x,t])
print m
a = m.subs(t=0)
#print a
b = diff(m, t).subs(t=0)
#print b
k1 = var('k1')
b = function('b',k1)
k1 = solve(b==2,k1)
#print k1
a.subs(k1=2)
#print a
k2 = var('k2')
a = function('a',k2)
k2 = solve(a==2,k2)
#print k2
x = m.subs(k1=2,k2=0,t=6)
v = diff(m,t).subs(k1 =2, k2=0, t=6)
print x
print v

SageDays/Kreyszig-2.1-14 (last edited 2010-08-09 14:11:52 by zubinmehta)