function y = grhs(x) % % y = x; y(1) = 2*x(1); y(2) = 2*x(2);