function [x] = L2_ball(rho) dt = 2*pi/500; t = 0:dt:2*pi; x(1,:) = rho*cos(t); x(2,:) = rho*sin(t);