## 1481 Reputation

8 years, 193 days

## BVP method only work....

You have a BVP type ODE ,RK methods dosen't work(in Maple) for yours ICS and BC.

 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >
 >

## Try this: >&nbs...

`After eliminating  syntax errors:`

 >
 >
 (1)
 >
 >
 (2)
 >
 (3)
 >
 (4)
 >
 (5)
 >
 (6)
 >
 >

## Try this:...

Try this:

 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >
 (5)
 >
 >
 >
 (6)
 >
 (7)
 >
 > derivative
 >
 >

 >
 >
 (1)
 >
 (2)
 >
 >
 (3)
 >
 (4)
 >
 (5)
 >
 (6)
 >
 (7)
 >
 (8)
 >
 (9)
 >
 (10)
 >
 (11)
 >

## useint......

Command useint in dsolve helps:

 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >

## Another workround:...

```simplify(int(sin(y)^(1/3)*cos(y)^3, y = 0 .. x, AllSolutions));

#piecewise(sin(x) < 0, FAIL, 0 <= sin(x), -3*(2*sin(x)^2 - 5)*sin(x)^(4/3)/20)```

## Pity......

Yes, Maple can't  handle elliptic PDEs,because last Updates(Enhanced) to solving PDE's equation was in Maple 9 in year 2003.

During this time  no improvements have been made.

What a pity.

## pdsolve can solve......

pdsolve can solve,on Maple 2019.1 with Physics:-Version(399):

PDE := diff(u(x, t), t) + 2*u(x, t)^2*diff(u(x, t), x) - diff(u(x, t), x)^2 - 1/2*diff(u(x, t), x \$ 2)*u(x, t) = 0;

pdsolve([PDE, u(x, 0) = -tanh(x)]);

#u(x, t) = tanh(t - x)

## Use simplify command:...

sol := solve({sin(9*x - 1/3*Pi) = sin(7*x - 1/3*Pi)}, [x], explicit, allsolutions);

seq(simplify(sol[j]), j = 1 .. nops(sol));

#[x = 2*Pi*_Z13], [x = Pi*(2*_Z14 + 1)], [x = -1/48*Pi + 2*Pi*_Z15], [x = 47/48*Pi + 2*Pi*_Z15], [x = #23/48*Pi + 2*Pi*_Z15], [x = #-25/48*Pi + 2*Pi*_Z15], [x = 5/48*Pi + 2*Pi*_Z15], [x = -43/48*Pi + #2*Pi*_Z15], [x = -19/48*Pi + 2*Pi*_Z15], [x = 29/48*Pi + 2*Pi*_Z15], #[x = -7/48*Pi + 2*Pi*_Z15], [x = #41/48*Pi + 2*Pi*_Z15], [x = 17/48*Pi + 2*Pi*_Z15], [x = -31/48*Pi + 2*Pi*_Z15], [x = 11/48*Pi + #2*Pi*_Z15], [x = -37/48*Pi + 2*Pi*_Z15], [x = -13/48*Pi + 2*Pi*_Z15], [x = 35/48*Pi + 2*Pi*_Z15]

 >
 >
 (1)
 >

## With MMA...

With MMA I was able to solve  double integral.

 >
 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >
 (5)
 >
 >
 (6)
 >
 >
 >
 (7)
 >
 (8)
 >
 >

## Workaround....

It's seems Maple dosen't know the answer.With workaround:

 >
 (1)
 >
 >
 (2)
 >
 (3)
 >
 (4)
 >
 (5)
 >
 >
 (6)
 >
 (7)
 >
 >
 (8)
 >

## Try Explore...

Try:

Explore(plots:-implicitplot(y^2 = a*x^3 + b*x^2 + c*x + d, x = -5 .. 5, y = -5 .. 5, gridrefine = 2, view = [-5 .. 5, -5 .. 5]), a = -5 .. 5, b = -5 .. 5, c = -5 .. 5, d = -5 .. 5);

## Try:...

```func := 2*omega^2*B[1]*Zeta*omega[0] - omega^3*B[1]*I + omega*B[1]*omega[0]^2*I - 2*I*B[0]*Zeta*omega*omega[0] - B[0]*omega^2 + B[0]*omega[0]^2;

(evalc(func) assuming (B[0] in real, omega[0] in real, omega in real, Zeta in real, B[1] in real));

#2*omega^2*B[1]*Zeta*omega[0] - B[0]*omega^2 + B[0]*omega[0]^2 + (-2*Zeta*omega*B[0]*omega[0] - omega^3*B[1] + omega*B[1]*omega[0]^2)*I```

## You can find x in terms of y(x) using so...

You can find x in terms of y(x) using fsolve command.Analytical solution is  impossible.

 >
 >
 (1)
 >
 >
 >
 (2)
 >
 (3)
 >
 >
 >
 >
 >
 >