function qin=Qin(h)
if 0<=h && h<=20
    qin=1-h/20;
else
    qin=0;
end