请教MATLAB编程？

展翅翱翔with

优化设计老师布置用MATLAB对变尺度法进行MATLAB编程（也叫拟牛顿法），我是用拟牛顿法中的BFD法进行编程，程序如下：（问题在程序后面）
function [r,n]=mulDFP(F,x0,B,eps)
1 q1 q9 F. f6 e3 {& n%用拟牛顿法中的DFP法中求极小点和极小值。
%Designed by GAO,Mechanical Engineering College
%Shenyang University of Technology
%November,2015
3 G) [4 M8 x3 z' q- z9 L- I9 x) r%调用格式：
%syms x1 x2；
%z=表达式;
- V! _' _: _, n: u, R% j%zx1=simple(diff(z,x1)),zx2=diff(z,x2)
5 w2 x3 A" a, o/ n* W$ a%X0=[给定初始值];
%f=[zx1;zx2];
%[n,r]=mulDFP(f,x0)
; E, t, E5 A  D4 D# K5 G0 K9 Jif nargin==2
* q9 m& G% n9 V# E; h) [% wl = length(x0);
, z9 s* P3 g, wB=eye(l); %A取为单位阵
0 }: Y; `& X% P+ b7 V/ F, m. ^eps=1.0e-4;
else
if nargin==3
eps=1.0e-4;
end
end
fx = subs(F,findsym(F),x0);
2 u, a8 i  v/ [, Rr=transpose(x0)-B*fx;
n=1;
* o! i1 W: u/ itol=1;
& S. j7 ?1 l3 l" ~2 A/ zwhile tol>eps
x0=r;
: ^5 n# E, [0 s2 r$ Hfx = subs(F,findsym(F),x0);
$ t" l5 Z4 n" Z) Ur=x0-B*fx;
( w9 V$ b8 e% `( l) ]5 Jy=r-x0;
fr = subs(F,findsym(F),r);
z = fr-fx;
) m( K& f3 y" K. A$ C$ X' iB1=B+ y*y'/(y'*z)-B*z*z'*B/(z'*B*z); %调整A
+ T! ]" |& |$ T5 l, n/ S: fB=B1;
n=n+1;
if(n>100000) %迭代步数控制
disp('迭代步数太多，可能不收敛！');
return;
end
& z% R/ N. \4 _& W1 [tol=norm(r-x0);
end
0 ^2 _$ ?* t0 N4 p5 [问题：这个程序我是通过其他方法改变来的，现在想在这个基础上增加一个求偏导的一个函数类似于“df=jacobian(f,[x1 x2]);%函数f的偏导”，可是加进去以后就会出现错误，改正好多次，都不行，后来实在不行我就加在了调用函数中，可是这么做，运行效率太低，唯恐老师会减分，所以希望各位前辈能给予实质性的指点，谢谢！