function plota_phi(tipok,toroidal,tensorial,dobrado,C,R,CD,TD,LD)
nx = TD(1)*4;
ny = TD(2)*4;
x = LD(1)*(0:nx)/nx + (CD(1)-LD(1)/2)*ones(1,nx+1);
y = LD(2)*(0:ny)/ny + (CD(2)-LD(2)/2)*ones(1,ny+1);

[x,y] = meshgrid(x,y);
for i = 1:size(x,1)
    for j = 1:size(x,2)
        z = [x(i,j), y(i,j)];
            f(i,j) = phi(z,tipok,toroidal,tensorial,dobrado,...
                C,R,TD,LD);
    end
end

surf(x,y,f)
set(gca,'Fontsize', 20)
