function ps33_6242_2004(n,m,M,R,K)
% function ps33_6242_2004(n,m,M,R,K)
%
% Solves Problem 3.3

if nargin<1, n=50; end
if nargin<2, m=10; end
if nargin<3, M=1; end
if nargin<4, R=0.2; end
if nargin<5, K=4; end

g=(2*(n^2)*K/M)*toeplitz([ 2 -1 zeros(1,2*n-1)]);
In=eye(2*n+1);
On=zeros(2*n+1);
b=[1;zeros(2*n,1)];
c=[zeros(1,n) 1 zeros(1,n)];
A=[On In; -g -(R/M)*In];
B=[zeros(2*n+1,1);b];
C=[c zeros(1,2*n+1)];

G=ss(A,B,C,0);
H=hsvd_6242(G,m);