function MM = MM(a, b, c)
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
am(1) = 0.325303;
am(2) = 0.4211071;
am(3) = 0.1334425;
am(4) = 0.006374323;
bm(01) = 0.1337764;
bm(2) = 0.6243247;
bm(3) = 1.3425378;
bm(4) = 2.2626645;
MM = 0;
aa = a/sqrt(2*(1 - c*c));
bb = b/sqrt(2*(1 - c*c));
for i=1:4
    for j=1:4
        MM = MM + am(i)*am(j)*f(bm(i), bm(j), aa, bb, c)
    end
end
MM = MM*sqrt(1 - c*c)/3.121592654;
end

