function [g,y]=cag(x,mu,nivel);
N=length(x);
h=zeros(N,1);
g=zeros(N,1);
y=zeros(N,1);
for k=2:N
   h(k)=h(k-1)+mu*(nivel-y(k-1).^2);
   if(h(k)<0)
      h(k)=0;
   end;
   g(k)=sqrt(h(k));
   y(k)=g(k)*x(k);
end;