How to solve a differential equation with fractional order with function fde12

6 visualizzazioni (ultimi 30 giorni)
Hi I am trying to solve this fractional order differential equation D^(1/5)y(t)+7*D^(1/10)y(t)=100 with using Caputo Derivative, when I tried to use the fuction Fde12 there is an error, Maybe you can help me I am new on this,
% Enfoque de caputo
%Resolución de ec. de orden no entero usando fde12
f=@(t,y)[y(2);100-7*y(2)];
h=0.00001;
gam=1/10;
t0=0.001;
tf=10;
y0=[0;0];
[t,y]=fde12(gam,f,t0,tf,y0,h);
plot(t,y(1,:)),grid,title('Función fde12')
Thanks for your help

Risposte (0)

Categorie

Scopri di più su Startup and Shutdown in Help Center e File Exchange

Prodotti


Release

R2018b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by