function err=optcontr2(n,m,alpha)
A = (n+1)*(n+1)*(sparse(2:n,1:n-1,-1,n,n) + sparse(1:n,1:n,2,n,n) + sparse(1:n-1,2:n,-1,n,n));
matobs = zeros(n,n);
for i=1:m
  index=int8((i-0.5)/m*(n+1));
  matobs(index,index)=1.;
end
mat1 = horzcat(???,-speye(n));
mat2 = horzcat(???,alpha*A');
mat  = vertcat(mat1,mat2);
xobs=A\ones(n,1);
b = ones(n,1);
rhs = vertcat(b,???);
sol=mat\rhs;
x=sol(1:end/2);
q=sol(end/2+1:end);
plot(x,'o');
hold on;
plot(xobs,'+');
plot(q,'x');
err=max(abs(x-xobs));
end
