首先是给出costFunction, 然后设置参数的优化选项，以及参数的初始值，最后是调用fminunc函数
代码如下：
function [jVal, gradient]=costFunction(theta)
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);
initialTheta=zeros(2,1)
[optTheta, functionVal, exitFlag]=fminunc(@costFunction, initialTheta, options)
