2. A common problem of applied mathematics is that of solving the integral equation f(x) + |K(x,t)y(t)dt -y(x), where the functions f(x) and K(x,t) are given and the problem is to compute y(x). If we approximate the integral by the quadrature formula then the integral equation becomes a system of linear algebraic equations: The solution y(C,),i -1,2,..., n, is the desired approximation to the function y(t). Using Simpson's rule, find an approximate solution of the integral equation 3 xy(t)dt =y(x) 0 (a) Write a MATLAB program to solve the equation above. The main program (script file testinteqn.m) calls a general integral equation solver inteqnsolve.m using the function y-inteqn(t, kernel, fun, coef) whose inputs are t - the vector of partition points, kernel- the kernel function K, fun - the function f(x), and coef- the vector of coefficients in the quadrature rule (Simpson's rule in this case) and whose output are the discrete solution values at t. Your program should be general enough to take any coefficient vector [a,,a,,..cn ]TTherefore you should first generate the coefficient vector for Simpson's rule before calling inteqn. You could use another function y = inteqnsolve (X, kernel, f) to generate the Simpson's rule coefficients and solve the integral equation by calling int eqn. In inteqn solve the inputs are x -the partition of [a, b], such that a = x(1)