function [tau] = prp_time(eta)
% Schwarzschild radial spacelike geodesic
% Proper time in Schwarzschild space-time
% for radial spacelike geodesic with initial
% condition r_i = M.
% eta is a variable.
tau = (eta+sinh(eta)) *sqrt(1^3/8);