請教MATLAB編程？

展翅翱翔with

優(yōu)化設(shè)計老師布置用MATLAB對變尺度法進行MATLAB編程（也叫擬牛頓法），我是用擬牛頓法中的BFD法進行編程，程序如下：（問題在程序后面）
function [r,n]=mulDFP(F,x0,B,eps)
3 V! x. w8 i) l2 q3 N%用擬牛頓法中的DFP法中求極小點和極小值。
%Designed by GAO,Mechanical Engineering College
%Shenyang University of Technology
%November,2015
%調(diào)用格式：
' M% o1 \  ~4 Y& a) |; N%syms x1 x2；
& E1 z* G+ K0 S: f* ?" p%z=表達式;
%zx1=simple(diff(z,x1)),zx2=diff(z,x2)
%X0=[給定初始值];
' g& |8 e% H3 I: i5 G/ P%f=[zx1;zx2];
%[n,r]=mulDFP(f,x0)
2 G# u. t* a# Pif nargin==2
l = length(x0);
. K7 r& O2 k2 ?. JB=eye(l); %A取為單位陣
eps=1.0e-4;
8 [/ j4 ], i" n% L% Kelse
if nargin==3
eps=1.0e-4;
2 j1 e- j  k1 S$ c4 gend
: K9 i3 T) V; m, @) ~end
0 h- m8 h1 f/ {1 mfx = subs(F,findsym(F),x0);
r=transpose(x0)-B*fx;
5 \  b0 O3 J# }# O9 Q) [n=1;
tol=1;
$ Q& J6 `; Y7 [0 s6 zwhile tol>eps
/ T" T; I2 m- \" N1 Y( Px0=r;
9 b. ]8 x8 m7 [% R. R3 nfx = subs(F,findsym(F),x0);
r=x0-B*fx;
y=r-x0;
fr = subs(F,findsym(F),r);
z = fr-fx;
B1=B+ y*y'/(y'*z)-B*z*z'*B/(z'*B*z); %調(diào)整A
B=B1;
8 [) W2 ]& p* p# T2 Vn=n+1;
if(n>100000) %迭代步數(shù)控制
disp('迭代步數(shù)太多，可能不收斂！');
return;
! j8 U7 q. K8 _) H  a9 k( iend
tol=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ì)性的指點，謝謝！