for u=-9999:1:9999,
x(u+10000)=450*(1+u-u*sign(u))/(1-u*sign(u));
end
plot(x,u)