clc; lamda=5893*1e-10;//in meter ue=1.55333;//given ue u0=1.5442;//given u0 t=lamda/(2*(ue-u0));//Thicknesss of half wave plate disp(+'meter',t,'Thicknesss of half wave plate =');