این چیزی هست که در کد شما اضافه کردم
clc
clear
close all
syms f x1 x2 x3 k m p l o
x0=[0, -1, 1]
tol=0.1;
f=(2*x1^2)-(3*x1*x3)+(4*x2^3)+(3*x2*x1)+3
m=eval(x0+k*(-gradient(f,[x1,x2,x3])')*(inv(hessian(f))))
z=subs(f,[x1,x2,x3],m)
y=subs(z,[x1,x2,x3],x0)
p=diff(y)
u=solve (p)
k=u...