%
% Reduction of BDFIM model from three loops to one loop
%
Lrl=[0.72 0.576 0.576;0.576 1.878 1.727;0.576 1.727 3.037]*1e-5;
Rrl=[1.056 0 0;0 1.209 0;0 0 1.361]*1e-4;
Mpl=[0.5793 1.6693 2.5533]*1e-3;  
Mcl=[0.5555 1.4137 1.6072]*1e-3;
[v,d]=eig(Lrl);
x=[0;0;1];
Lr=x'*d*x
Mp=Mpl*v*x
Mc=Mcl*v*x
Rr=x'*(v'*Rrl*v)*x