There is no closed-form solution for the error function: erf(a) 2 ra e dx Numerically evaluate erf(1.5) using Romberg integration to an accuracy of Es 0.1% (i.e. until Ea Es). Present your results in the format of Fig. 20.1. You may not use the romberg M-file provided by the textbook (I want you to become familiar with the Romberg formulas and their application), but you are welcome to use the trap M-file or some other aid to speed up the trap rule calculations. Also, compute the true percent relative error based on the true value, which can be determined using MATLAB's built-in function erf. There is no closed-form solution for the error function: erf(a) 2 ra e dx Numerically evaluate erf(1.5) using Romberg integration to an accuracy of Es 0.1% (i.e. until Ea Es). Present your results in the format of Fig. 20.1. You may not use the romberg M-file provided by the textbook (I want you to become familiar with the Romberg formulas and their application), but you are welcome to use the trap M-file or some other aid to speed up the trap rule calculations. Also, compute the true percent relative error based on the true value, which can be determined using MATLAB's built-in function erf