About Neville's Algorithm

So how exactly can I number my functions if Matlab doesn't support function numbering?

For example:

I got n nodes with x=-1:(2/n):1

But since there's no function numbering, I can't create P0,1(x)~Pn-1,n(x) with loops

The seemingly doable way I can come up now is to define function p=((input point)-x(i)) and go for a double loop, which avoids numbering the functions in the first place