請教MATLAB編程,？

展翅翱翔with

優(yōu)化設(shè)計老師布置用MATLAB對變尺度法進行MATLAB編程（也叫擬牛頓法），我是用擬牛頓法中的BFD法進行編程，程序如下：（問題在程序后面）
: z% y" z5 x, j" M' E5 Y: l. n. _6 Jfunction [r,n]=mulDFP(F,x0,B,eps)
%用擬牛頓法中的DFP法中求極小點和極小值。
%Designed by GAO,Mechanical Engineering College
$ i8 i& g7 z. N8 j7 k" B%Shenyang University of Technology
* g& H/ z' I, x! u4 [- c%November,2015
%調(diào)用格式：
%syms x1 x2；
3 d! X) ^+ f' }1 Y%z=表達式;
%zx1=simple(diff(z,x1)),zx2=diff(z,x2)
2 ]& V7 U) w. u# b# W. b%X0=[給定初始值];
%f=[zx1;zx2];
$ V, c) B4 V5 L4 ?. L%[n,r]=mulDFP(f,x0)
/ G. @3 a3 B9 i: [. a, `if nargin==2
l = length(x0);
2 J' J% X% q7 L) }1 F$ aB=eye(l); %A取為單位陣
eps=1.0e-4;
else
if nargin==3
3 l# p5 q% V' a. oeps=1.0e-4;
6 j( C1 i7 E3 f* Y5 d6 s' U0 Kend
  r4 }1 c; N: n% S# f6 A( c0 ?( |end
fx = subs(F,findsym(F),x0);
" a- S0 S. r% U9 v. A1 H, @r=transpose(x0)-B*fx;
n=1;
tol=1;
while tol>eps
x0=r;
fx = subs(F,findsym(F),x0);
1 Y5 b# ~( X( o# ir=x0-B*fx;
3 T- n* h3 N* @; K4 v4 Fy=r-x0;
fr = subs(F,findsym(F),r);
) d; J% ^' h3 W0 h" g# qz = fr-fx;
9 E" R: ?% Q! KB1=B+ y*y'/(y'*z)-B*z*z'*B/(z'*B*z); %調(diào)整A
( v) u' J& Z- l. YB=B1;
n=n+1;
$ s8 Z+ x" K. qif(n>100000) %迭代步數(shù)控制
disp('迭代步數(shù)太多,，可能不收斂,！');
return;
end
& L+ H2 c: v; ptol=norm(r-x0);
end
問題：這個程序我是通過其他方法改變來的，現(xiàn)在想在這個基礎(chǔ)上增加一個求偏導(dǎo)的一個函數(shù)類似于“df=jacobian(f,[x1 x2]);%函數(shù)f的偏導(dǎo)”,，可是加進去以后就會出現(xiàn)錯誤,，改正好多次，都不行,，后來實在不行我就加在了調(diào)用函數(shù)中,，可是這么做，運行效率太低,，唯恐老師會減分,，所以希望各位前輩能給予實質(zhì)性的指點，謝謝,！