Skip to content
Snippets Groups Projects
Commit 32de0fd6 authored by Marco Ratto's avatar Marco Ratto
Browse files

fixed typos

parent f4555573
No related branches found
No related tags found
No related merge requests found
function x=sylvester3mr(a,b,c,d)
function x=sylvester3(a,b,c,d)
% solves a*x+b*x*c=d where d is [n x m x p]
% Copyright (C) 2005-2009 Dynare Team
......@@ -21,6 +21,7 @@ function x=sylvester3mr(a,b,c,d)
n = size(a,1);
m = size(c,1);
p = size(d,3);
x=zeros(n,m,p);
if n == 1
for j=1:p,
x(:,:,j)=d(:,:,j)./(a*ones(1,m)+b*c);
......@@ -33,7 +34,6 @@ if m == 1
end
return;
end
x=zeros(n,m,p);
[u,t]=schur(c);
if exist('OCTAVE_VERSION')
[aa,bb,qq,zz]=qz(full(a),full(b));
......
......@@ -32,7 +32,7 @@ for j=1:size(dd,3),
iter = iter + 1;
end
if iter == 500
sprintf('sylvester3amr : Only accuracy of %10.8f is achieved after 500 iterations',e);
sprintf('sylvester3a : Only accuracy of %10.8f is achieved after 500 iterations',e);
flag=1;
end
end
\ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment