clear?all
clc
tic;
t=toc;
format?long?g;
A=load（‘C:\Users\WYG704\Desktop\stage35\r1.dat‘）;??????%讀取數據文件
x=A（:1）;??????????????????????????????????????
y=A（:2）;??????????????????????????????????????????????????????????????????????????????????
z=A（:3）;
A=[x?y?z];
n?=size（x）;???%求控制點總數
%判斷開曲線或是閉曲線并作相應處理
d=sqrt（（x（1）-x（n））^2+（y（1）-y（n））^2+（z（1）-z（n））^2）;
if??d>0.01
??%開曲線端點增加處理
????????x（n+2）=x（2）;y（n+2）=y（2）;z（n+2）=z（2）;
??????for?i=n+1:-1:2
????????x（i）=x（i-1）;y（i）=y（i-1）;z（i）=z（i-1）;
??????end
??else
????%閉曲線端點增加處理
???????x（n+2）=x（2）;y（n+2）=y（2）;z（n+2