clear; 
Fq:=GF(131);
a:=0; b:=5;
A2<x,y>:=AffineSpace(Fq,2);
Eaff:=Curve(A2,[y^2-(x^3+a*x+b)]);
P:=Eaff![2,12];
Q:=Eaff![67,56];
lambda:=(Q[2]-P[2])/(Q[1]-P[1]);
nu:=Q[2]-lambda*Q[1];
l:=Curve(A2,[y-(lambda*x+nu)]);
IntersectionPoints(Eaff,l);

P2<X,Y,Z>:=ProjectiveSpace(Fq,2);
E:=EllipticCurve([Fq|a,b]);
F<x,y>:=FunctionField(E);
l:=y-(lambda*x+nu);
D:=Divisor(l);
Support(D);
//or
Zeros(l);
Poles(l);
