Source: http://rims.ruforum.org/CABB640/fsolve-matlab-examples.pdf