Skip to content
Snippets Groups Projects
Verified Commit c5c13077 authored by Sébastien Villemot's avatar Sébastien Villemot
Browse files

Minor simplification

parent 36ab8d00
No related merge requests found
Pipeline #1177 passed
......@@ -29,7 +29,7 @@ steady;
stoch_simul(order=2,irf=0);
sigma2=M_.Sigma_e;
i = linspace(1,800,800);
i = 1:800;
a = theta*xbar*i+(theta^2*sigma2)/(2*(1-rho)^2)*(i-2*rho*(1-rho.^i)/(1-rho)+rho^2*(1-rho.^(2*i))/(1-rho^2));
a1 = theta^2*sigma2/(2*(1-rho)^2)*(i-2*rho*(1-rho.^i)/(1-rho)+rho^2*(1-rho.^(2*i))/(1-rho^2));
b = theta*rho*(1-rho.^i)/(1-rho);
......@@ -46,4 +46,4 @@ x2 = [
];
if any(abs(x1-x2) > 1e-14);
error('burnside_1 doesn''t reproduce the analytical solution');
end;
\ No newline at end of file
end;
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment