function [ri] =  i_radius(xi, l, e)
% Schwarzschild spacelike geodesic
% Radius in Schwarzschild space-time
% for spacelike geodesic with immaginary eccentricity
% l and e are parameters
% xi is a variable.
ri = l./(1+e*tan(xi/2));