x,c=var('x,c')
y = function('y', x)
y=c*exp(2*x)
d_y=derivative(y,x)
# d_y is the first derivative of y wrt x
y = var('y')
c_val=solve([c==y*exp(-2*x)], c)

d_y=d_y.subs(c_val[0])
print d_y

SageDays/Kreyszig-1.8-5 (last edited 2010-08-12 10:06:08 by msunyuma)