## 20 Reputation

0 years, 204 days

 >
 >

 (1)
 >
 >
 >

 (2)
 >

 (3)

 (4)

## @tomleslie       T...

Thank you for taking time out of your busy day.
Sincerely,

delvin

## @rcorless    ok...

The equation is correct, but I didn't pay attention, I wrote some in my own language,sorry.

## @Rouben Rostamian      S...

Sorry, I was not paying attention.

Thanks ocean.

## @tomleslie   Ok, I want to ch...

Ok, I want to check the different solutions of the equation.

Thanks a million, You made my day.

## @tomleslie    Thank osean...

Thank osean very much for your time.

"Because the answers are not what I want , The answers obtained are obvious."

m4.mw

## @Joe Riel  This is the code I ...

This is the code I wrote, but it didn't get a good answer

 >
 >

 (1)
 >
 (2)
 >
 (3)

 (4)
 >
 (5)

 (6)

 (7)

 (8)

 (9)

okey

## @dharr  thanks an ocean...

@dharr thanks an ocean

## How can I compare the values of A_0, A_1...

@dharr Thanks

 > restart; with(PDEtools):
 > pde2:=-lambda*diff(y(xi),xi)+y(xi)*diff(y(xi),xi)+diff(y(xi),xi\$2)=0;
 (1)
 >

 > phi(xi):=(p[1]*exp(q[1]*xi)+p[2]*exp(q[2]*xi))/(p[3]*exp(q[3]*xi)+p[4]*exp(q[4]*xi)); eval(phi(xi),[p[1]=I,p[2]=-I,p[3]=1,p[4]=1,q[1]=I,q[2]=-I,q[3]=I,q[4]=-I]); simplify(%); phi(xi):=%; U(xi):=A[0]+sum(A[d]*phi(xi)^d,d=1)+sum(B[d]*phi(xi)^(-d),d=1);
 (2)
 > subs(y(xi)=U(xi),pde2);
 (3)

## @Rouben Rostamian Thanks an ocean...

@Rouben Rostamian

Thanks an ocean

## Maple code that runs correctly...

@Rouben Rostamian

Dear Rouben Rostamian,

thank you for your explanation. You can see where this last part of the code is wrong. It is for balancing and obtaining the balance number.

## @Rouben Rostamian   hi,thank you ...

hi,thank you

```Thank you very much for your help
Where is this code wrong?

>
Error, `>` unexpected
restart;
with(PDEtools):
#local gamma:
get_order:=proc(expr,fun,var)
local the_list,i;
the_list := convert( indets(expr, function),list);
PDETools:-difforder(select( i-> has(i, diff(fun,var)),the_list))
end proc:
DGP:=proc(ff,var) local ee;
ee:=combine(frontend(expand,[ff]),':-power');
max(map(proc(u,v)
if u::`*` then
op(2,combine(select(type,u,
identical(v)^anything),
':-power'));
else
if type(u,identical(v)^anything) then
op(2,u);
else return FAIL end if;
end if;
end proc,
[`if`(ee::`+`,op(ee),ee)],var));
end proc:

fracdiff(u(x,t), x, alpha) assuming alpha > 0, alpha < 1:
FPDE := % + u(x,t)*diff(u(x,t),x) + diff(u(x,t),x,x) = 0;

/   /         (-alpha) /  d            \              \\
|   |(x - tau)         |----- u(tau, t)|              ||
|   |                  \ dtau          /              ||
|int|-----------------------------------, tau = 0 .. x||
\   \         GAMMA(1 - alpha)                        //

/ d         \   / d  / d         \\
+ u(x, t) |--- u(x, t)| + |--- |--- u(x, t)|| = 0
\ dx        /   \ dx \ dx        //
pde1 := D^alpha*u(x,y,t)+u(x,y,t)*diff(u(x,y,t),x)+diff(u(x,y,t),x\$2)=0;

tr1 := {u(x,t) = y(xi),  xi = x - v*t^alpha,  v = lambda*t^alpha / GAMMA(alpha+1)};

alpha                         / d            \
D      u(x, y, t) + u(x, y, t) |--- u(x, y, t)|
\ dx           /

/ d  / d            \\
+ |--- |--- u(x, y, t)|| = 0
\ dx \ dx           //
/             alpha                                      \
|     lambda t                     alpha                 |
< v = ----------------, xi = x - v t     , u(x, t) = y(xi) >
|    GAMMA(alpha + 1)                                    |
\                                                        /

pde2:=diff(y(xi), xi, xi) + y(xi)*diff(y(xi),xi) - lambda*diff(y(xi),xi) = 0;

/ d   / d        \\         / d        \          / d        \
|---- |---- y(xi)|| + y(xi) |---- y(xi)| - lambda |---- y(xi)| = 0
\ dxi \ dxi      //         \ dxi      /          \ dxi      /
n[1]:=get_order(pde2, U(xi),xi);
el1:=select(has, expand(lhs(pde2)), diff(U(xi), xi\$n[1]));
n[2]:=degree(expand(eval(lhs(pde2),diff=0)),U(xi));
el2:=select(has, expand(lhs(pde2)), U(xi)^(n[2]));
eq3:=expand(eval(el1,[U(xi)=a[h]*xi^(-h)]));
eq4:=expand(eval(el2,[U(xi)=a[h]*xi^h]));
d1:=DGP(eq3, xi) assuming h>0;
d2:=DGP(eq4, xi) assuming h>0;
assign(solve(-d1=d2,{h}));
print("N"=h);
0
Error, invalid input: diff expects 2 or more arguments, but received 1
-infinity
0
el1
0
FAIL
FAIL
"N" = h

```

 Page 1 of 1
﻿