1. check problem 24: converges in C implementation but not in MATLAB
