matlabpool?open;
tic;
clear;
i?=?sqrt?（-1）;
theta=0/180*pi;
%lamda=2142;
lamda=418;
%lamda=699.2393853;
%lamda=2097.718156;
%lamda=2097.718156;
%lamda=1048.859078;
%lamda=1048859078;
a=800;
h1=0.62*a;
h2=0.38*a;
%###########################
n0=1;
%###########################
n1=sqrt（1.7）;
theta1=asin（n0*sin（theta）/n1）;
%###########################
n2=sqrt（1.75）;
?
theta2=asin（n1*sin（theta1）/n2）;
%###########################
?
%##########################################################################
ng=1;
theta_g=asin（n2*sin（theta2）/ng）;
%###########################
?
?
del1=2*pi/lamda*n1*h1*cos（theta1）;
del2=2*pi/lamda*n2*h2*cos（theta2）;
phi0=n0*cos（theta）;
phi1=n1*cos（th