数值计算二维拉格朗日插值的求解

怎么求二维的拉格朗日插值

很简单代入已有的公式求，公式在课件里面给了

function val = linquadref(xi,eta)
val=[1/4*(1-xi)*(1-eta);1/4*(1+xi)*(1-eta);1/4*(1+xi)*(1+eta);1/4*(1-xi)*(1+eta)];
end

function deriv = linquadderivref(xi,eta)
deriv=[-1/4*(1-eta),-1/4*(1-xi)
        1/4*(1-eta),-1/4*(1+xi)
        1/4*(1+eta),1/4*(1+xi)
       -1/4*(1+eta),1/4*(1-xi)];
end

f=[0 1 3 1];
wert1=f*linquadref(0,0)
wert2=f*linquadref(0.577,-0.577)
dwert1=f*linquadderivref(0,0)
dwert2=f*linquadderivref(0.577,-0.577)