y=(10^2*(1-((x-1000)/1000)^2))^0.5, in [0, 1000]