请教MATLAB编程？

展翅翱翔with

优化设计老师布置用MATLAB对变尺度法进行MATLAB编程（也叫拟牛顿法），我是用拟牛顿法中的BFD法进行编程，程序如下：（问题在程序后面）
1 C5 \' t8 O5 z8 x+ d, sfunction [r,n]=mulDFP(F,x0,B,eps)
' x. `5 x$ |7 L7 ]! J$ M%用拟牛顿法中的DFP法中求极小点和极小值。
. |% Q% w! D6 R) m%Designed by GAO,Mechanical Engineering College
) I3 b( K; O6 Q* [* `%Shenyang University of Technology
) A9 A9 O8 N3 W4 b3 w7 w& C%November,2015
8 Z7 M3 ]/ F, C" h4 f$ ^0 A" ~%调用格式：
%syms x1 x2；
% I3 m5 v, d1 Y# K3 W%z=表达式;
8 {( F5 c, z5 L%zx1=simple(diff(z,x1)),zx2=diff(z,x2)
; T( \4 [1 S4 P9 Z6 I9 y6 S%X0=[给定初始值];
2 @" V! r: c9 i( C; _( R& L* y: \; d  m%f=[zx1;zx2];
- x" K( U4 T) D4 T6 a/ X%[n,r]=mulDFP(f,x0)
if nargin==2
l = length(x0);
B=eye(l); %A取为单位阵
eps=1.0e-4;
4 c4 i( U2 J7 X. s. V$ X* ielse
if nargin==3
eps=1.0e-4;
. \: h/ m' n0 {+ k& b: eend
end
fx = subs(F,findsym(F),x0);
# @4 I2 c/ R) l" \+ j/ m& {! o- l# ]r=transpose(x0)-B*fx;
n=1;
tol=1;
while tol>eps
x0=r;
fx = subs(F,findsym(F),x0);
r=x0-B*fx;
y=r-x0;
fr = subs(F,findsym(F),r);
3 C  x' _" Q! Q4 Z# ez = fr-fx;
B1=B+ y*y'/(y'*z)-B*z*z'*B/(z'*B*z); %调整A
8 ]; O$ @9 Q7 ]$ y( XB=B1;
n=n+1;
" o% x5 y) D# |( X" l- yif(n>100000) %迭代步数控制
disp('迭代步数太多，可能不收敛！');
return;
$ w& I- v! l' Fend
  z. J! |& T. U7 E& t+ |+ Stol=norm(r-x0);
end
问题：这个程序我是通过其他方法改变来的，现在想在这个基础上增加一个求偏导的一个函数类似于“df=jacobian(f,[x1 x2]);%函数f的偏导”，可是加进去以后就会出现错误，改正好多次，都不行，后来实在不行我就加在了调用函数中，可是这么做，运行效率太低，唯恐老师会减分，所以希望各位前辈能给予实质性的指点，谢谢！