x=var('x')   #define a variable x
y=function('y',x)   #define y to be a function of that variable
t= diff(y,x)-(2*y)*cos(2*x)/sin(2*x)
desolve(t,[y,x],ics=(pi/4,3.8))

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