%****讀入數據****%
Pt=GetData;??%調用GetData函數，讀入數據
ptn=num2str（Pt（:1））;?%取點號
ptx=Pt（:2）;?%取x
pty=Pt（:3）;?%取y
ptz=Pt（:4）;?%取z
step=input（‘在此輸入步長：\n‘）;?%得到網格間距
x_max=max（ptx）;?y_max=max（pty）;?%計算xy最大值
x_min=min（ptx）;?y_min=min（pty）;?%計算xy最小值
x0=floor（x_min/step）*step;???y0=floor（y_min/step）*step;?%Grid起始點
nx=ceil（（x_max-x0）/step）;???ny=ceil（（y_max-y0）/step）;???%網格數量
l_x=nx*step;??l_y=ny*ste