function n=refract(x,y)
%% indice de refractionw
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n=ones(size(y)) ;