function x_new=int_point_quiz(n,eps)
c=ones(n,1);
c(n)=-1;
x_old=max(eps,zeros(n,1));
lambda_old=max(eps,zeros(n,1));
mu_old=0;
for iter=1:10
  mat1 = horzcat(sparse(n,n),-speye(n,n),-???);
  mat2 = horzcat(ones(1,n),sparse(1,n),sparse(1,1));
  mat3 = horzcat(sparse(1:n,1:n,lambda_old,n,n),???,sparse(n,1));
  mat  = vertcat(mat1,mat2,mat3);
  rhs1 = ???;
  rhs2 = ones(1,n)*x_old-1;
  rhs3 = sparse(n,1);
  for i=1:n
   rhs3(i)=???;
  end
  rhs = vertcat(rhs1,rhs2,rhs3);
  sol=mat\rhs;
  x_new=max(eps,x_old-sol(1:n));
  lambda_new=max(eps,lambda_old-sol(n+1:2*n));
  mu_new=mu_old-sol(2*n+1:2*n+1);
  discrep=norm(x_new-x_old)/norm(x_new);
  fprintf ("iter: %d Discrepancy: %f \n",iter,discrep);
  %fprintf ("iter: %d Discrepancy: %f \n",iter,discrep); %matlab
  x_old=x_new;
  lambda_old=lambda_new;
  mu_old=mu_new;
  if (discrep<0.001)
   break
  end
end
end


