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".