function [t,uT1,uT2,uT4]=prTaylor(N)
h=1/N;
t=[0:h:1];
uT1(1)=1;
uT2(1)=1;
uT4(1)=1;
for n=1:N
    uT1(n+1)=uT1(n)*(1-h)+t(n)*h+h;
    uT2(n+1)=uT2(n)*(1-h+h^2/2)+t(n)*(h-h^2/2)+h;
    uT4(n+1)=uT4(n)*(1-h+h^2/2-h^3/6+h^4/24)+t(n)*(h-h^2/2+h^3/6-h^4/24)+h;
end
