%Setpoint calculation for QEH
function sp=spQel(u)
%par_el=[p_el, mv0, Cp, kT]';
Tdes = u(1); 
p_el = u(2); 
mv0  = u(3);
Cp   = u(4); 
kT   = u(5);

sp   = Tdes - p_el*mv0*Cp/(2*kT);

end 