\% 6: What's happening on line 90? mytri = DelaunayTri(xobs,yobs); XY=mytri.X; [NXY, i] = size(XY); TRI=mytri.Triangulation; [NTRI, i] = size(TRI); for i = [1:NTRI] v1=TRI(i,1); v2=TRI(i,2); v3=TRI(i,3); plot ([XY(v1,2),XY(v2,2)],[XY(v1,1),XY(v2,1)], 'k-', 'LineWidth', 2); plot ([XY(v2,2),XY(v3,2)],[XY(v2,1),XY(v3,1)], 'k-', 'LineWidth', 2); plot ([XY(v3,2),XY(v1,2)],[XY(v3,1),XY(v1,1)],k, 'LineWidth', 2); end