Advanced Search

Question: Schurform call from within a procedure

Question:Schurform call from within a procedure

Johan Sjoberg 40 Maple
complex package
February 19 2007
0
Hi, I have a small problem. I want to use the Schurform computation in the LinearAlgebra-package. The main objective is to compute the stabilizing solution to a Algebraic Matrix Riccati equation. (perhaps someone knows some code for this?) I call the SchurForm function as (T,Z) := SchurForm(A,output=['T','Z']) where A has the correct type (ie complex(sfloat)). This works well in a Maple sheet, but in a procedure (where T and Z are defined local) the kernel says that the output chosen is invalid. If I don't specify any output, it also works good (but then I get no Z). The Maple version is 10 and I am using Linux. Thanks in advance! Best regards Johan Sjöberg
Please Wait...