function y = RBF1out( x, a, mu, sigma )

%function y = RBF1out( x, a, mu, sigma )
%  output by RBF 
%  with parameter 'a_j',  'mu', 'sigma'

[d,M] = size( a ); % d; target vector dimension
                   % M; number of base centers
[q,T] = size( x ); % q; data vector dimension

y = zeros( d, T );
for t = 1:T
  for j = 1:M
    x_mu = x(:,t) - mu(:,j);
    z = exp( -x_mu'*x_mu/(2*sigma^2) );
    y(:,t) = y(:,t) + a(:,j)*z;
  end
end