function [ jval,gradient ] = CostFunction( theta )
%UNTITLED2 Summary of this function goes here
% Detailed explanation goes here
jval=(theta(1)-5)^2+(theta(2)-5)^2;
gradient=zeros(2,1);
gradient(1)=2*(theta(1)-5);
gradient(2)=2*(theta(2)-5);
end
options=optimset('GradObj','on','MaxIter',100);
initiatheta=zeros(2,1);
[optheta,functionVal,exitFlag]=fminunc(@CostFunction,initiatheta,options);