function  out=Michalewicz(X)
 


%I=(X>0);
%X=I.*X;
n=size(X,2)-1;
r=[]; t=2/pi;
for i=1:n
r=[   sin(X(:,i+1)).*sin(t*X(:,i+1).^2).^20 +  sin(X(:,i)).*sin( (X(:,i+1).^2)/pi ).^20  ,r];
end

out=sum(r,2);


    