Qm=function(T){ # # Parameters Tmu=-1;Qmu=0.0042; # # Variable source if(T<=Tmu){ Qm=0; }else if(T>Tmu){ Qm=Qmu; } return(c(Qm)); }