Numeric routines for optically stimulated luminescence dating