function Pytha( nMax )
   fprintf( 'a^2 + b^2 = c^2: ganzzahlig \n' );
   fprintf( '   a   |   b   |   c  \n' );
   fprintf( '----------------------\n' );
   for( a = 1:nMax )
     bMax = sqrt( nMax^2 - a^2 );  
     for( b = a:bMax )         
        c = sqrt( a^2 + b^2 );
        % c ist ganzzahlig 
        if( (c - floor(c)) < eps )
           fprintf( ' %3d  |  %3d  |  %3d \n', a, b, int32(c) );  
        end
     end
   end
end