%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc;clear?all;close?all;
v=xlsread（‘v.xlsx‘‘v-p‘‘A3:C20‘）;%分層速度，包括橫波與縱波
p=xlsread（‘v.xlsx‘‘v-p‘‘D3:D20‘）;%分層密度表示該深度至下一深度之間的密度
X0=xlsread（‘v.xlsx‘‘炮檢點(diǎn)‘‘A4:C4‘）;
X1=xlsread（‘v.xlsx‘‘炮檢點(diǎn)‘‘G3:I41‘）;%D3:F43

[n_sensorser]=size（X1）;
[n_sourcescr]=size（X0）;

%n0=550;%格林函數(shù)及震源時(shí)間函數(shù)采樣點(diǎn)數(shù);
nm=100;%循環(huán)次數(shù)！




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%??????畫(huà)圖????????%%%%%%%%%%%%%%%%%%%%%%%%%%%%
?[nnx0nny0]=size（v）;


xxx=[X0（:1）;X1（:1）];
yyy=[X0（:2）;X1（:2）];
zzz=[X0（:3）;X1（:3）;v（:1）];


x0max=max（xxx）;
y0max=max（yyy）;
z0max=max（zzz）;
x0m