How to use a formula for all elements of matrix?

I have a matrix the

Mij =[ NaN    2.5000         0
  2.5000       NaN   -0.8333
       0   -0.8333       NaN]

And I have a general formula::

phi=atand(abs((m2-m1)/(1+m1*m2)))

I need to use this formula for all elements of my matrix.

For example, I want to use as m2=M12(element in first row and second column)and

as m1=M32 (element in third row and second column). and according the formula

phi=atand(abs((Mij(1,2)-Mij(3,2))/(1+Mij(3,2)*Mij(1,2)));

And I need for all cases, which for all elements in matrix.

Can anyone help me?

I have tried as below, but the result is wrong:

for u=1:E1
    for v=1:E1
          p(v,u)=abs((Mij(v,u)-Mij(v,u))/(1+Mij(v,u)*Mij(v,u)));
    end
end

Result is wrong:

p =
 NaN     0     0
   0   NaN     0
   0     0   NaN

[Azzi Abdelmalek]

This is not clear,

Mij(v,u)-Mij(v,u)) is equal to 0

Answer 1

[Torsten]

for u1=1:size(M,1)
  for v1=1:size(M,2)
    m1 = M(u1,v1);
    for u2=1:size(M,1)
      for v2=1:size(M,2)
        m2 = M(u2,v2);
        phi(u1,v1,u2,v2)=atand(abs((m1-m2)/(1+m1*m2)));
      end 
    end 
  end
end

Best wishes

Torsten.

Answer 2

[Andrei Bobrov]

M=[NaN    2.5000         0
    2.5000       NaN   -0.8333
         0   -0.8333       NaN];

phi=@(m1,m2)atand(abs((m2-m1)./(1+m1.*m2)));

out1 = bsxfun(phi,M(:),M(:)');

out2 = reshape(out1,[size(M),size(M)]);

[Stephen23]

+1 very neat.