function [tau] = tm_prp_t(eta)
% Schwarzschild radial time-like geodesic
% Proper time in Schwarzschild space-time
% for radial time-like geodesic with initial
% condition r_i = 6M. (M = 1)
% eta is a variable.
tau = (eta+sin(eta)) *3*sqrt(3);