Shiul's MATLAB solution added