Items tagged with fractional



I hope every one is ok.

I am running this code (see below)

m := 2;

X[0] := 14;
Y[0] := 18;
a := 1; b := 1; c := .1; d := 1;

alpha := 1;

for k from 0 to m do X[k+1] := GAMMA(k*alpha+1)*(a*X(k)-b*(sum(X(s)*Y(k-s), s = 0 .. k)))/GAMMA(k*alpha+1+1); Y[k+1] := GAMMA(k*alpha+1)*(-c*Y(k)+d*(sum(X(s)*Y(k-s), s = 0 .. k)))/GAMMA(k*alpha+1+1) end do

x := 0; y := 0

The following message pop out.

PLease HELP! HELP!.....



I have recently been working on a problem using fractional calculus and have come across something in Maple's fracdiff  command that makes no sense to me.

fracdiff(1, x, 1/2) = 0

It should be:     1/(sqrt(x)*sqrt(Pi))




I want to draw  phase plane of system of three fractional order equations. 


Note that 

Also want the  phase portrait when the values of alpha are not same....


 i can pdsolve this equation numerically or analyticlly?

this equation is time-fractional  equation with generalized Cattaneo model



 is the fractional derivative operator considered in the
Caputo sense.


k := 1; -1; rho := 1; -1; h := 1; -1; alpha := 2-Upsilon; -1; 0 < Upsilon and Upsilon <= 1

0 < Upsilon and Upsilon <= 1


k*(diff(T(z, t), z, z)) = rho*(diff(T(z, t), [`$`(t, alpha)]))

diff(diff(T(z, t), z), z) = diff(T(z, t), [`$`(t, 2-Upsilon)])


k*(diff(T((1/2)*h, t), z)) = 1:

k*(diff(T((-h)*(1/2), t), z)) = 0:

T(z, 0) = 0

T(z, 0) = 0





i want a scheme of fractional differential equation so that i solve my questions and make a code of it.

please provide me the scheme


I am simulate the code for fractional differential equation. But the out put is not wright...


S[0] := .8;



V[0] := .2;



R[0] := 0;



alpha := 1;







gamma = 0.3e-1


q := .9;



T := 1;



N := 5;



h := T/N;




for i from 0 to N do for j from 0 to 0 do a[j, i+1] := i^(alpha+1)-(i-alpha)*(i+1)^alpha; b[j, i+1] := h^alpha*((i+1-j)^alpha-(i-j)^alpha)/alpha end do end do;

for n from 0 to N do Sp[n+1] = S[0]+(sum(b[d, n+1]*(mu*(1-q)-beta*S[d]*V[d]-mu*S[d]), d = 0 .. n))/GAMMA(alpha); Vp[n+1] = V[0]+(sum(b[d, n+1]*(beta*S[d]*V[d]-(mu+gamma)*S[d]), d = 0 .. n))/GAMMA(alpha); Rp[n+1] = R[0]+(sum(b[d, n+1]*(mu*q-mu*R[d]+gamma*V[d]), d = 0 .. n))/GAMMA(alpha); S[n+1] = S[0]+h^alpha*(mu*(1-q)-beta*Sp[n+1]*Vp[n+1]-mu*Sp[n+1])/GAMMA(alpha+2)+h^alpha*(sum(a[e, n+1]*(mu*(1-q)-beta*S[e]*V[e]-mu*S[e]), e = 0 .. n))/GAMMA(alpha+2); V[n+1] = V[0]+h^alpha*(beta*Sp[n+1]*Vp[n+1]-(mu+gamma)*Sp[n+1])/GAMMA(alpha+2)+h^alpha*(sum(a[e, n+1]*(beta*S[e]*V[e]-(mu+gamma)*S[e]), e = 0 .. n))/GAMMA(alpha+2); R[n+1] = R[0]+h^alpha*(mu*q-mu*Rp[n+1]-gamma*Vp[n+1])/GAMMA(alpha+2)+h^alpha*(sum(a[e, n+1]*(mu*q-mu*R[e]-gamma*V[e]), e = 0 .. n))/GAMMA(alpha+2) end do;

Sp[1] = .7184000000


Vp[1] = 0.692454936e-1


Rp[1] = 0.9508862660e-1


S[1] = .7632000000-0.8000000000e-1*Sp[1]*Vp[1]-0.4000000000e-1*Sp[1]


V[1] = .1346227468+0.8000000000e-1*Sp[1]*Vp[1]-(1/10)*(.4+gamma)*Sp[1]


R[1] = 0.6045568670e-1-(1/10)*gamma*Vp[1]-0.4000000000e-1*Rp[1]


Sp[2] = .7264000000-.1600000000*S[1]*V[1]-0.8000000000e-1*S[1]


Vp[2] = 0.692454936e-1+.1600000000*S[1]*V[1]-.1954431330*S[1]


Rp[2] = .1670886266+.1154431330*V[1]-0.8000000000e-1*R[1]


S[2] = .7712000000-0.8000000000e-1*Sp[2]*Vp[2]-0.4000000000e-1*Sp[2]-.1600000000*S[1]*V[1]-0.8000000000e-1*S[1]


V[2] = .1346227468+0.8000000000e-1*Sp[2]*Vp[2]-(1/10)*(.4+gamma)*Sp[2]+.1600000000*S[1]*V[1]-.1954431330*S[1]


R[2] = .1324556867-(1/10)*gamma*Vp[2]-0.4000000000e-1*Rp[2]-.1154431330*V[1]-0.8000000000e-1*R[1]


Sp[3] = .7344000000-.1600000000*S[1]*V[1]-0.8000000000e-1*S[1]-.1600000000*S[2]*V[2]-0.8000000000e-1*S[2]


Vp[3] = 0.692454936e-1+.1600000000*S[1]*V[1]-.1954431330*S[1]+.1600000000*S[2]*V[2]-.1954431330*S[2]


Rp[3] = .2390886266+.1154431330*V[1]-0.8000000000e-1*R[1]+.1154431330*V[2]-0.8000000000e-1*R[2]


S[3] = .7792000000-0.8000000000e-1*Sp[3]*Vp[3]-0.4000000000e-1*Sp[3]-.1600000000*S[1]*V[1]-0.8000000000e-1*S[1]-.1600000000*S[2]*V[2]-0.8000000000e-1*S[2]


V[3] = .1346227468+0.8000000000e-1*Sp[3]*Vp[3]-(1/10)*(.4+gamma)*Sp[3]+.1600000000*S[1]*V[1]-.1954431330*S[1]+.1600000000*S[2]*V[2]-.1954431330*S[2]


R[3] = .2044556867-(1/10)*gamma*Vp[3]-0.4000000000e-1*Rp[3]-.1154431330*V[1]-0.8000000000e-1*R[1]-.1154431330*V[2]-0.8000000000e-1*R[2]


Sp[4] = .7424000000-.1600000000*S[1]*V[1]-0.8000000000e-1*S[1]-.1600000000*S[2]*V[2]-0.8000000000e-1*S[2]-.1600000000*S[3]*V[3]-0.8000000000e-1*S[3]


Vp[4] = 0.692454936e-1+.1600000000*S[1]*V[1]-.1954431330*S[1]+.1600000000*S[2]*V[2]-.1954431330*S[2]+.1600000000*S[3]*V[3]-.1954431330*S[3]


Rp[4] = .3110886266+.1154431330*V[1]-0.8000000000e-1*R[1]+.1154431330*V[2]-0.8000000000e-1*R[2]+.1154431330*V[3]-0.8000000000e-1*R[3]


S[4] = .7872000000-0.8000000000e-1*Sp[4]*Vp[4]-0.4000000000e-1*Sp[4]-.1600000000*S[1]*V[1]-0.8000000000e-1*S[1]-.1600000000*S[2]*V[2]-0.8000000000e-1*S[2]-.1600000000*S[3]*V[3]-0.8000000000e-1*S[3]


V[4] = .1346227468+0.8000000000e-1*Sp[4]*Vp[4]-(1/10)*(.4+gamma)*Sp[4]+.1600000000*S[1]*V[1]-.1954431330*S[1]+.1600000000*S[2]*V[2]-.1954431330*S[2]+.1600000000*S[3]*V[3]-.1954431330*S[3]


R[4] = .2764556867-(1/10)*gamma*Vp[4]-0.4000000000e-1*Rp[4]-.1154431330*V[1]-0.8000000000e-1*R[1]-.1154431330*V[2]-0.8000000000e-1*R[2]-.1154431330*V[3]-0.8000000000e-1*R[3]


Sp[5] = .7504000000-.1600000000*S[1]*V[1]-0.8000000000e-1*S[1]-.1600000000*S[2]*V[2]-0.8000000000e-1*S[2]-.1600000000*S[3]*V[3]-0.8000000000e-1*S[3]-.1600000000*S[4]*V[4]-0.8000000000e-1*S[4]


Vp[5] = 0.692454936e-1+.1600000000*S[1]*V[1]-.1954431330*S[1]+.1600000000*S[2]*V[2]-.1954431330*S[2]+.1600000000*S[3]*V[3]-.1954431330*S[3]+.1600000000*S[4]*V[4]-.1954431330*S[4]


Rp[5] = .3830886266+.1154431330*V[1]-0.8000000000e-1*R[1]+.1154431330*V[2]-0.8000000000e-1*R[2]+.1154431330*V[3]-0.8000000000e-1*R[3]+.1154431330*V[4]-0.8000000000e-1*R[4]


S[5] = .7952000000-0.8000000000e-1*Sp[5]*Vp[5]-0.4000000000e-1*Sp[5]-.1600000000*S[1]*V[1]-0.8000000000e-1*S[1]-.1600000000*S[2]*V[2]-0.8000000000e-1*S[2]-.1600000000*S[3]*V[3]-0.8000000000e-1*S[3]-.1600000000*S[4]*V[4]-0.8000000000e-1*S[4]


V[5] = .1346227468+0.8000000000e-1*Sp[5]*Vp[5]-(1/10)*(.4+gamma)*Sp[5]+.1600000000*S[1]*V[1]-.1954431330*S[1]+.1600000000*S[2]*V[2]-.1954431330*S[2]+.1600000000*S[3]*V[3]-.1954431330*S[3]+.1600000000*S[4]*V[4]-.1954431330*S[4]


R[5] = .3484556867-(1/10)*gamma*Vp[5]-0.4000000000e-1*Rp[5]-.1154431330*V[1]-0.8000000000e-1*R[1]-.1154431330*V[2]-0.8000000000e-1*R[2]-.1154431330*V[3]-0.8000000000e-1*R[3]-.1154431330*V[4]-0.8000000000e-1*R[4]


Sp[6] = .7584000000-.1600000000*S[1]*V[1]-0.8000000000e-1*S[1]-.1600000000*S[2]*V[2]-0.8000000000e-1*S[2]-.1600000000*S[3]*V[3]-0.8000000000e-1*S[3]-.1600000000*S[4]*V[4]-0.8000000000e-1*S[4]-.1600000000*S[5]*V[5]-0.8000000000e-1*S[5]


Vp[6] = 0.692454936e-1+.1600000000*S[1]*V[1]-.1954431330*S[1]+.1600000000*S[2]*V[2]-.1954431330*S[2]+.1600000000*S[3]*V[3]-.1954431330*S[3]+.1600000000*S[4]*V[4]-.1954431330*S[4]+.1600000000*S[5]*V[5]-.1954431330*S[5]


Rp[6] = .4550886266+.1154431330*V[1]-0.8000000000e-1*R[1]+.1154431330*V[2]-0.8000000000e-1*R[2]+.1154431330*V[3]-0.8000000000e-1*R[3]+.1154431330*V[4]-0.8000000000e-1*R[4]+.1154431330*V[5]-0.8000000000e-1*R[5]


S[6] = .8032000000-0.8000000000e-1*Sp[6]*Vp[6]-0.4000000000e-1*Sp[6]-.1600000000*S[1]*V[1]-0.8000000000e-1*S[1]-.1600000000*S[2]*V[2]-0.8000000000e-1*S[2]-.1600000000*S[3]*V[3]-0.8000000000e-1*S[3]-.1600000000*S[4]*V[4]-0.8000000000e-1*S[4]-.1600000000*S[5]*V[5]-0.8000000000e-1*S[5]


V[6] = .1346227468+0.8000000000e-1*Sp[6]*Vp[6]-(1/10)*(.4+gamma)*Sp[6]+.1600000000*S[1]*V[1]-.1954431330*S[1]+.1600000000*S[2]*V[2]-.1954431330*S[2]+.1600000000*S[3]*V[3]-.1954431330*S[3]+.1600000000*S[4]*V[4]-.1954431330*S[4]+.1600000000*S[5]*V[5]-.1954431330*S[5]


R[6] = .4204556867-(1/10)*gamma*Vp[6]-0.4000000000e-1*Rp[6]-.1154431330*V[1]-0.8000000000e-1*R[1]-.1154431330*V[2]-0.8000000000e-1*R[2]-.1154431330*V[3]-0.8000000000e-1*R[3]-.1154431330*V[4]-0.8000000000e-1*R[4]-.1154431330*V[5]-0.8000000000e-1*R[5]





Download sir_(2).mw


How to find the determining equation for a system of fractional differential equation using Maple 15?


I have been looking at some new models of Casio Scientific Calculators and came across with "Fx-115es Plus" Model which seem to have a some sort of simple CAS(Computer Algebra System) built into it.

Two new features which i really liked were

(i) Ability to make any part of the expression inert and simplying the rest.

(ii) Fully Integrated Repeated decimal display for fractions.


I want to ask if there is any builtin commands that can achieve these two effects in maple.

I will give some example for each of these

(i) simplifying say 2^3*2^4 in maple gives 32.

but forexample if i want to make 2 in the bases inert then simplifying the result should give 2^7

if i make 3 inert then the result is 16*2^3

if i make 4 inert then the result is 8*2^4

another example say (2^3)^4 in maple gives 4096

but if i make 2 inert then the result should be 2^12

if i make 3 inert then the result is 16^3

if i make 4 inert then the result is 8^4

In this way it is possible to keep any interesting part of large complex expression unevaluated and simplifying the rest across it to maintain focus on the interesting part.

I know i can try to achieve this effect by using unevaluation quotes but they get messy and harder to track in large nested forms.

Another approach might be to replace the inert parts by explicit undeclared symbols with required assumptions and simplifying, but this is not it.

I know in Maple 18 they have introduced some package called InertForm or something, can it achieve this effect and also mark inert parts of the expression as grey like it is possible for some operators.

(ii) the example for the second is quite obvious, say given the fraction 237/14, evalf of this gives 16.92857143 but a result like 16.9Overscript[285714, _] is more closer to differentiation it from a irrational expansion. Sorry i donot know how to pretty print this here.

Another advantage is when i want to give some large repeating decimal expansion and have maple convert it to fractional form. Currently i have no idea how many times to repeat the decimals explicitly to make maple understand that it is a repeating decimal expansion.

When Benoit Mandelbrot was still alive I sent him an email:

I am writing to you because I have trouble understanding why
the covariance function in the Fractional Brownian Motion (FBM) is given by:


I would be very grateful if you could please explain this to me in simple and step by step terms.
I understand simple ARIMA models (P(t)=P(t-1)+E(t) where E(t)=p*E(t-1)+r(t)) which have the

Page 1 of 1