How to solve a differential equation with fractional order with function fde12
6 views (last 30 days)
Show older comments
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
1 Comment
Answers (0)
See Also
Categories
Find more on Startup and Shutdown in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!