The ODE diff(y(x),x) = sec(x)^2*sec(y(x))^3  can be solved as separable. So the answer should be 

simple_answer:=sin(y(x))*(cos(y(x))^2+2)=C_1+3*tan(x);

as can be seen by direct integration of each side of the differential equation. I am trying to make Maple give the same answer, but not having any luck. 

restart;
ode:=diff(y(x),x) = sec(x)^2*sec(y(x))^3;
sol:=dsolve(ode, y(x),implicit);

I tried simplify(sol,trig) and tried simplify(sol,size).  Both Maple answer, and the simple answer solve the ODE.

Is there a way to make Maple dsolve give the simpler answer, or simplify/convert the answer it gives to the simpler one? I am newbie in Maple.

Hello; I found another ODE which Maple gives division by zero on.  Is this also a bug? 

dsolve(x*(a^2*x+(x^2-y(x)^2)*y(x))*diff(y(x),x)^2-(2*a^2*x*y(x)-(x^2-y(x)^2)^2)*diff(y(x),x)+a^2*y(x)^2-x*y(x)*(x^2-y(x)^2) = 0, y(x));

Error, (in dsolve) numeric exception: division by zero

This is from a book. Using Maple 2016.1 on windows.

Maple 2016.1 on windows. This ODE from a book, and Maple gives division by zero. Is this a bug or expected?

ode:=(a^2*x+(x^2-y(x)^2)*y(x))*diff(y(x),x)+x*(x^2-y(x)^2) = a^2*y(x);
dsolve(ode, y(x));

Error, (in dsolve) numeric exception: division by zero

Mathematica gives this to same ODE, but no division by zero.

DSolve[x*(x^2 - y[x]^2) + (a^2*x + y[x]*(x^2 - y[x]^2))*y'[x] == a^2*y[x], y[x], x]

Solve[x^2/2 - (1/2)*a^2*Log[x - y[x]] + (1/2)*a^2*Log[x + y[x]] +y[x]^2/2 == C[1], y[x]]

Where is the division by zero coming from in Maple?

2016.1 on windows.

I am learning how to use Maple with boundary value ODE. Given this ODE

y''''(x)+ lam* y(x) =0

with some B.C., say  y(0)=0,y'(0)=0,y''(L)=0,y'''(L)=0, where L is length.

I can't figure the correct syntax to use. It seems Maple do not like the syntax I am using, but it works on a second order ODE?

Here is my attempt:

restart;
assume(lam>0); assume(L>0);
bc:=y(0)=0,D[1](y)(0)=0,D[2](y)(L)=0,D[3](y)(L)=0;
dsolve({diff(y(x),x$4)+lam*y(x)=0,bc},y(x));

Error is 

Error, (in evalapply) too few variables for the derivative with respect to the 2nd variable
Error, (in dsolve) found the following equations not depending on the unknowns of the input system: {bc}

But on a simpler second order ODE, the syntax works

restart;
assume(lam>0); assume(L>0);
bc:=y(0)=0,D[1](y)(0)=0;
dsolve([diff(y(x),x$2)+lam*y(x)=0,bc],y(x));

No error. 

Is the syntax I am using in first example wrong? what would be the correct syntax? I googled for long time, and can't find one example that shows how to use BVP with higher order ODE. I am Maple newbie.

How does one obtain all solutions from dsolve? I see an option called Allsolutions, but this seems to only apply to solve and other functions. It does not work with dsolve.

For example, maple gives one solution for the following first order non-linear ODE. But the ODE has another solution y(x)=0 as well. How does one tell Maple to return all solutions? I am interested in this when using the 'implicit' option mainly. Here is an example

restart;
num:=-(exp(x)*sin(y(x))-2*y(x)*sin(x)):
den:=(exp(x)*cos(y(x))+2*cos(x)):
eq:=diff(y(x),x)=num/den;

r0:=dsolve(eq,y(x));

But when I tried y(x)=0, it turned out to also be a solution

odetest(y(x)=0,eq);
                          0

But dsolve did not return this solution on its own along with the first one.  But on another example, Maple did well, and returned all solutions. Here is the other example

eq:=(2*x*y(x)^2+2*y(x))+(2*x^2*y(x)+2*x)*diff(y(x),x);
dsolve(eq=0,y(x),'implicit');

In the above, Maple returned the two solutions. 

Is there a correct way to tell Maple dsolve to return all solutions all the time? Why did it return both solution in the above example, but not in the first example?

I am maple newbie. Thank you.