9. Anexos.

Código. Preparación de las restricciones.
9.1.3.2. La llamada a fmincon y la comprobación de resultados.
Anexos.
El siguiente fragmento de código establece todas las opciones a utilizar en la llamada a la función
fmincon, la propia llamada y la comprobación de que la solución cumple las restricciones.
options=optimset('GradObj','On','GradConstr','On','LargeScale','on','Display','iter');
options = optimset(options,'MaxIter',500,'MaxFunEvals',10000);
[x,fval,exitflag,output,lambda,grad]=fmincon(@fun_objetivo,Xij0,[],[],
Aeq,beq,lb,ub,@nlincon,options,VectorDir(:,3),NumTerminales)
[C,Ceq,DC,DCeq]=nlincon(x,VectorDir(:,3),NumTerminales)
for i=1:1:size(x)/2
end
if x(i)>0 & x(i)0 & x(i+size(x,1)/2) 0 & x(i+size(x,1)/2)>0)
pos_error=i;
error=1;
elseif x(i)>0 | x(i+size(x,1)/2)>0
end
CosteSolucion=CosteSolucion + VectorDir(i,3);
if error==1
end
fprintf(1,'Se ha incumplido una restriccion: Xij*Xji=0-->VectorDir(%u)',pos_error);
fid = fopen('resultados.txt','a+');
fprintf(fid,'%s\t\t%u\t%f\t%u\t\t%u \n',archivo,origen,fval,exitflag,CosteSolucion);
fclose(fid);
Código. Opciones, llamada a fmincon y comprobación del resultado.
13