diff --git a/matlab/autoregressive_process_specification.m b/matlab/autoregressive_process_specification.m
index d802439e0..ec1cbc90d 100644
--- a/matlab/autoregressive_process_specification.m
+++ b/matlab/autoregressive_process_specification.m
@@ -75,17 +75,6 @@ function [InnovationVariance,AutoregressiveParameters] = autoregressive_process_
         AutoregressiveParameters(1) = t2-t3-Rho(1);
         AutoregressiveParameters(2) = (Rho(2)*Rho(2)-Rho(3)*Rho(1)-Rho(1)*Rho(1)+Rho(2))*t1 ;
         AutoregressiveParameters(3) = t3-Rho(1)+t2;
-      case 4
-        AutoregressiveParameters(1) = (Rho(1)*Rho(1)*Rho(1)*(Rho(4)+2)-Rho(1)*(Rho(3)*Rho(3)+1)+Rho(3)*Rho(4)-Rho(2)*Rho(2)*(Rho(3)-Rho(1)*Rho(4))-Rho(1)*Rho(1)*(Rho(3)+Rho(3)*Rho(4))+Rho(2)*(2*Rho(3)*Rho(1)*Rho(1)-3*Rho(1)*Rho(1)*Rho(1)+(Rho(3)*Rho(3)-2*Rho(4)+1)*Rho(1)+Rho(3))+Rho(2)*Rho(2)*Rho(2)*(Rho(1)-Rho(3)))/(2*Rho(1)*Rho(1)*Rho(1)*Rho(3)-Rho(1)*Rho(1)*Rho(1)*Rho(1)+2*Rho(1)*Rho(1)*Rho(2)*Rho(2)-4*Rho(1)*Rho(1)*Rho(2)-Rho(1)*Rho(1)*Rho(3)*Rho(3)+3*Rho(1)*Rho(1)+2*Rho(1)*Rho(2)*Rho(2)*Rho(3)-4*Rho(1)*Rho(2)*Rho(3)-Rho(2)*Rho(2)*Rho(2)*Rho(2)+2*Rho(2)*Rho(2)+Rho(3)*Rho(3)-1);
-        AutoregressiveParameters(2) = ((Rho(2)*Rho(2)-Rho(3)*Rho(1)-Rho(1)*Rho(1)+Rho(2))*(Rho(1)*Rho(1)-2*Rho(1)*Rho(3)+Rho(3)*Rho(3)+Rho(2)+Rho(4)-Rho(2)*Rho(4)-1))/(2*Rho(1)*Rho(1)*Rho(1)*Rho(3)-Rho(1)*Rho(1)*Rho(1)*Rho(1)+2*Rho(1)*Rho(1)*Rho(2)*Rho(2)-4*Rho(1)*Rho(1)*Rho(2)-Rho(1)*Rho(1)*Rho(3)*Rho(3)+3*Rho(1)*Rho(1)+2*Rho(1)*Rho(2)*Rho(2)*Rho(3)-4*Rho(1)*Rho(2)*Rho(3)-Rho(2)*Rho(2)*Rho(2)*Rho(2)+2*Rho(2)*Rho(2)+Rho(3)*Rho(3)-1);
-        AutoregressiveParameters(3) = (Rho(1)*Rho(1)*(2*Rho(2)*Rho(3)+Rho(3)*Rho(4))-Rho(3)+Rho(2)*Rho(2)*Rho(3)+Rho(2)*Rho(2)*Rho(2)*Rho(3)+Rho(1)*((Rho(4)-2)*Rho(2)*Rho(2)-Rho(2)*Rho(2)*Rho(2)+(2-Rho(4)-3*Rho(3)*Rho(3))*Rho(2)+Rho(4))+Rho(3)*Rho(3)*Rho(3)-Rho(1)*Rho(1)*Rho(1)*(Rho(4)-Rho(2)+1)-Rho(2)*Rho(3)*Rho(4))/(2*Rho(1)*Rho(1)*Rho(1)*Rho(3)-Rho(1)*Rho(1)*Rho(1)*Rho(1)+2*Rho(1)*Rho(1)*Rho(2)*Rho(2)-4*Rho(1)*Rho(1)*Rho(2)-Rho(1)*Rho(1)*Rho(3)*Rho(3)+3*Rho(1)*Rho(1)+2*Rho(1)*Rho(2)*Rho(2)*Rho(3)-4*Rho(1)*Rho(2)*Rho(3)-Rho(2)*Rho(2)*Rho(2)*Rho(2)+2*Rho(2)*Rho(2)+Rho(3)*Rho(3)-1);
-        AutoregressiveParameters(4) = (Rho(1)+Rho(3)/2-Rho(4)/2-(3*Rho(1)*Rho(2))/2+(Rho(1)*Rho(4))/2-(Rho(2)*Rho(3))/2+Rho(2)*Rho(2)-1/2)/(Rho(3)-Rho(1)*(2*Rho(2)+Rho(3)-1)+Rho(1)*Rho(1)+Rho(2)*Rho(2)-1)+(Rho(1)+Rho(3)/2+Rho(4)/2-(3*Rho(1)*Rho(2))/2+(Rho(1)*Rho(4))/2-(Rho(2)*Rho(3))/2-Rho(2)*Rho(2)+1/2)/(Rho(3)+Rho(1)*(Rho(3)-2*Rho(2)+1)-Rho(1)*Rho(1)-Rho(2)*Rho(2)+1)-1;
-      case 5
-           AutoregressiveParameters(1) = (Rho(1)*Rho(1)*Rho(1)*Rho(1)*Rho(1) - 3*Rho(1)*Rho(1)*Rho(1)*Rho(1)*Rho(3) + Rho(5)*Rho(1)*Rho(1)*Rho(1)*Rho(1) - 5*Rho(1)*Rho(1)*Rho(1)*Rho(2)*Rho(2) + 2*Rho(1)*Rho(1)*Rho(1)*Rho(2)*Rho(4) + 6*Rho(1)*Rho(1)*Rho(1)*Rho(2) + 2*Rho(1)*Rho(1)*Rho(1)*Rho(3)*Rho(3) - 2*Rho(5)*Rho(1)*Rho(1)*Rho(1)*Rho(3) - Rho(1)*Rho(1)*Rho(1)*Rho(4)*Rho(4) - Rho(1)*Rho(1)*Rho(1)*Rho(4) - 3*Rho(1)*Rho(1)*Rho(1) + Rho(1)*Rho(1)*Rho(2)*Rho(2)*Rho(3) + 2*Rho(5)*Rho(1)*Rho(1)*Rho(2)*Rho(2) + 2*Rho(1)*Rho(1)*Rho(2)*Rho(3)*Rho(4) + 3*Rho(1)*Rho(1)*Rho(2)*Rho(3) - 2*Rho(5)*Rho(1)*Rho(1)*Rho(2)*Rho(4) - 3*Rho(5)*Rho(1)*Rho(1)*Rho(2) + Rho(1)*Rho(1)*Rho(3)*Rho(3)*Rho(3) + Rho(5)*Rho(1)*Rho(1)*Rho(3)*Rho(3) + Rho(1)*Rho(1)*Rho(3)*Rho(4)*Rho(4) - Rho(1)*Rho(1)*Rho(3)*Rho(4) + Rho(1)*Rho(1)*Rho(3) + 2*Rho(5)*Rho(1)*Rho(1)*Rho(4) + 2*Rho(1)*Rho(2)*Rho(2)*Rho(2)*Rho(2) - Rho(1)*Rho(2)*Rho(2)*Rho(3)*Rho(3) + 2*Rho(5)*Rho(1)*Rho(2)*Rho(2)*Rho(3) + Rho(1)*Rho(2)*Rho(2)*Rho(4)*Rho(4) - 4*Rho(1)*Rho(2)*Rho(2)*Rho(4) - Rho(1)*Rho(2)*Rho(2) - 2*Rho(1)*Rho(2)*Rho(3)*Rho(3)*Rho(4) - 4*Rho(1)*Rho(2)*Rho(3)*Rho(3) - 2*Rho(5)*Rho(1)*Rho(2)*Rho(3) - Rho(1)*Rho(2)*Rho(4)*Rho(4) + 2*Rho(1)*Rho(2)*Rho(4) - Rho(1)*Rho(2) - Rho(1)*Rho(3)*Rho(3)*Rho(3)*Rho(3) - Rho(1)*Rho(3)*Rho(3)*Rho(4) + 2*Rho(5)*Rho(1)*Rho(3) + Rho(1)*Rho(4)*Rho(4) + Rho(1) - Rho(2)*Rho(2)*Rho(2)*Rho(2)*Rho(3) - Rho(5)*Rho(2)*Rho(2)*Rho(2)*Rho(2) + Rho(2)*Rho(2)*Rho(3)*Rho(3)*Rho(3) + Rho(2)*Rho(2)*Rho(3)*Rho(4) + Rho(2)*Rho(2)*Rho(3) + Rho(5)*Rho(2)*Rho(2)*Rho(4) + Rho(5)*Rho(2)*Rho(2) + Rho(2)*Rho(3)*Rho(3)*Rho(3) - Rho(5)*Rho(2)*Rho(3)*Rho(3) - Rho(2)*Rho(3)*Rho(4)*Rho(4) + 2*Rho(2)*Rho(3)*Rho(4) - Rho(2)*Rho(3) + Rho(3)*Rho(3)*Rho(3)*Rho(4) - Rho(3)*Rho(4) - Rho(5)*Rho(4))/((Rho(1)*Rho(1) - 2*Rho(1)*Rho(3) + Rho(3)*Rho(3) + Rho(2) + Rho(4) - Rho(2)*Rho(4) - 1)*(2*Rho(1)*Rho(3) - Rho(4) - Rho(2) - Rho(2)*Rho(4) - 4*Rho(1)*Rho(1)*Rho(2) + 2*Rho(1)*Rho(1)*Rho(4) + 3*Rho(1)*Rho(1) + 2*Rho(2)*Rho(2) + 2*Rho(2)*Rho(2)*Rho(2) + Rho(3)*Rho(3) - 4*Rho(1)*Rho(2)*Rho(3) - 1));
-           AutoregressiveParameters(2) =  (3*Rho(1)*Rho(1)*Rho(1)*Rho(1)*Rho(2) - Rho(1)*Rho(1)*Rho(1)*Rho(1)*Rho(4) - 2*Rho(1)*Rho(1)*Rho(1)*Rho(1) - Rho(1)*Rho(1)*Rho(1)*Rho(2)*Rho(3) - Rho(5)*Rho(1)*Rho(1)*Rho(1)*Rho(2) + Rho(5)*Rho(1)*Rho(1)*Rho(1)*Rho(4) + Rho(5)*Rho(1)*Rho(1)*Rho(1) - 2*Rho(1)*Rho(1)*Rho(2)*Rho(2)*Rho(2) - 2*Rho(1)*Rho(1)*Rho(2)*Rho(2) - 2*Rho(1)*Rho(1)*Rho(2)*Rho(3)*Rho(3) - 2*Rho(5)*Rho(1)*Rho(1)*Rho(2)*Rho(3) - 3*Rho(1)*Rho(1)*Rho(2)*Rho(4)*Rho(4) + 4*Rho(1)*Rho(1)*Rho(2)*Rho(4) + Rho(1)*Rho(1)*Rho(2) + 2*Rho(1)*Rho(1)*Rho(3)*Rho(3)*Rho(4) + 3*Rho(1)*Rho(1)*Rho(3)*Rho(3) - Rho(5)*Rho(1)*Rho(1)*Rho(3)*Rho(4) + Rho(1)*Rho(1)*Rho(4)*Rho(4)*Rho(4) + Rho(1)*Rho(1)*Rho(4)*Rho(4) - Rho(1)*Rho(1)*Rho(4) + Rho(1)*Rho(1) + 3*Rho(1)*Rho(2)*Rho(2)*Rho(2)*Rho(3) + Rho(5)*Rho(1)*Rho(2)*Rho(2)*Rho(2) + 3*Rho(1)*Rho(2)*Rho(2)*Rho(3)*Rho(4) - Rho(5)*Rho(1)*Rho(2)*Rho(2)*Rho(4) + 2*Rho(5)*Rho(1)*Rho(2)*Rho(2) - Rho(1)*Rho(2)*Rho(3)*Rho(3)*Rho(3) + 3*Rho(5)*Rho(1)*Rho(2)*Rho(3)*Rho(3) - Rho(1)*Rho(2)*Rho(3)*Rho(4)*Rho(4) + Rho(1)*Rho(2)*Rho(3)*Rho(4) - 3*Rho(1)*Rho(2)*Rho(3) + Rho(5)*Rho(1)*Rho(2)*Rho(4) - 2*Rho(5)*Rho(1)*Rho(2) - Rho(1)*Rho(3)*Rho(3)*Rho(3)*Rho(4) - Rho(1)*Rho(3)*Rho(3)*Rho(3) - Rho(1)*Rho(3)*Rho(4)*Rho(4) - 3*Rho(1)*Rho(3)*Rho(4) + Rho(1)*Rho(3) - Rho(5)*Rho(1)*Rho(4) - Rho(2)*Rho(2)*Rho(2)*Rho(3)*Rho(3) - Rho(5)*Rho(2)*Rho(2)*Rho(2)*Rho(3) + 2*Rho(2)*Rho(2)*Rho(2)*Rho(4)*Rho(4) - 4*Rho(2)*Rho(2)*Rho(2)*Rho(4) + 2*Rho(2)*Rho(2)*Rho(2) - 2*Rho(2)*Rho(2)*Rho(3)*Rho(3)*Rho(4) + Rho(2)*Rho(2)*Rho(3)*Rho(3) - Rho(5)*Rho(2)*Rho(2)*Rho(3) + Rho(2)*Rho(3)*Rho(3)*Rho(3)*Rho(3) + Rho(2)*Rho(3)*Rho(3)*Rho(4) + Rho(5)*Rho(2)*Rho(3)*Rho(4) - Rho(2)*Rho(4)*Rho(4)*Rho(4) + Rho(2)*Rho(4)*Rho(4) + Rho(2)*Rho(4) - Rho(2) - Rho(5)*Rho(3)*Rho(3)*Rho(3) + Rho(3)*Rho(3)*Rho(4)*Rho(4) + Rho(5)*Rho(3))/((Rho(2) + Rho(4) - 2*Rho(1)*Rho(3) + Rho(2)*Rho(4) + 4*Rho(1)*Rho(1)*Rho(2) - 2*Rho(1)*Rho(1)*Rho(4) - 3*Rho(1)*Rho(1) - 2*Rho(2)*Rho(2) - 2*Rho(2)*Rho(2)*Rho(2) - Rho(3)*Rho(3) + 4*Rho(1)*Rho(2)*Rho(3) + 1)*(Rho(1)*Rho(1) - 2*Rho(1)*Rho(3) + Rho(3)*Rho(3) + Rho(2) + Rho(4) - Rho(2)*Rho(4) - 1));
-          AutoregressiveParameters(3) = (Rho(1)*Rho(1)*Rho(1) + Rho(5)*Rho(1)*Rho(1) - 2*Rho(1)*Rho(2) - 2*Rho(1)*Rho(3)*Rho(3) + Rho(5)*Rho(1)*Rho(3) - Rho(1)*Rho(4)*Rho(4) - Rho(1)*Rho(4) + Rho(2)*Rho(2)*Rho(3) - Rho(5)*Rho(2)*Rho(2) + 2*Rho(2)*Rho(3)*Rho(4) + Rho(2)*Rho(3) - Rho(5)*Rho(2) - Rho(3)*Rho(3)*Rho(3) + Rho(3)*Rho(4) + Rho(3))/(Rho(2) + Rho(4) - 2*Rho(1)*Rho(3) + Rho(2)*Rho(4) + 4*Rho(1)*Rho(1)*Rho(2) - 2*Rho(1)*Rho(1)*Rho(4) - 3*Rho(1)*Rho(1) - 2*Rho(2)*Rho(2) - 2*Rho(2)*Rho(2)*Rho(2) - Rho(3)*Rho(3) + 4*Rho(1)*Rho(2)*Rho(3) + 1);
-          AutoregressiveParameters(4) = (Rho(1)*Rho(1)*Rho(1)*Rho(1)*Rho(2) - Rho(1)*Rho(1)*Rho(1)*Rho(1)*Rho(4) - Rho(1)*Rho(1)*Rho(1)*Rho(1) + Rho(1)*Rho(1)*Rho(1)*Rho(2)*Rho(3) - 3*Rho(5)*Rho(1)*Rho(1)*Rho(1)*Rho(2) + 2*Rho(1)*Rho(1)*Rho(1)*Rho(3)*Rho(4) + Rho(1)*Rho(1)*Rho(1)*Rho(3) + Rho(5)*Rho(1)*Rho(1)*Rho(1)*Rho(4) + 2*Rho(5)*Rho(1)*Rho(1)*Rho(1) + 4*Rho(1)*Rho(1)*Rho(2)*Rho(2)*Rho(4) - 4*Rho(1)*Rho(1)*Rho(2)*Rho(2) - 2*Rho(1)*Rho(1)*Rho(2)*Rho(3)*Rho(3) + 2*Rho(5)*Rho(1)*Rho(1)*Rho(2)*Rho(3) - 3*Rho(1)*Rho(1)*Rho(2)*Rho(4)*Rho(4) + 3*Rho(1)*Rho(1)*Rho(2) - 2*Rho(1)*Rho(1)*Rho(3)*Rho(3)*Rho(4) - 2*Rho(1)*Rho(1)*Rho(3)*Rho(3) - Rho(5)*Rho(1)*Rho(1)*Rho(3)*Rho(4) - Rho(5)*Rho(1)*Rho(1)*Rho(3) + Rho(1)*Rho(1)*Rho(4)*Rho(4)*Rho(4) - Rho(1)*Rho(1)*Rho(4) - Rho(1)*Rho(2)*Rho(2)*Rho(2)*Rho(3) + Rho(5)*Rho(1)*Rho(2)*Rho(2)*Rho(2) + Rho(1)*Rho(2)*Rho(2)*Rho(3)*Rho(4) - 2*Rho(1)*Rho(2)*Rho(2)*Rho(3) + Rho(5)*Rho(1)*Rho(2)*Rho(2)*Rho(4) + Rho(1)*Rho(2)*Rho(3)*Rho(3)*Rho(3) + Rho(5)*Rho(1)*Rho(2)*Rho(3)*Rho(3) - 3*Rho(1)*Rho(2)*Rho(3)*Rho(4)*Rho(4) + 4*Rho(1)*Rho(2)*Rho(3) - 2*Rho(5)*Rho(1)*Rho(2)*Rho(4) + Rho(5)*Rho(1)*Rho(2) + Rho(1)*Rho(3)*Rho(3)*Rho(3)*Rho(4) + 2*Rho(1)*Rho(3)*Rho(3)*Rho(3) - Rho(5)*Rho(1)*Rho(3)*Rho(3) + 3*Rho(1)*Rho(3)*Rho(4)*Rho(4) - 2*Rho(1)*Rho(3) - Rho(5)*Rho(1) - 2*Rho(2)*Rho(2)*Rho(2)*Rho(2)*Rho(4) + 2*Rho(2)*Rho(2)*Rho(2)*Rho(2) + Rho(2)*Rho(2)*Rho(2)*Rho(3)*Rho(3) - Rho(5)*Rho(2)*Rho(2)*Rho(2)*Rho(3) + 2*Rho(2)*Rho(2)*Rho(3)*Rho(3)*Rho(4) - Rho(2)*Rho(2)*Rho(3)*Rho(3) - Rho(5)*Rho(2)*Rho(2)*Rho(3) + 3*Rho(2)*Rho(2)*Rho(4)*Rho(4) - 2*Rho(2)*Rho(2)*Rho(4) - Rho(2)*Rho(2) - Rho(2)*Rho(3)*Rho(3)*Rho(3)*Rho(3) - 2*Rho(2)*Rho(3)*Rho(3)*Rho(4) + Rho(2)*Rho(3)*Rho(3) + Rho(5)*Rho(2)*Rho(3) - Rho(3)*Rho(3)*Rho(4) + Rho(5)*Rho(3)*Rho(4) - Rho(4)*Rho(4)*Rho(4) + Rho(4))/((2*Rho(1)*Rho(3) - Rho(4) - Rho(2) - Rho(2)*Rho(4) - 4*Rho(1)*Rho(1)*Rho(2) + 2*Rho(1)*Rho(1)*Rho(4) + 3*Rho(1)*Rho(1) + 2*Rho(2)*Rho(2) + 2*Rho(2)*Rho(2)*Rho(2) + Rho(3)*Rho(3) - 4*Rho(1)*Rho(2)*Rho(3) - 1)*(Rho(1)*Rho(1) - 2*Rho(1)*Rho(3) + Rho(3)*Rho(3) + Rho(2) + Rho(4) - Rho(2)*Rho(4) - 1));
-          AutoregressiveParameters(5) = (Rho(1)*Rho(1)*Rho(1)*Rho(1)*Rho(1) - 3*Rho(1)*Rho(1)*Rho(1)*Rho(1)*Rho(3) + Rho(5)*Rho(1)*Rho(1)*Rho(1)*Rho(1) + 3*Rho(1)*Rho(1)*Rho(1)*Rho(2)*Rho(2) - 6*Rho(1)*Rho(1)*Rho(1)*Rho(2)*Rho(4) - 4*Rho(1)*Rho(1)*Rho(1)*Rho(2) + 2*Rho(1)*Rho(1)*Rho(1)*Rho(3)*Rho(3) - 2*Rho(5)*Rho(1)*Rho(1)*Rho(1)*Rho(3) + Rho(1)*Rho(1)*Rho(1)*Rho(4)*Rho(4) + 4*Rho(1)*Rho(1)*Rho(1)*Rho(4) + 5*Rho(1)*Rho(1)*Rho(2)*Rho(2)*Rho(3) - 2*Rho(5)*Rho(1)*Rho(1)*Rho(2)*Rho(2) + 4*Rho(1)*Rho(1)*Rho(2)*Rho(3)*Rho(4) - 2*Rho(1)*Rho(1)*Rho(2)*Rho(3) + 4*Rho(5)*Rho(1)*Rho(1)*Rho(2) + Rho(1)*Rho(1)*Rho(3)*Rho(3)*Rho(3) + Rho(5)*Rho(1)*Rho(1)*Rho(3)*Rho(3) - Rho(1)*Rho(1)*Rho(3)*Rho(4)*Rho(4) - 2*Rho(1)*Rho(1)*Rho(3)*Rho(4) + 3*Rho(1)*Rho(1)*Rho(3) - 3*Rho(5)*Rho(1)*Rho(1) - 2*Rho(1)*Rho(2)*Rho(2)*Rho(2)*Rho(2) + 2*Rho(1)*Rho(2)*Rho(2)*Rho(2)*Rho(4) - 2*Rho(1)*Rho(2)*Rho(2)*Rho(2) - 5*Rho(1)*Rho(2)*Rho(2)*Rho(3)*Rho(3) - 2*Rho(5)*Rho(1)*Rho(2)*Rho(2)*Rho(3) + Rho(1)*Rho(2)*Rho(2)*Rho(4)*Rho(4) + 3*Rho(1)*Rho(2)*Rho(2) + 2*Rho(1)*Rho(2)*Rho(3)*Rho(3)*Rho(4) - 4*Rho(1)*Rho(2)*Rho(3)*Rho(3) + 4*Rho(5)*Rho(1)*Rho(2)*Rho(3) - 2*Rho(1)*Rho(2)*Rho(4)*Rho(4) + 2*Rho(1)*Rho(2)*Rho(4) - Rho(1)*Rho(3)*Rho(3)*Rho(3)*Rho(3) - 2*Rho(1)*Rho(3)*Rho(3)*Rho(4) + Rho(1)*Rho(3)*Rho(3) - 2*Rho(1)*Rho(4) + Rho(2)*Rho(2)*Rho(2)*Rho(2)*Rho(3) + Rho(5)*Rho(2)*Rho(2)*Rho(2)*Rho(2) - 2*Rho(2)*Rho(2)*Rho(2)*Rho(3)*Rho(4) + 2*Rho(2)*Rho(2)*Rho(2)*Rho(3) + Rho(2)*Rho(2)*Rho(3)*Rho(3)*Rho(3) - 2*Rho(2)*Rho(2)*Rho(3)*Rho(4) - 2*Rho(5)*Rho(2)*Rho(2) + 2*Rho(2)*Rho(3)*Rho(3)*Rho(3) + 2*Rho(2)*Rho(3)*Rho(4) - 2*Rho(2)*Rho(3) - Rho(5)*Rho(3)*Rho(3) + Rho(3)*Rho(4)*Rho(4) + Rho(5))/((Rho(1)*Rho(1) - 2*Rho(1)*Rho(3) + Rho(3)*Rho(3) + Rho(2) + Rho(4) - Rho(2)*Rho(4) - 1)*(2*Rho(1)*Rho(3) - Rho(4) - Rho(2) - Rho(2)*Rho(4) - 4*Rho(1)*Rho(1)*Rho(2) + 2*Rho(1)*Rho(1)*Rho(4) + 3*Rho(1)*Rho(1) + 2*Rho(2)*Rho(2) + 2*Rho(2)*Rho(2)*Rho(2) + Rho(3)*Rho(3) - 4*Rho(1)*Rho(2)*Rho(3) - 1));
       otherwise
         AutocorrelationMatrix = eye(p);
         for i=1:p-1