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

6 次查看(过去 30 天)
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

回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Startup and Shutdown 的更多信息

产品


版本

R2018b

Community Treasure Hunt

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

Start Hunting!

Translated by