Two-FORTRAN 77 function subprograms to calculate gravity anomalies of bodies of finite and infinite strike length with the density contrast differing with depth