function [f,g] = exafun(u);
% Example function for testing the optimization routine grad_method.m
f = 0.5*norm(u-1)^2;
g = u - 1;
