Exercise and solutions from Ogata's book in MATLAB