Numerical recipes in C