%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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‘‘炮檢點‘‘A4:C4‘）;
X1=xlsread（‘v.xlsx‘‘炮檢點‘‘G3:I41‘）;%D3:F43

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

%n0=550;%格林函數及震源時間函數采樣點數;
nm=100;%循環次數！




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%??????畫圖????????%%%%%%%%%%%%%%%%%%%%%%%%%%%%
?[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