% Solution proposal Task 13
% 

% generate some data matrices Y and U to be used for testing purposes
A=0.9; B=0.5; D=1; E=-1;     % Model parameters

N=200;
U=prbs1(N,10,40);
Y=dsrsim(A,B,D,E,U);
ny=size(Y,2);

% Parameters for DSR
L=2; g=1; J=2;

[a,b,d,e,c,f,x0]=dsr(Y,U,L,g,J);

%[a,b,d,e,x0]=mydsr(Y,U,L,g,J);