How is the Riemann-Liouville fractional order integral code in MATLAB?

How is the Riemann-Liouville fractional order integral code in MATLAB?

[iraj mohammadpour]

Riemann-Liouville fractional order integral formula

Answer 1

[Walter Roberson]

syms alpha t s real
syms f(s)
Iaf(t) = 1/gamma(alpha) * int((t-s)^(alpha-1) * f, s, 0, t)

Iaf(t) = 

%example
subs(Iaf, {f}, {sin(s)})

ans(t) = 

subs(Iaf, {alpha, f}, {sym(5)/4, sin(s)})

ans(t) = 

%example 2
subs(Iaf, {f}, {s*(s-5)})

ans(t) = 

subs(Iaf, {alpha, f}, {sym(5)/4, s*(s-5)})

ans(t) =