Question: How I can pdsolve this partial fractional equation?

How I can pdsolve this partial fractional  equation?

1.mw
 

restart; CK := .3; Z := 10; L := 1; alpha := .95

###The Riemann–Liouville fractional integral H^alpha  is written in a convolution-form as:

Parse:-ConvertTo1D, "first argument to _Inert_ASSIGN must be assignable"

 

#where GAMMA(...) represents the Gamma function

 

 
PDE := diff(theta(xi, beta), beta, beta)+L*(diff(theta(xi, beta), xi, beta, beta))+diff(theta(xi, beta), beta, beta, beta)+(1/2)*(diff(theta(xi, beta), beta, beta, beta, beta)) = H^(alpha-1)*CK*(diff(theta(xi, beta), xi, xi))+H(CK*Z+1)^(alpha-1)*(diff(theta(xi, beta), xi, xi, beta))+H^(alpha-1)*Z*(diff(theta(xi, beta), xi, xi, beta, beta))

Init := {theta(xi, 0) = 0, (D[2](theta))(xi, 0) = 0}

{theta(xi, 0) = 0, (D[2](theta))(xi, 0) = 0}

(1)

Bdry := {theta(0, beta) = 1, theta(10, beta) = 0}

{theta(0, beta) = 1, theta(10, beta) = 0}

(2)

smod1 := pdsolve(PDE, Init, Bdry, numeric)

 

 

``


 

Download 1.mw

 

Please Wait...