function TT = evalLagrangeNewton (Lx,Ly,T) 
[U,DV] = DiffDiv(Lx,Ly); 
N = length(DV); [m,n] = size(T); 
TT = DV(1)*ones(m,n);  
for i = 2:N
    TT = TT.*(T -Lx(N-i+1))+DV(i); 
end
endfunction
