Please use MATLAB code

3. (It is required to use a program to finish this question. Only programs written in MATLAB is acceptable. Hand in your code via through Ilearn (please put all of them in one ZIP file, in case you have more than one file). Print your tabulated results/graphs if the question asks you to. Please also print the code out and hand in the hard copy of the code.) Consider the following initial value problem: Sz'(t) = (x +t) 1.(0) = -1 on (0,1] (a) Approximate the solution r(t) by the modified Euler method with a mesh size h = 0.01. You only need to hand in the graph of the approximate solution. (b) Approximate the solution z(t) by the classical Fourth-order Runge-Kutta method (RK4) Wit1 = w; + (F1+2F2 +2F3 + Fy) where (F1 = hf(ti, wi) F2 = hf(ti + h, Wi+F1) F3 = hf (ti+th, W; + F2) F1 = hf(t; +h, Wi+F3) with a mesh size h = 0.01. Again, you only need to hand in the graph of the approximate solution