colorImage = imread('1.jpg');


R = colorImage(:, :, 1);
G = colorImage(:, :, 2);
B = colorImage(:, :, 3);

R = double(R);
G = double(G);
B = double(B);
L = R + G + B;
[H,W] = size(L)
r = R./L;
g = G./L;

r = int16(r*255)+1;
g = int16(g*255)+1;


h = zeros(256,256);
for j=1:H
    for i=1:W
        wr = r(j,i);
        wg = g(j,i);
        h(wg,wr) = h(wg,wr)+1;
    end
end
mesh(h)
S=0;mi_r=0;mi_g=0;delta_rr=0;delta_gg=0;delta_rg=0;
for j=1:256
    for i=1:256
        S=S+h(j,i);
        mi_r=mi_r+h(j,i)*j;
    end
end
mi_r=mi_r/S;
for j=1:256
    for i=1:256
        mi_g=mi_g+h(j,i)*i;
    end
end
mi_g=mi_g/S;
for j=1:256
    for i=1:256
        delta_rr=delta_rr+h(j,i)*(j-mi_r)^2;
    end
end
delta_rr=delta_rr/S;
for j=1:256
    for i=1:256
        delta_gg=delta_gg+h(j,i)*(i-mi_g)^2;
    end
end
delta_gg=delta_gg/S;
for j=1:256
    for i=1:256
        delta_rg=delta_rg+h(j,i)*(j-mi_r)^2*(i-mi_g)^2;
    end
end
delta_rg=delta_rg/S;