Составить подпрограммы для выполнения операций над векторами

Program vector;
uses crt;
type
Vec=record
 x:integer;
 y:integer;
 end;
 
var
a,b:Vec;
 
procedure SUM(const a1,b1:Vec); {сложение векторов}
var c:Vec;
begin
c.x:=a1.x+b1.x;
c.y:=a1.y+b1.y;
writeln('Результат суммы векторов А и В: вектор С(',c.x,',',c.y,').');
end;
 
procedure SUB(const a1,b1:Vec); {вычитание векторов}
var c:Vec;
begin
c.x:=a1.x-b1.x;
c.y:=a1.y-b1.y;
writeln('Результат разности векторов А и В: вектор С(',c.x,',',c.y,').');
end;
 
procedure MUL(const a1,b1:Vec); {скалярное произведение векторов}
begin
writeln('Результат скалярного произведения векторов А и В: ', a1.x*b1.x+a1.y+b1.y, '.');
end;
 
procedure MUL2(const a1,b1:Vec); {умножение векторов на число}
var m:integer;
begin
write('Введите число: ');readln(m);
writeln('Результат умножения вектора А на',m,' : вектор С (',a1.x*m,',',a1.y*m,').');
writeln('Результат умножения вектора B на',m,' : вектор С (',b1.x*m,',',b1.y*m,').');
end;
 
procedure LEN(const a1,b1:Vec); {длина векторов}
begin
writeln('Длина вектора А: ',sqrt(a1.x*a1.x+a1.y*a1.y):5:4,'.');
writeln('Длина вектора B: ',sqrt(b1.x*b1.x+b1.y*b1.y):5:4,'.');
end;
 
BEGIN
clrscr;
writeln('ВЕКТОР А');
write('X: ');readln(a.x); {x для А}
write('Y: ');readln(a.y); {y для А}
writeln;
writeln('ВЕКТОР В'); 
write('X: ');readln(b.x); {x для B}
write('Y: ');readln(b.y); {y для В}
writeln;
 
SUM(a,b); {вызов соответствующих процедур}
SUB(a,b);
MUL(a,b);
MUL2(a,b);
LEN(a,b);
 
readkey;
END.