Sbar    = [3.017e-03,2.374e-01,0.3,log(.17)]';
P       = zeros(4);
Q       = zeros(4);
P(2,2)  = 9.213e-01;
P(4,4)  = 9.726e-01;
Q(1,1)  = 9.264e-03;
Q(2,2)  = 9.447e-03;
Q(4,4)  = 1.994e-02;
P0      = (eye(4)-P)*Sbar;

Sbar    = [5.155e-03,2.432e-01,0.3,log(.17)]';
P       = zeros(4);
Q       = zeros(4);
P(2,2)  = 9.523e-01;
P(4,4)  = 9.686e-01;
Q(1,1)  = 1.306e-02;
Q(2,2)  = 1.361e-02;
Q(4,4)  = 2.057e-02;
P0      = (eye(4)-P)*Sbar;


