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

6 Ansichten (letzte 30 Tage)
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

Antworten (0)

Kategorien

Mehr zu Startup and Shutdown finden Sie in Help Center und File Exchange

Produkte


Version

R2018b

Community Treasure Hunt

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

Start Hunting!

Translated by