function [y1,T,y0] = prob1(t,y,jaco)
%
% [alpha,T,y0] = prob1
%
% or
%            f = prob1(t,y)
%
% or
%           f' = prob1(t,y,1)
%
% or
%            y = prob1(t)
%
% Problem1 1 in:
%
%  L.Brugnano, G.Gurioli, F.Iavernaro, M.Vikerpuur. 
%  FDE-testset: comparing Matlab© codes for solving fractional
%  differential equations of Caputo type
%  Fractal Fract. 2025, 9(5), 312; https://doi.org/10.3390/fractalfract9050312
%
alfa = 0.5;
if nargin==0    % problem data
    y1   = alfa;
    T    = 1000;
    y0   = 10;
elseif nargin==1 % solution
    ybar = 1;
    y0   = 10;
    y1   = ml(-t.^alfa,alfa)*(y0-ybar)+ybar;
elseif nargin==2 % vector field
    y1   = 1-y;
else             % Jacobian
    y1   = -1;
end
return
