#example 2 - x^2y" + 7xy' + 13y = 0
var('f,m')
f = m**2 + 6*m + 13
print f

var('m1,m2,y1,y2,c1,c2')
m1,m2 = solve(f==0, m)
print m1,m2

var('y,A,B')
y = (x**-3)*(A*cos(2*ln(x)) + B*sin(2*ln(x)))
print y

SageDays/Kreyszig-Example-4 (last edited 2010-12-16 11:13:34 by shantipriya)