Numerical algorithms with C