codice matlab

[barco_]

Salve a tutti, sono da poco iscritto a questo forum. Oltre ad essere nuovo del forum ho cominciato da poco a studiare matlab, benchè il mio background sia di ricerca operativa. Vorrei chiedervi un aiuto nel decifrare il seguente codice in un modello LP. Ringrazio in anticipo chiunque voglia darmi una mano, anche solo per parte di esso. Ciao!!

function [x,y1,y2,z]=P18(T1)

c1=T1(1:3,3);
c2=T1(1:3,4);
c3=zeros(1,2);
c=[c1',c1',c2',c2',c3];
d=ones(1,14);

O6=ones(1,6);
Z6=zeros(1,6);
Z12=zeros(4,12);
I2=eye(2);
NI2=I2*-1;
I6=eye(6);
Z62=zeros(6,2);
Min=NI2*T1(4,1);
Max=NI2*T1(5,1);
O3=ones(1,3);
Z3=zeros(1,3);

A1=[I6,I6,Z62];
A2=[O3,Z3,Z6;Z6,O3,Z3];
A3=[A2,Min;A2,Max];
A4=[I2;I2];
A5=[Z12,A4];
A6=[O6,Z6;Z6,O6];
A7=[A6,NI2];
A=[A7;A5;A3;A1];

b1=zeros(2,1);
for i=1:2
for j=1:2
b2=[b2;T1(3+j,2+i)];
end
end
b3=zeros(4,1);
for i=1:2
for j=1:3
b4=[b4;T1(j,i)];
end
end
b=[b1;b2;b3;b4];
e3=ones(2,1)*3;
e2=ones(2,1)*2;
e1=ones(2,1);
e=[e3;e2;e1;e2;e1;e3;e3;e3];

[status,x,y,z]=lp_solver(A,b,c,'min',e,d);