function u=SS(q2,c1,c2,yss) 
u = q2 * (yss-c2)/(c1-yss);
end