plot (((1.01815748 * x) / (0.23133484 * x^2 + 1.)) - sin(x))/sin(x) for x in [0, pi/2]