Matlab - Ejercicio ternas

c=nchoosek(1:20,3);

c=c(sum(c(:,1:2).^2,2)-c(:,3).^2==0,:);

disp(c)

clear all

V=[];

for a=1:20

    for b=1:20

        for c=1:20

            if a^2+b^2==c^2

                V=[V; [a b c]];

               m= find(V(:,3)==c);

               if length(m)==2

                   V(end,:)=[];

               end

            end

        end

    end

end

V

>> ternas

V =

     3     4     5

     5    12    13

     6     8    10

     8    15    17

     9    12    15

    12    16    20