请教MATLAB编程？

展翅翱翔with

优化设计老师布置用MATLAB对变尺度法进行MATLAB编程（也叫拟牛顿法），我是用拟牛顿法中的BFD法进行编程，程序如下：（问题在程序后面）
function [r,n]=mulDFP(F,x0,B,eps)
) }  @7 w! C8 z3 e! D7 `%用拟牛顿法中的DFP法中求极小点和极小值。
%Designed by GAO,Mechanical Engineering College
  E) a- L1 D/ _3 |%Shenyang University of Technology
%November,2015
" Z( k1 d8 o% K! O4 h%调用格式：
+ G7 Q$ C/ w0 B4 v7 v%syms x1 x2；
%z=表达式;
%zx1=simple(diff(z,x1)),zx2=diff(z,x2)
%X0=[给定初始值];
( ~& P2 \9 c3 W%f=[zx1;zx2];
; Q0 S9 Z  i" M: e4 O4 e8 v1 A%[n,r]=mulDFP(f,x0)
if nargin==2
l = length(x0);
* y4 C8 y: W4 g3 ]4 AB=eye(l); %A取为单位阵
eps=1.0e-4;
else
, @, }6 E7 Y* x4 Q- ^, pif nargin==3
eps=1.0e-4;
4 N4 M: \% `+ \4 S8 Fend
end
fx = subs(F,findsym(F),x0);
. q5 o' i, L& ~! s6 a/ l* z7 `r=transpose(x0)-B*fx;
( \) h: x7 d# I  x- i  hn=1;
tol=1;
while tol>eps
  s( K9 t' j* I4 I0 A% ?/ U1 Hx0=r;
# X0 B( R4 m/ k& B; J# W$ ?fx = subs(F,findsym(F),x0);
r=x0-B*fx;
y=r-x0;
fr = subs(F,findsym(F),r);
6 p+ d& d2 y  r0 D' m) G* b4 Cz = fr-fx;
% g) }) `+ i- h" g% C. h% {B1=B+ y*y'/(y'*z)-B*z*z'*B/(z'*B*z); %调整A
B=B1;
n=n+1;
, G( Z- }9 J: o+ N7 c) ?if(n>100000) %迭代步数控制
disp('迭代步数太多，可能不收敛！');
return;
end
4 {$ [5 ^0 B. X2 S( G$ xtol=norm(r-x0);
end
问题：这个程序我是通过其他方法改变来的，现在想在这个基础上增加一个求偏导的一个函数类似于“df=jacobian(f,[x1 x2]);%函数f的偏导”，可是加进去以后就会出现错误，改正好多次，都不行，后来实在不行我就加在了调用函数中，可是这么做，运行效率太低，唯恐老师会减分，所以希望各位前辈能给予实质性的指点，谢谢！