The code solves the Sachs equations in a similar fashion as the Dyer-Roeder approximation but
in addition it allows to account for the null shearing.
To use the code download sDR.f and
param.ini files.
To compile the code type "gfortran -o sDR sDR.f".