Matlab - Algoritmo multiplicacion de matrices Matlab

clc

a=[-0.1473   -0.0775   -0.0233;-0.0388   -0.1783    0.0465;-0.0620   -0.0853   -0.3256]

b=[1 2 3;4 5 6;7 8 9]

[m r]=size(a)

[s n]=size(b)

ab=zeros(m,n)

for j=1:m

for k=1:n

    abo=0

    for i=1:r

    abo=abo+a(j,i)*b(i,k)

    end

    ab(j,k)=abo

end

end

    ab(j,k)=abo

end

end

for i=0:255

    [x , y] = find(C4(:,:,1)==i);Pr(i+1)= length(x)/380/640;

    [x , y] = find(C4(:,:,2)==i);Pg(i+1)= length(x)/380/640;

    [x , y] = find(C4(:,:,3)==i);Pb(i+1)= length(x)/380/640;

end

for k = 1:255

    for i=1:k

        w1r = sum(Pr(1:i)); w1g=sum(Pg(1:i)); w1b=sum(Pb(1:i));

        u1r = sum([1:i].*Pr(1:i)/w1r);

        u1g = sum([1:i].*Pg(1:i)/w1g);

        u1b = sum([1:i].*Pb(1:i)/w1b);

    end

    for i=k+1:255

        w2r = sum(Pr(i+1:255));w2g = sum(Pg(i+1:255)); w2b = sum(Pb(i+1:255));

        u2r = sum([i+1:255].*Pr(i+1:255)/w2r);

        u2g = sum([i+1:255].*Pg(i+1:255)/w2g);

        u2b = sum([i+1:255].*Pb(i+1:255)/w2b);

    end

    uTr = w1r*u1r + w2r*u2r;

    uTg = w1g*u1g + w2g*u2g;

    uTb = w1b*u1b + w2b*u2b;

    maxVARr(k) = w1r*(u1r-uTr)^2 + w2r*(u2r-uTr)^2;

    maxVARg(k) = w1g*(u1g-uTg)^2 + w2g*(u2g-uTg)^2;

    maxVARb(k) = w1b*(u1b-uTb)^2 + w2b*(u2b-uTb)^2;

end