y=[ ];
for x=[0:100];
    y1=sqrt(22500-224*x+(7.2E9)*x+(-72E6)*x.^2);
    y=[y y1];
end
  plot([0:100],y)