sdp p1

clc
clear

x = [1 2 3 4; 0.48 0.27 0.13 0.12];
if(sum(x(2, :)) == 1)
    med = sum(x(1, :) .* x(2, :)) / 4;
    med
    disp = (x(1, :) - med).^2 .* x(2, :);
    disp
end