3.3 How to find a particular solution to ODE?
restart;
ode:=diff(y(x),x)+y(x)^2*sin(x)-2*sin(x)/cos(x)^2 = 0;
yp:=DETools:-particularsol(ode);
To step into the code, do
restart;
ode:=diff(y(x),x)+y(x)^2*sin(x)-2*sin(x)/cos(x)^2 = 0;
stopat(`DEtools/particularsol`);
DETools:-particularsol(ode);
To print it do
print(`DEtools/particularsol`);