   function [X,L,G,rr]=rikk;
   A=[0 1;-1 1];
   B=[0 0;0 1];
   R=[2 0;0 2];
   Q=[1 0;0 1];
   [X,L,G,rr]=dare(A,B,Q,R);
   end;