x=var('x')   #define a variable x
w=var('w')    #define a constant w
y=function('y',x)   #define y to be a function of that variable
t= diff(y,x)+(w/2)*cot(w*y)
desolve(t,[y,x],ics=(0,pi/4))

SageDays/Kreyszig-1.5-18 (last edited 2010-08-09 14:04:04 by rachitjha)