这个函数计算具有网格点x和网格值y的拉格朗日多项式在给定点的值。这个功能基于重心法实现,并且经过了矢量化优化以提高运算速度。网格点x是一个行向量,网格值y是一个矩阵,每一行代表一个多项式在不同网格点处的值。返回的矩阵p中,p(i,j)表示第i个多项式在点t(j)处的值。