Chapter 5.4, Problem 9CP

Write a MATLAB function called myerf.m that uses Adaptive Simpsons Rule to calculate the value of

$erf\left(x\right)=\frac{2}{\sqrt{\pi }}{\displaystyle {\int }_0^x{e}^{-s^2}ds}$

within eight correct decimal places for arbitrary input $x$ . Test your program for

$x=1$ and

$x=3$ by comparing with MATLAB’s function erf.