#example 3 - 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

2.6-3 (last edited 2010-12-16 12:20:13 by shantipriya)