function nouburjaold(N)
Mm=2;
for x1=1:N
    for y1=1:N
        x(1,:)=[x1,y1];
        for x2=1:N
            for y2=1:N
                x(2,:)=[x2,y2];
                for x3=1:N
                    for y3=1:N
                        x(3,:)=[x3,y3];
                        for x4=1:N
                            for y4=1:N
                                x(4,:)=[x4,y4];
                                for x5=1:N
                                    for y5=1:N
                                        x(5,:)=[x5,y5];
                                        for x6=1:N
                                            for y6=1:N
                                                x(6,:)=[x6,y6];
                                                m=100; M=0;
                                                for i=1:5
                                                    for j=i+1:6
                                                        d=norm(x(i,:)-x(j,:));
                                                        if d<m
                                                            m=d;
                                                        end
                                                        if d>M
                                                            M=d;
                                                        end
                                                    end
                                                end
                                                if m>=N/10
                                                    if M/m<Mm
                                                        Mm=M/m;
                                                        y=x;
                                                    end
                                                end
                                            end
                                        end
                                    end
                                end
                            end
                        end
                    end
                end
            end
        end
    end
end
fprintf('\nEl mínim valor de M/m és %f\n',Mm);
plot(y(:,1),y(:,2),'r');