# I want to find a generator set of the $l$-torsion points of an elliptic curve defined over a finite field $F_q$

I have the elliptic curve $E: y^2=x^3+x$ defined over a field $F_q$ with $q=p^2$ for a certain prime $p$, and want to calculate the $l$-torsion points, or in other words $E[l]$, for an integer $l$ (probably a prime different from $p$.)

Welcome to Ask Sage. Thank you for your question.

Ideally, provide a minimal concrete example, with code to define the objects, that can be copied and pasted in a fresh Sage session, leaving only the final step to answer.

`p=18446744073709551667`

`E = EllipticCurve(GF(p^2),[1,0])`

E is an Elliptic Curve defined by y^2 = x^3 + x over Finite Field in z2 of size 18446744073709551667^2 and I want to find generators $P$ and $Q$ for $E[4999]$, but I can't find a suitable code.use

`sage: p=E.torsion_polynomial(4999)`

maybe ?