clc;
clear;
close?all;

rand（‘state‘?0）;
node?=?8;
x?=?rand（1node）;
y?=?rand（1node）;
%?delaunay是MATLAB中三角剖分的函數，返回的TRI是三角形的矩陣
%?TRI的每一行表示三角形的三個點
TRI?=?delaunay（xy）;

%?繪圖
figure;?
xmin?=?min（x（:））;?xmax?=?max（x（:））;
ymin?=?min（y（:））;?ymax?=?max（y（:））;
xl?=?xmax?-?xmin;?yl?=?ymax?-?ymin;
axis（[xmin-xl*0.1?xmax+xl*0.1...
????ymin-yl*0.1?ymax+yl*0.1]）;
hold?on;

n?=?size（TRI?1）;
for?i?=?1?:?n
????t1?=?TRI（i?:）;
????for?j?=?1?:?length（t1）-1
????????xt?=?[x（t1（j））?x（t1（j+1））];
????????yt?=?[y（t1（j））?y（t1（j+1））];
??????