? e=ellinit(Mod(1,101)*[0,0,0,0,7])
%1 = [Mod(0, 101), Mod(0, 101), Mod(0, 101), Mod(0, 101), Mod(7, 101), Mod(0, 101), Mod(0, 101), Mod(28, 101), Mod(0, 101), Mod(0, 101), Mod(12, 101), Mod(42, 101), Mod(0, 101), Vecsmall([3]), [101, [0, 59, [6, 0, 0, 0]]], [0, 0, 0, 0]]
? issquare(Mod(1,101)^3+Mod(7,101),&c)
%2 = 0
? issquare(Mod(2,101)^3+Mod(7,101),&c)
%3 = 0
? issquare(Mod(3,101)^3+Mod(7,101),&c)
%4 = 0
? issquare(Mod(4,101)^3+Mod(7,101),&c)
%5 = 1
? c
%6 = Mod(24, 101)
? G=[Mod(4,101),c]
%7 = [Mod(4, 101), Mod(24, 101)]
? baby=vector(12,j,ellpow(e,G,j))
%8 = [[Mod(4, 101), Mod(24, 101)], [Mod(94, 101), Mod(88, 101)], [Mod(26, 101), Mod(3, 101)], [Mod(50, 101), Mod(75, 101)], [Mod(52, 101), Mod(15, 101)], [Mod(62, 101), Mod(50, 101)], [Mod(34, 101), Mod(74, 101)], [Mod(77, 101), Mod(90, 101)], [Mod(91, 101), Mod(44, 101)], [Mod(74, 101), Mod(76, 101)], [Mod(75, 101), Mod(56, 101)], [Mod(99, 101), Mod(10, 101)]]
? lift(%)
%9 = [[4, 24], [94, 88], [26, 3], [50, 75], [52, 15], [62, 50], [34, 74], [77, 90], [91, 44], [74, 76], [75, 56], [99, 10]]
? giant=vector(12,j,ellpow(e,G,-12*j))
%10 = [[Mod(99, 101), Mod(91, 101)], [Mod(69, 101), Mod(93, 101)], [Mod(56, 101), Mod(61, 101)], [Mod(63, 101), Mod(68, 101)], [Mod(98, 101), Mod(92, 101)], [Mod(6, 101), Mod(18, 101)], [Mod(84, 101), Mod(12, 101)], [Mod(62, 101), Mod(50, 101)], [Mod(62, 101), Mod(51, 101)], [Mod(84, 101), Mod(89, 101)], [Mod(6, 101), Mod(83, 101)], [Mod(98, 101), Mod(9, 101)]]
? lift(%)
%11 = [[99, 91], [69, 93], [56, 61], [63, 68], [98, 92], [6, 18], [84, 12], [62, 50], [62, 51], [84, 89], [6, 83], [98, 9]]
? for(i=1,12,if(baby[i]=Mod(1,101)*[62,50],print(i))
  ***   syntax error, unexpected $end, expecting )-> or ',' or ')': 
  ***   ...d(1,101)*[62,50],print(i))
  ***                               ^-
? for(i=1,12,if(baby[i]=Mod(1,101)*[62,50],print(i)))
1
2
3
4
5
6
7
8
9
10
11
12
? for(i=1,12,if(baby[i]==Mod(1,101)*[62,50],print(i)))
1
2
3
4
5
6
7
8
9
10
11
12
? baby=vector(12,j,ellpow(e,G,j))
%14 = [[Mod(4, 101), Mod(24, 101)], [Mod(94, 101), Mod(88, 101)], [Mod(26, 101), Mod(3, 101)], [Mod(50, 101), Mod(75, 101)], [Mod(52, 101), Mod(15, 101)], [Mod(62, 101), Mod(50, 101)], [Mod(34, 101), Mod(74, 101)], [Mod(77, 101), Mod(90, 101)], [Mod(91, 101), Mod(44, 101)], [Mod(74, 101), Mod(76, 101)], [Mod(75, 101), Mod(56, 101)], [Mod(99, 101), Mod(10, 101)]]
? for(i=1,12,if(baby[i]==Mod(1,101)*[62,50],print(i)))
6
? for(i=1,12,if(giant[i]==Mod(1,101)*[62,50],print(i)))
8
? ellpow(e,G,102)
%17 = [0]
? ellcard(G)
  ***   at top-level: ellcard(G)
  ***                 ^----------
  *** ellcard: incorrect type in checkell (t_VEC).
  ***   Break loop: type 'break' to go back to GP prompt
break> break

? ellcard(e)
%18 = 102
? for(j=1,101,if(ellpow(e,G,j)==[0],print(j)))
? issquare(Mod(5,101)^3+Mod(7,101),&c)
%20 = 1
? c
%21 = Mod(43, 101)
? Q=[Mod(5,101),c]
%22 = [Mod(5, 101), Mod(43, 101)]
? giant=vector(12,j,elladd(e,Q,ellpow(e,G,-12*j)))
%23 = [[Mod(77, 101), Mod(90, 101)], [Mod(50, 101), Mod(26, 101)], [Mod(23, 101), Mod(16, 101)], [Mod(89, 101), Mod(81, 101)], [Mod(15, 101), Mod(94, 101)], [Mod(8, 101), Mod(32, 101)], [Mod(48, 101), Mod(2, 101)], [Mod(25, 101), Mod(52, 101)], [Mod(14, 101), Mod(78, 101)], [Mod(94, 101), Mod(88, 101)], [Mod(74, 101), Mod(25, 101)], [Mod(35, 101), Mod(82, 101)]]
? for(i=1,12,for(j=1,12,if(baby[i]==giant[j],print(i," ",j))))
2 10
8 1
? ellpow(e,G,20)
%25 = [Mod(5, 101), Mod(43, 101)]
? Q
%26 = [Mod(5, 101), Mod(43, 101)]
? \q