Implement the procedures coef and eval shown on above picture in C++, NOT MatLab please. Then find the Newton form of the interpolating polynomial p for sin(x) at 10 equidistant points in the interval [0,1.6875]. If only the pseudocode can be done, that will still be very helpful.

If

Coef Pseudocode procedure Coef (n, (xi), (yi), (a;)) integer i, j, n; real array (Xi)0:n, (yi)0:n, (aj)o:n for i = 0 ton a; y end for for j = 1 ton for i = n to j step-1 a; + (a; 4;-1)/(x; xi-;) end for end for end procedure Coef The second is function Eval. It requires as input the array (xi) from the original table and the array (a;), which is output from Coef. The array (a;) contains the coefficients for the Newton form of the interpolation polynomial. Finally, as input, a single real value for t is given. The function then returns the value of the interpolating polynomial at t. Eval Pseudocode real function Evalin, (x;), (a;), t) integer i, n; real t, temp; real array (x;)o:n, (a;)o:n temp + an for i = n - 1 to 0 step -1 temp + (temp)(t - x;) +a; end for Eval temp end function Eval