function [y,Temp]=Gp(u,h,Temp)
T=10;K=2;Tao=3;
Ytemp=Temp(Tao/h);
K1=(-Ytemp+K*u)/T;
K2=(-(Ytemp+0.5*K1*h)+K*u)/T;
K3=(-(Ytemp+0.5*K2*h)+K*u)/T; 
K4=(-(Ytemp+K3*h)+K*u)/T;
 Ytemp=Ytemp+h*(K1+2*K2+2*K3+K4)/6;
 y=Temp(1);
for i=1:Tao/h-1
    Temp(i)=Temp(i+1);
end
Temp(Tao/h)=Ytemp;