funcaotransicao( x, y, frequencia)

function [saida, s_saida] = funcaotransicao( x, y, frequencia)
%
% Gerao da funo de transio complexa F(S) 
% que torna o coeficiente de difraco 
% Dt(teta)=D(teta).F(S) contnua e tata = teta(x,y).
%
% function [saida, s] = funcaotransicao( x, y, frequencia)
%
% saida       : Vector de sada com o funo de transio F(S);
% s_saida     : Vector de sada com os valores do parametro 'S';
% x	        : Distncia horizontal da descontinuidade espacial;
% y	        : Vector com as elevaes em relao ao plano horizontal
%		          que contm a descontinuidade espacial (edge);
% frequencia  : Valor da frequencia de trabalho (MHz).
%

saida = y*0;
s_saida = y*0;
s = 0;
tgteta2 = 0;
ro = 0;
k = 0;
fs = 0;
lambda = 300/frequencia;
beta = 2*pi/lambda;
csi = 0;
k1 = 0;
k2 = 0;
k3 = -1/sqrt(2*pi*beta);
f = 0;
g = 0;
num1 = 0;
den1 = 0;
num2 = 0;
den2 = 0;
cdaux = 0;
num3 = 0;
den3 = 0;
costeta = 0;
senteta = 0;
for k = 1 : length(y),
  ro = power( (x*x + y(k)*y(k)), 0.5);
  tgteta2 = power( (y(k)/x), 2);   
  s = beta*ro*tgteta2/2;
  k1 = power( (2*pi*s), 0.5);
  csi = power( (2*s/pi), 0.5);
  num1 = 1 + 0.926*csi;
  den1 = 2 + 1.792*csi + 3.104*power( csi, 2);
  f = num1/den1;
  num2 = 1;
  den2 = 2 + 4.142*csi + 3.492*power( csi, 2) + 6.670*power( csi, 3);
  g = num2/den2;
  k2 = f + j*g;
  fs = k1*k2;
  saida(k) = fs;
  s_saida(k) = s;
end