function [y1,T,y0] = exe2(t,y,jaco)
%
%  [alpha,T,y0] = exe2
% 
%  or
%             f = exe2(t,y)
% 
%  or 
%            f' = exe2(t,y,1)
% 
%  or
%             y = exe2(t)
%
%  Example 2 in  https://arxiv.org/abs/2407.11460
%
%  Taken from Example 2, pag. 51, in:
%
%  Z.Satmari. Iterative Bernstein spline technique applied to fractional 
%  order differential equations. Mathematical Foundations of Computing 6
%  (2023) 41-53. https://doi.org/10.3934/mfc.2021039
%
%
alfa = 3/2;
if nargin==0    % problem data
    y1   = alfa;
    T    = 1;
    y0   = [-1;0];
elseif nargin==1 % solution
    y1   = t.^(1.9)-1;
elseif nargin==2 % vector field
    y1   = y.^2/2 -(t.^(1.9)-1).^2/2 +(gamma(2.9)/gamma(1.4))*t.^(0.4);
%    y1   = y.^2/2 -(t.^(1.9)-1).^2/2 + ...
%         t.^(0.4)*(0.9*1.9*2^(4/5)*gamma(9/10)^2)/(gamma(9/5)*gamma(1/2));
else             % Jacobian
    y1   = y; 
end
return