function f = fnctlin(x); % Evaluates 1/2 x^T my_A x + x^T my_c. global my_A my_c f = .5 * (x' * (my_A*x)) + x'*my_c;