clear all;
close all;
clc;
h0=0.01;
beta=0.002589*10^(-8);
loss=0.01601;
density=1.428;
v=1.1020;
delta2=0.002;
gam=3.25;
c=1500;
f=40*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.0061636*10^(-8);
loss=0.01856;
density=1.966;
v=1.133;
delta2=0.004;
gam=2.92;
c=1500;
f=40*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;