simpson39;s 13 rule matlab code