function [Y] = Spline1(P0, P1, T)

   t = [0.01:0.01:T];
   tau = (1/2) * (1 - cos(pi*t/T));
   Y = [];
   
   for i=1:length(t)
       a = tau(i);
       Y = [Y, (1-a)*P0 + a*P1];
   end
   
end