Steepest Descent

Ska beräkna lokala extrempunkter till funktionen:

$f={x_1}^3+0.3{x_1}^2+3x_1{x}_2+0.9{x_2}^3-x_2$

Använder steepest descent-metoden på -f med startapproximationen (-1,-1) för att hitta den maxpunkt som ligger nära (-1,-1).

Har fått en del hjälpkod som jag måste justera/lägga till i och sen avläsa vissa värden. Jag får inte själva koden att köra ens dock, är väldigt förvirrad. Infogar det jag har.

sd.m

function [x,fmin,k]=sd(f,x0,tol)

x=x0; kmax=750;

for k=1;kmax;

    grad=jacobi(f,x); 

    g=@(s)f(x-s*grad'); 

    alpha=fminbnd(g,0,1); 

    x=x-alpha*grad'; %stegar fram till nästa punkt

    if norm(alpha*grad)<tol, break, end 

end

fmin=f(x); 

jacobi.m

function A = jacobi(f,x)

h=1e-6; n=length(x);

for i=1:n

    x1=x; x1(i)=x1(i)+h;

    x2=x; x2(i)=x2(i)-h;

    A(:,i)=(f(x1)-f(x2))/(2*h);

end

f.m

function y = f(x)

y = x(1).^3+0.3*x(1).^2+3*x(1).*x(2)+0.9*x(2).^3-x(2);

......

Förstår inte hur alla dessa ska gå ihop, eller hur jag ska hämta min funktion f, i samtliga scripts får jag felmeddelanden för att input arguments saknas.