Numerical Recipes in fortran 90