Как построить эллиптические кривые над конечным полем с помощью Matlab

Мне нужно нарисовать эллиптическую кривую над конечным полем F17 (другими словами, я хочу нарисовать некоторые конкретные точки на кривой), но почему-то я этого не делаю. я не понял.

Кривая определяется уравнением:

y^2 = x^3 +x + 1 (mod 17)

Я пробовал описанный ниже способ, но он не сработал.

для x = 0:16, plot (x, mod (sqrt (x ^ 3 + x + 1), 16), 'r') ', end

Кто-нибудь может помочь?

[Обновление]

Согласно предложениям Натана и Билла, вот немного измененная версия.

 x = 0:18
 plot(mod(x,16), mod(sqrt(x.^3+x+1), 16),'ro')

Однако мне кажется, что цифра НЕПРАВИЛЬНАЯ , например, y не является целым числом, когда x = 4.