Значит вы не так считаете. Вот листинг для
Maple 14Код:
restart;
p:=proc(n)
local L,a,b,c:
L:=[]:
for a from 1 to n do
for b from 1 to n do
c:=sqrt(a^2+b^2):
if is(c,integer) then L:=[op(L),[1,binomial(a-1,1)+binomial(b-1,1)-binomial(c-1,1),binomial(a-1,2)+binomial(b-1,2)-binomial(c-1,2)]]:fi:
od:od:
plots[pointplot3d](L,axes=normal,color=red,symbol=point,orientation=[-40,50],labels=[x,y,z]);
end proc:
p(100);
Берем тройки

для которых выполняется

и строим для каждой точку с координатами

где



а

- биномиальные коэффициенты.
Проверьте.
-- Ср сен 15, 2010 12:03:51 --Кстати, попробуйте сделать то же самое для

, это любопытно. Там, действительно, не точка. Т.е. точка если положить

. А вообще, получается вертикальная прямая, но точки на ней расположены не равномерно - расстояние между ними увеличивается. Тогда как для

точки на (уже не вертикальной) прямой расположены равномерно. Особенно интересны графики "два в одном" - и тройки и их отображения.
Вот листинг для

Код:
restart;
p:=proc(n)
local L,a,b,c:
L:=[]:
for a from 1 to n do
for b from 1 to n do
c:=a+b:
if is(c,integer) then L:=[op(L),([a,b,c],[1,binomial(a-1,1)+binomial(b-1,1)-binomial(c-1,1),binomial(a-1,2)+binomial(b-1,2)-binomial(c-1,2)])]:fi:
od:od:
plots[pointplot3d](L,color=blue,symbol=point,orientation=[-40,50],labels=[x,y,z]);
end proc:
p(10);
а вот для

Код:
restart;
p:=proc(n)
local L,a,b,c:
L:=[]:
for a from 1 to n do
for b from 1 to n do
c:=sqrt(a^2+b^2):
if is(c,integer) then L:=[op(L),([a,b,c],[1,binomial(a-1,1)+binomial(b-1,1)-binomial(c-1,1),binomial(a-1,2)+binomial(b-1,2)-binomial(c-1,2)])]:fi:
od:od:
plots[pointplot3d](L,color=blue,symbol=point,labels=[x,y,z]);
end proc:
p(100);
Посмотрите, это красиво

Обратите внимание на рисунок пифагоровых троек на втором графике.