exec('/mnt/media/work/vsu/Информационная защита/code/ellipse.sci');

p = 751;
a = -1;
P = [56, 332];
Q = [69, 241];
R = [83, 373];
out = sumpq(sumpq(multkp(2,P,p,a),multkp(3,Q,p,a),p,a),minusp(R,p),p,a);
mprintf('2P + 3Q – R = (%i, %i)\n', out(1), out(2));

P = [43,527];
k = 107;
out = multkp(k,P,p,a);
mprintf('kP = (%i, %i)\n', out(1), out(2));