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

mertebe := proc( x1, y1p , , k)
y2 := int;
local xyx2y2nm
, , , , , ; n := 2;
x := int;
if y1 ≠ 0 then
y := int;
n := n + 1;
x2 := int;
m := ( 12 / × ( 3× x1^2− k^2)/ y1) mod p;
x := ( m^2 − 2× x1) mod p;
y := ( m × ( x1 − x) − y1) mod p;
x2 := x;
y2 := y;
while x≠ x1 or y ≠ p− y1do
n := n + 1;
y2 := y
m := ( y− y1 )/( x− x1) mod p; end do ;
x := ( m^2 − x− x1) mod p;
print ( "mertebe:" , n)
y := ( m × ( x2 − x) − y2) mod p; else print ( "mertebe:" , n)
x2 := x;
end proc
end if
> hesapla:=proc(p,k);
xgec:=int;ygec:=int;l:=int;
> for xgec from 0 to p-1 do
> with(numtheory);
> l:=xgec^3-(k^2)*xgec;
> ygec:=msqrt(l,p);
> if ygecFAIL then
if ygec0 then
print([xgec,ygec],[xgec,-ygec]);
else
print([xgec,ygec]);
end if;
mertebe(xgec,ygec,p,k);
> end if;
> end do;
> end proc;
73