Do in MATLAB or Maple

$7\mathrm{.}1$ Vandermonde$-$based collocation

Consider the function $f(x)=\frac{30}{1+3x^2},$ defined over the interval $-3,+3.$ Using the $7$

data points

$(-3,f(-3)),(-2,f(-2)),(-1,f(-1)),(0,f(0)),(1,f(1)),(2,f(2)),(3,f(3)),$

compute the corresponding interpolation polynomial of degree $6,$ using Vandermonde$-$

based collocation.