function z=phi(i,tri,xy)
% Quadratic triangular basis function associeted to vertex i
% Uses L, same parameters.
z=2*L(i,tri,xy).*(L(i,tri,xy)-1/2);
