plot(r/(r+1)^2,r = 0 to 3)