How do I write the matlab code for implementing the Bjorck-Pereyra algorithm?

function: a=BjorckPereyra(x,f)

Input: The function f(x) is given by its values stored in the vector f at points stored in the vector x.

Both arrays are of the same size n.

Output: The program 'BjorckPereyra' computes the vector of the coefficients a in: p_n-1(x)=a1+a2x+a3x^2+...+anx^n-1

where P_n-1(x) is the interpolating polynomial.

Any help is greatly appreciated!