задача подпрограмма-процедура в Delphi

[Ксюшка 76]

Заданы экспериментальные значения четырёх случайных величин a={a1,a2,a3}, b={b1,b2,b3,b4}, c={c1,c2,c3,c4,c5}, p={p1,p2,p3,p4}
Вычислить d(x)=1/n*(сумма от i=1 до n)*xi^2-M^2 , где M(x)-среднее значение каждой случайной величины.Вычисление среднего значения случайной величины оформить в виде подпрограммы- процедуры.
Исходные данные:a={1,2,3},b={-5;6;1.2;2.4},c={-3.2;0;1;-1;2},p={2;3;1;4}.