function [r1] =  radius1(chi, l, e)
% Schwarzschild spacelike geodesic
% Radius in Schwarzschild space-time
% for spacelike geodesic of the first kind
% l and e are parameters
% chi is a variable.
r1 = l./(1+e*cos(chi));