clear;
k=0;
xstp=1000;
ystp=1000;?%沒步的大小??
for?i=0:xstp:2000
??ymax=1000;?%?ymax=-0.12*i+150;
??ymin=0;????%?ymin=0.12*i;
????
????for?j=0:ystp:1000
????????if?j>=ymin?&&?j<=ymax
????????????k=k+1;
????????????coord（k:）=[ij];
????????????plot（ij?‘.‘）;
????????????text（i?j?num2str（k）?‘fontsize‘?8）;
????????????hold?on;
????????end
????????if?abs（j-ymin）????????????k=k+1;
????????????coord（k:）=[iymin];
????????????plot（iymin?‘.‘）;
????????????text（iyminnum2str（k）?‘fontsize‘?8）;
?????????????hold?on;
????????end
?????????if?abs（j-ymax）ymax
????????????k=k+1;
????????????coord（k:）=[iymax];
????????????plot（iymax?‘.‘）;
????????????text（iymaxnum2str（k）?‘fontsize‘?8）;
?????????????hold?on;
????????end
????end
end
length（coord）
plot（[02000][00]）;
plot（[02000][10001000]）;
axis?equal;??????????
pp=0;
%{
????*

??????1st?column?—?The?unique?values?of?x
????*

??????2nd?column?—?The?number?of?instances?of?each?value
????*

??????3rd?column?—?The?percentage?of?each?value
%}
kk=tabulate（coord（:1））

lnods=setelment（kkystp）
disp（‘element?number?is‘）
bb=size（lnods）;
ofile1=fopen（‘Meshout.txt‘‘w‘）;
fprintf（ofile1‘Results\n‘）;
for?i=1:bb（1）?
fprintf（ofile1‘%d?%d?%d?%d?%d\n‘ilnods（i:））;
end
fprintf（ofile1‘coords\n‘）;
dd=size（coord）;
for?i=1:dd（1）
????fprintf（ofile1‘%d?%9.2f?%9.2f\n‘icoord（i:））;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1441??2020-09-20?08:36??meshcoord.m
?????文件?????????167??2020-09-20?08:36??Meshout.txt