Here we see how the lines of a projective plane PG(2, q2 ) meet a hermitian curve. It is well known that every line meets in either 1 or q+1 points.
gap> h := HermitianVariety(2,49);H(2, 7^2)
gap> pg := AmbientSpace(h);
PG(2, 49)
gap> lines := Lines(pg);
<lines of PG(2, 49)>
gap> curve := AsList(Points(h));;
#I Computing nice monomorphism...
gap> Size(curve);
344
gap> Collected(List(lines,t -> Number(curve,c -> c in t)));
[ [ 1, 344 ], [ 8, 2107 ] ]
