tc balıkesir üniversitesi fen bilimleri enstitüsü matematik anabilim ...

hesapla := proc( p, k)
for xgec from 0 to p − 1 do
local xgec, ygec, l;
with( numtheory ) ;
xgec := int;
l := xgec^3 − k^2× xgec;
ygec := int;
ygec := msqrt ( l, p ) ;
l := int;
if ygec ≠ FAIL then
if ygec ≠ 0 then print ( [ xgec, ygec ] , [ xgec, −ygec]
)
else print ( [ xgec, ygec]
)
end if ;
mertebe ( xgec, ygec, p, k)
end if
hesapla(5,1);
[ 00 , ]
"mertebe:" , 2
[ 10 , ]
"mertebe:" , 2
[ 21 , ] , [ 2-1 , ]
"mertebe:" , 4
[ 32 , ] , [ 3-2 , ]
"mertebe:" , 4
[ 40 , ]
"mertebe:" , 2
74
end do
end proc