求椭圆曲线上的整点

function ecc(a,b,p)
%求椭圆曲线的整点
%例：考虑椭圆曲线上E11(1,6),即y^2=x^3+x+6和一个模数p=11定义的曲线，确定E11(1,6)中所有的点
%对应的命令为：ecc(1,6,11)
for x=0:p-1
    for y=0:p-1
        for k=-p^2:p^2
            if (x^3+a*x+b+k*p)==y^2
                [x,y]
            end        
        end
    end
end