clc ;
clear ;
b =2.5;
a =5;
c =3*10^10;
w0 =4.5;
disp ( ' For a TE10 mode which is the dominant mode ' );
wc =2* a ;
wg =( w0 / sqrt (1 -( w0 / wc ) ^2) ) ;
Vp =( c/ sqrt (1 -( w0 / wc ) ^2) ) /10^10;
B =((2*( %pi ) * sqrt ( wc ^2 - w0 ^2) ) /( w0 * wc ) ) ;
mprintf ('\n Guide wavelength is wg=%1.5f cm \nPhase constant is B=%1.3f radians \nPhase Velocity is Vp=%1.2f ∗10ˆ10 cm/sec',wg,B ,Vp);