FO logistic map in MATLAB

Matlab program for solution of the Caputo-like discrete fractional difference

Matlab program for solution of the Caputo-like discrete fractional difference Initial Value Problem

\Delta_*^q u(t)=f(t+q-1),u(t+q-1)), t\in N_{1-q}, u(0)=u_0

q\in(0,1), f continuous function depending on a real parameter p

The program can be simply adapted for other unidmensional map, or maps of higher dimensions

Example of use:

[u]=FO_discrete(0.8,2.3,.1,1000);

Regarding the paramer, see the ref.

The line

u(1)=u0+f(u0);

is designed to avoid the problem of the zero index in matlab