clc
clear?all
vp1=1600;
vp2=2800;
vs1=800;
vs2=1600;
p1=1;
p2=1.5;
c11=-0.1;
c12=-0.1;

?for?i=1:80
rs（1i）=i;
x1=i;
y1=asind（sin（x1*pi/180）*vs1/vp1）;
y2=asind（sin（x1*pi/180）*vs2/vp1）;
x2=asind（sin（x1*pi/180）*vp2/vp1）;
A=[（sin（x1*pi/180））（cos（y1*pi/180））（-sin（x2*pi/180））（cos（y2*pi/180））];
B=[cos（x1*pi/180）-sin（y1*pi/180）cos（x2*pi/180）sin（y2*pi/180）];
C=[sin（2*x1*pi/180）（vp1/vs1）*cos（2*y1*pi/180）（p2*（vs2）^2*vp1）/（p1