clear all;
close all;
clc;
h0=0.01;
beta=0.005999*10^(-8);
loss=0.01643;
density=2.176;
v=1.2323;
delta2=0.002;
gam=3.25;
c=1500;
f=20*10^3;
lamda=c/f;
k=2*pi/lamda;
theta=(10:1:90)*pi/180;
[delta_whole,sb_kr,sb_cr,sb_mr,sb_v,sb]=Sb(h0,beta,loss,density,v,delta2,gam,k,theta);
plot(theta*180/pi,sb);
hold on;

h0=0.01;
beta=0.005874*10^(-8);
loss=0.04534;
density=1.993;
v=1.113;
delta2=0.005;
gam=2.47;
c=1500;
f=20*10^3;
lamda=c/f;
k=2*pi/lamda;
theta=(10:1:90)*pi/180;
[delta_whole,sb_kr,sb_cr,sb_mr,sb_v,sb]=Sb(h0,beta,loss,density,v,delta2,gam,k,theta);
plot(theta*180/pi,sb,'--');
grid on;