frank_wolfe algorithm for order 3

frank_wolfe algorithm for order 3

e=0.005;
X0=[0.1;0.2;0.7];
syms x1 x2 x3;
f=0.5*(x1^2+x2^2+0.1*x3^2)+0.55*x3;
[X,fval]=frank_wolfe(f,X0,e,[],[],[1 1 1],[1],[0;0;0],[1;1;1])