Cp=function(T){ # # Parameters Tml=-8;Tmu=-1;Cf=1.8;Cu=3.6;Ql=250; # # Variable heat capacity if(TTmu){ Cp=Cu; } return(c(Cp)); }