h=150
H=10
y=[ ];
for x=[0:100];
   y1=sqrt(((H^2-h^2)/100)*x+h^2);
   y=[y y1];
end
   plot([0:100],y)