%punti_di_equilibrio con disturbi costanti e matrici linearizzate

ueq=[volume1/(60*6) volume2/(60*6 ) volume3/(60*6) volume4/(60*6) volume5/(60*6) volume6/(60*6) volume7/(60*6) volume8/(60*6) volume9/(60*6) volume10/(60*6) volume11/(60*6) volume12/(60*6) volume13/(60*6) volume14/(60*6) volume15/(60*6) volume16/(60*6) volume17/(60*6) volume18/(60*6) volume19/(60*6) volume20/(60*6) volume21/(60*6) volume22/(60*6) volume23/(60*6) volume24/(60*6) volume25/(60*6) volume26/(60*6) volume27/(60*6) volume28/(60*6) volume29/(60*6) volume30/(60*6) volume31/(60*6) volume32/(60*6) volume33/(60*6) %0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 %5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 %tasso di ventilazione in metri cubi s u minuto
      0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1 0.24/volume1%0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 %(22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) (22*10^-3)/(60*2) %tasso di umidifiaczione in kilogrammi su minuto
      10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 ];%1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 ]; %riscaldamento in watt
Ueq=[ueq(1,1) ueq(2,1) ueq(3,1) ueq(1,2) ueq(2,2) ueq(3,2) ueq(1,3) ueq(2,3) ueq(3,3) ueq(1,4) ueq(2,4) ueq(3,4) ueq(1,5) ueq(2,5) ueq(3,5) ueq(1,6) ueq(2,6) ueq(3,6) ueq(1,7) ueq(2,7) ueq(3,7) ueq(1,8) ueq(2,8) ueq(3,8) ueq(1,9) ueq(2,9) ueq(3,9) ueq(1,10) ueq(2,10) ueq(3,10) ueq(1,11) ueq(2,11) ueq(3,11) ueq(1,12) ueq(2,12) ueq(3,12) ueq(1,13) ueq(2,13) ueq(3,13) ueq(1,14) ueq(2,14) ueq(3,14) ueq(1,15) ueq(2,15) ueq(3,15) ueq(1,16) ueq(2,16) ueq(3,16) ueq(1,17) ueq(2,17) ueq(3,17) ueq(1,18) ueq(2,18) ueq(3,18) ueq(1,19) ueq(2,19) ueq(3,19) ueq(1,20) ueq(2,20) ueq(3,20) ueq(1,21) ueq(2,21) ueq(3,21) ueq(1,22) ueq(2,22) ueq(3,22) ueq(1,23) ueq(2,23) ueq(3,23) ueq(1,24) ueq(2,24) ueq(3,24) ueq(1,25) ueq(2,25) ueq(3,25) ueq(1,26) ueq(2,26) ueq(3,26) ueq(1,27) ueq(2,27) ueq(3,27) ueq(1,28) ueq(2,28) ueq(3,28) ueq(1,29) ueq(2,29) ueq(3,29) ueq(1,30) ueq(2,30) ueq(3,30) ueq(1,31) ueq(2,31) ueq(3,31) ueq(1,32) ueq(2,32) ueq(3,32) ueq(1,33) ueq(2,33) ueq(3,33)]';
%xeq=[10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 10 0.105332	0.1053333	0.053333	0.053326	0.053333	0.053333	0.05332	0.053438	0.053351	0.053369	0.053439	0.053404 0.053369 0.053333 0.053333 0.053326 0.053333 0.053281 0.053473	0.053369 0.053333 0.053475 0.053403 0.053333 0.053311 0.053322 0.053333	0.050086 0.053333 0.053422 0.053333 0.053333 0.053374]; %vettore stati di equilibrio  
xeq=[293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 293.5 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3 22*10^-3];
Xeq=xeq';
% matrice A
  
A = [-ueq(1,1)/volume1-K_est*(somma1-pareti_esterne1(1,1))/cost1-K_int*(sum(area1))/cost1 ad(1,2)*K_int/cost1 ad(1,3)*K_int/cost1 ad(1,4)*K_int/cost1 ad(1,5)*K_int/cost1 ad(1,6)*K_int/cost1 ad(1,7)*K_int/cost1 ad(1,8)*K_int/cost1 ad(1,9)*K_int/cost1 ad(1,10)*K_int/cost1 ad(1,11)*K_int/cost1 ad(1,12)*K_int/cost1 ad(1,13)*K_int/cost1 ad(1,14)*K_int/cost1 ad(1,15)*K_int/cost1 ad(1,16)*K_int/cost1 ad(1,17)*K_int/cost1 ad(1,18)*K_int/cost1 ad(1,19)*K_int/cost1 ad(1,20)*K_int/cost1 ad(1,21)*K_int/cost1 ad(1,22)*K_int/cost1 ad(1,23)*K_int/cost1 ad(1,24)*K_int/cost1 ad(1,25)*K_int/cost1 ad(1,26)*K_int/cost1 ad(1,27)*K_int/cost1 ad(1,28)*K_int/cost1 ad(1,29)*K_int/cost1 ad(1,30)*K_int/cost1 ad(1,31)*K_int/cost1 ad(1,32)*K_int/cost1 ad(1,33)*K_int/cost1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(2,1)*K_int/cost2 -ueq(1,2)/volume2-K_est*(somma2-pareti_esterne2(1,1))/cost2-K_int*(sum(area2))/cost2 ad(2,3)*K_int/cost2 ad(2,4)*K_int/cost2 ad(2,5)*K_int/cost2 ad(2,6)*K_int/cost2 ad(2,7)*K_int/cost2 ad(2,8)*K_int/cost2 ad(2,9)*K_int/cost2 ad(2,10)*K_int/cost2 ad(2,11)*K_int/cost2 ad(2,12)*K_int/cost2 ad(2,13)*K_int/cost2 ad(2,14)*K_int/cost2 ad(2,15)*K_int/cost2 ad(2,16)*K_int/cost2 ad(2,17)*K_int/cost2 ad(2,18)*K_int/cost2 ad(2,19)*K_int/cost2 ad(2,20)*K_int/cost2 ad(2,21)*K_int/cost2 ad(2,22)*K_int/cost2 ad(2,23)*K_int/cost2 ad(2,24)*K_int/cost2 ad(2,25)*K_int/cost2 ad(2,26)*K_int/cost2 ad(2,27)*K_int/cost2 ad(2,28)*K_int/cost2 ad(2,29)*K_int/cost2 ad(2,30)*K_int/cost2 ad(2,31)*K_int/cost2 ad(2,32)*K_int/cost2 ad(2,33)*K_int/cost2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(3,1)*K_int/cost3 ad(3,2)*K_int/cost3 -ueq(1,3)/volume3-K_est*(somma3-pareti_esterne3(1,1))/cost3-K_int*(sum(area3))/cost3 ad(3,4)*K_int/cost3 ad(3,5)*K_int/cost3 ad(3,6)*K_int/cost3 ad(3,7)*K_int/cost3 ad(3,8)*K_int/cost3 ad(3,9)*K_int/cost3 ad(3,10)*K_int/cost3 ad(3,11)*K_int/cost3 ad(3,12)*K_int/cost3 ad(3,13)*K_int/cost3 ad(3,14)*K_int/cost3 ad(3,15)*K_int/cost3 ad(3,16)*K_int/cost3 ad(3,17)*K_int/cost3 ad(3,18)*K_int/cost3 ad(3,19)*K_int/cost3 ad(3,20)*K_int/cost3 ad(3,21)*K_int/cost3 ad(3,22)*K_int/cost3 ad(3,23)*K_int/cost3 ad(3,24)*K_int/cost3 ad(3,25)*K_int/cost3 ad(3,26)*K_int/cost3 ad(3,27)*K_int/cost3 ad(3,28)*K_int/cost3 ad(3,29)*K_int/cost3 ad(3,30)*K_int/cost3 ad(3,31)*K_int/cost3 ad(3,32)*K_int/cost3 ad(3,33)*K_int/cost3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(4,1)*K_int/cost4 ad(4,2)*K_int/cost4 ad(4,3)*K_int/cost4 -ueq(1,4)/volume4-K_est*(somma4-pareti_esterne4(1,1))/cost4-K_int*(sum(area4))/cost4 ad(4,5)*K_int/cost4 ad(4,6)*K_int/cost4 ad(4,7)*K_int/cost4 ad(4,8)*K_int/cost4 ad(4,9)*K_int/cost4 ad(4,10)*K_int/cost4 ad(4,11)*K_int/cost4 ad(4,12)*K_int/cost4 ad(4,13)*K_int/cost4 ad(4,14)*K_int/cost4 ad(4,15)*K_int/cost4 ad(4,16)*K_int/cost4 ad(4,17)*K_int/cost4 ad(4,18)*K_int/cost4 ad(4,19)*K_int/cost4 ad(4,20)*K_int/cost4 ad(4,21)*K_int/cost4 ad(4,22)*K_int/cost4 ad(4,23)*K_int/cost4 ad(4,24)*K_int/cost4 ad(4,25)*K_int/cost4 ad(4,26)*K_int/cost4 ad(4,27)*K_int/cost4 ad(4,28)*K_int/cost4 ad(4,29)*K_int/cost4 ad(4,30)*K_int/cost4 ad(4,31)*K_int/cost4 ad(4,32)*K_int/cost4 ad(4,33)*K_int/cost4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(5,1)*K_int/cost5 ad(5,2)*K_int/cost5 ad(5,3)*K_int/cost5 ad(5,4)*K_int/cost5 -ueq(1,5)/volume5-K_est*(somma5-pareti_esterne5(1,1))/cost5-K_int*(sum(area5))/cost5 ad(5,6)*K_int/cost5 ad(5,7)*K_int/cost5 ad(5,8)*K_int/cost5 ad(5,9)*K_int/cost5 ad(5,10)*K_int/cost5 ad(5,11)*K_int/cost5 ad(5,12)*K_int/cost5 ad(5,13)*K_int/cost5 ad(5,14)*K_int/cost5 ad(5,15)*K_int/cost5 ad(5,16)*K_int/cost5 ad(5,17)*K_int/cost5 ad(5,18)*K_int/cost5 ad(5,19)*K_int/cost5 ad(5,20)*K_int/cost5 ad(5,21)*K_int/cost5 ad(5,22)*K_int/cost5 ad(5,23)*K_int/cost5 ad(5,24)*K_int/cost5 ad(5,25)*K_int/cost5 ad(5,26)*K_int/cost5 ad(5,27)*K_int/cost5 ad(5,28)*K_int/cost5 ad(5,29)*K_int/cost5 ad(5,30)*K_int/cost5 ad(5,31)*K_int/cost5 ad(5,32)*K_int/cost5 ad(5,33)*K_int/cost5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(6,1)*K_int/cost6 ad(6,2)*K_int/cost6 ad(6,3)*K_int/cost6 ad(6,4)*K_int/cost6 ad(6,5)*K_int/cost6 -ueq(1,6)/volume6-K_est*(somma6-pareti_esterne6(1,1))/cost6-K_int*(sum(area6))/cost6 ad(6,7)*K_int/cost6 ad(6,8)*K_int/cost6 ad(6,9)*K_int/cost6 ad(6,10)*K_int/cost6 ad(6,11)*K_int/cost6 ad(6,12)*K_int/cost6 ad(6,13)*K_int/cost6 ad(6,14)*K_int/cost6 ad(6,15)*K_int/cost6 ad(6,16)*K_int/cost6 ad(6,17)*K_int/cost6 ad(6,18)*K_int/cost6 ad(6,19)*K_int/cost6 ad(6,20)*K_int/cost6 ad(6,21)*K_int/cost6 ad(6,22)*K_int/cost6 ad(6,23)*K_int/cost6 ad(6,24)*K_int/cost6 ad(6,25)*K_int/cost6 ad(6,26)*K_int/cost6 ad(6,27)*K_int/cost6 ad(6,28)*K_int/cost6 ad(6,29)*K_int/cost6 ad(6,30)*K_int/cost6 ad(6,31)*K_int/cost6 ad(6,32)*K_int/cost6 ad(6,33)*K_int/cost6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(7,1)*K_int/cost7 ad(7,2)*K_int/cost7 ad(7,3)*K_int/cost7 ad(7,4)*K_int/cost7 ad(7,5)*K_int/cost7 ad(7,6)*K_int/cost7 -ueq(1,7)/volume7-K_est*(somma7-pareti_esterne7(1,1))/cost7-K_int*(sum(area7))/cost7 ad(7,8)*K_int/cost7 ad(7,9)*K_int/cost7 ad(7,10)*K_int/cost7 ad(7,11)*K_int/cost7 ad(7,12)*K_int/cost7 ad(7,13)*K_int/cost7 ad(7,14)*K_int/cost7 ad(7,15)*K_int/cost7 ad(7,16)*K_int/cost7 ad(7,17)*K_int/cost7 ad(7,18)*K_int/cost7 ad(7,19)*K_int/cost7 ad(7,20)*K_int/cost7 ad(7,21)*K_int/cost7 ad(7,22)*K_int/cost7 ad(7,23)*K_int/cost7 ad(7,24)*K_int/cost7 ad(7,25)*K_int/cost7 ad(7,26)*K_int/cost7 ad(7,27)*K_int/cost7 ad(7,28)*K_int/cost7 ad(7,29)*K_int/cost7 ad(7,30)*K_int/cost7 ad(7,31)*K_int/cost7 ad(7,32)*K_int/cost7 ad(7,33)*K_int/cost7 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(8,1)*K_int/cost8 ad(8,2)*K_int/cost8 ad(8,3)*K_int/cost8 ad(8,4)*K_int/cost8 ad(8,5)*K_int/cost8 ad(8,6)*K_int/cost8 ad(8,7)*K_int/cost8 -ueq(1,8)/volume8-K_est*(somma8-pareti_esterne8(1,1))/cost8-K_int*(sum(area8))/cost8 ad(8,9)*K_int/cost8 ad(8,10)*K_int/cost8 ad(8,11)*K_int/cost8 ad(8,12)*K_int/cost8 ad(8,13)*K_int/cost8 ad(8,14)*K_int/cost8 ad(8,15)*K_int/cost8 ad(8,16)*K_int/cost8 ad(8,17)*K_int/cost8 ad(8,18)*K_int/cost8 ad(8,19)*K_int/cost8 ad(8,20)*K_int/cost8 ad(8,21)*K_int/cost8 ad(8,22)*K_int/cost8 ad(8,23)*K_int/cost8 ad(8,24)*K_int/cost8 ad(8,25)*K_int/cost8 ad(8,26)*K_int/cost8 ad(8,27)*K_int/cost8 ad(8,28)*K_int/cost8 ad(8,29)*K_int/cost8 ad(8,30)*K_int/cost8 ad(8,31)*K_int/cost8 ad(8,32)*K_int/cost8 ad(8,33)*K_int/cost8 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(9,1)*K_int/cost9 ad(9,2)*K_int/cost9 ad(9,3)*K_int/cost9 ad(9,4)*K_int/cost9 ad(9,5)*K_int/cost9 ad(9,6)*K_int/cost9 ad(9,7)*K_int/cost9 ad(9,8)*K_int/cost9 -ueq(1,9)/volume9-K_est*(somma9-pareti_esterne9(1,1))/cost9-K_int*(sum(area9))/cost9 ad(9,10)*K_int/cost9 ad(9,11)*K_int/cost9 ad(9,12)*K_int/cost9 ad(9,13)*K_int/cost9 ad(9,14)*K_int/cost9 ad(9,15)*K_int/cost9 ad(9,16)*K_int/cost9 ad(9,17)*K_int/cost9 ad(9,18)*K_int/cost9 ad(9,19)*K_int/cost9 ad(9,20)*K_int/cost9 ad(9,21)*K_int/cost9 ad(9,22)*K_int/cost9 ad(9,23)*K_int/cost9 ad(9,24)*K_int/cost9 ad(9,25)*K_int/cost9 ad(9,26)*K_int/cost9 ad(9,27)*K_int/cost9 ad(9,28)*K_int/cost9 ad(9,29)*K_int/cost9 ad(9,30)*K_int/cost9 ad(9,31)*K_int/cost9 ad(9,32)*K_int/cost9 ad(9,33)*K_int/cost9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(10,1)*K_int/cost10 ad(10,2)*K_int/cost10 ad(10,3)*K_int/cost10 ad(10,4)*K_int/cost10 ad(10,5)*K_int/cost10 ad(10,6)*K_int/cost10 ad(10,7)*K_int/cost10 ad(10,8)*K_int/cost10 ad(10,9)*K_int/cost10 -ueq(1,10)/volume10-K_est*(somma10-pareti_esterne10(1,1))/cost10-K_int*(sum(area10))/cost10 ad(10,11)*K_int/cost10 ad(10,12)*K_int/cost10 ad(10,13)*K_int/cost10 ad(10,14)*K_int/cost10 ad(10,15)*K_int/cost10 ad(10,16)*K_int/cost10 ad(10,17)*K_int/cost10 ad(10,18)*K_int/cost10 ad(10,19)*K_int/cost10 ad(10,20)*K_int/cost10 ad(10,21)*K_int/cost10 ad(10,22)*K_int/cost10 ad(10,23)*K_int/cost10 ad(10,24)*K_int/cost10 ad(10,25)*K_int/cost10 ad(10,26)*K_int/cost10 ad(10,27)*K_int/cost10 ad(10,28)*K_int/cost10 ad(10,29)*K_int/cost10 ad(10,30)*K_int/cost10 ad(10,31)*K_int/cost10 ad(10,32)*K_int/cost10 ad(10,33)*K_int/cost10 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(11,1)*K_int/cost11 ad(11,2)*K_int/cost11 ad(11,3)*K_int/cost11 ad(11,4)*K_int/cost11 ad(11,5)*K_int/cost11 ad(11,6)*K_int/cost11 ad(11,7)*K_int/cost11 ad(11,8)*K_int/cost11 ad(11,9)*K_int/cost11 ad(11,10)*K_int/cost11 -ueq(1,11)/volume11-K_est*(somma11-pareti_esterne11(1,1))/cost11-K_int*(sum(area11))/cost11 ad(11,12)*K_int/cost11 ad(11,13)*K_int/cost11 ad(11,14)*K_int/cost11 ad(11,15)*K_int/cost11 ad(11,16)*K_int/cost11 ad(11,17)*K_int/cost11 ad(11,18)*K_int/cost11 ad(11,19)*K_int/cost11 ad(11,20)*K_int/cost11 ad(11,21)*K_int/cost11 ad(11,22)*K_int/cost11 ad(11,23)*K_int/cost11 ad(11,24)*K_int/cost11 ad(11,25)*K_int/cost11 ad(11,26)*K_int/cost11 ad(11,27)*K_int/cost11 ad(11,28)*K_int/cost11 ad(11,29)*K_int/cost11 ad(11,30)*K_int/cost11 ad(11,31)*K_int/cost11 ad(11,32)*K_int/cost11 ad(11,33)*K_int/cost11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(12,1)*K_int/cost12 ad(12,2)*K_int/cost12 ad(12,3)*K_int/cost12 ad(12,4)*K_int/cost12 ad(12,5)*K_int/cost12 ad(12,6)*K_int/cost12 ad(12,7)*K_int/cost12 ad(12,8)*K_int/cost12 ad(12,9)*K_int/cost12 ad(12,10)*K_int/cost12 ad(12,11)*K_int/cost12 -ueq(1,12)/volume12-K_est*(somma12-pareti_esterne12(1,1))/cost12-K_int*(sum(area12))/cost12 ad(12,13)*K_int/cost12 ad(12,14)*K_int/cost12 ad(12,15)*K_int/cost12 ad(12,16)*K_int/cost12 ad(12,17)*K_int/cost12 ad(12,18)*K_int/cost12 ad(12,19)*K_int/cost12 ad(12,20)*K_int/cost12 ad(12,21)*K_int/cost12 ad(12,22)*K_int/cost12 ad(12,23)*K_int/cost12 ad(12,24)*K_int/cost12 ad(12,25)*K_int/cost12 ad(12,26)*K_int/cost12 ad(12,27)*K_int/cost12 ad(12,28)*K_int/cost12 ad(12,29)*K_int/cost12 ad(12,30)*K_int/cost12 ad(12,31)*K_int/cost12 ad(12,32)*K_int/cost12 ad(12,33)*K_int/cost12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(13,1)*K_int/cost13 ad(13,2)*K_int/cost13 ad(13,3)*K_int/cost13 ad(13,4)*K_int/cost13 ad(13,5)*K_int/cost13 ad(13,6)*K_int/cost13 ad(13,7)*K_int/cost13 ad(13,8)*K_int/cost13 ad(13,9)*K_int/cost13 ad(13,10)*K_int/cost13 ad(13,11)*K_int/cost13 ad(13,12)*K_int/cost13 -ueq(1,13)/volume13-K_est*(somma13-pareti_esterne13(1,1))/cost13-K_int*(sum(area13))/cost13 ad(13,14)*K_int/cost13 ad(13,15)*K_int/cost13 ad(13,16)*K_int/cost13 ad(13,17)*K_int/cost13 ad(13,18)*K_int/cost13 ad(13,19)*K_int/cost13 ad(13,20)*K_int/cost13 ad(13,21)*K_int/cost13 ad(13,22)*K_int/cost13 ad(13,23)*K_int/cost13 ad(13,24)*K_int/cost13 ad(13,25)*K_int/cost13 ad(13,26)*K_int/cost13 ad(13,27)*K_int/cost13 ad(13,28)*K_int/cost13 ad(13,29)*K_int/cost13 ad(13,30)*K_int/cost13 ad(13,31)*K_int/cost13 ad(13,32)*K_int/cost13 ad(13,33)*K_int/cost13 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(14,1)*K_int/cost14 ad(14,2)*K_int/cost14 ad(14,3)*K_int/cost14 ad(14,4)*K_int/cost14 ad(14,5)*K_int/cost14 ad(14,6)*K_int/cost14 ad(14,7)*K_int/cost14 ad(14,8)*K_int/cost14 ad(14,9)*K_int/cost14 ad(14,10)*K_int/cost14 ad(14,11)*K_int/cost14 ad(14,12)*K_int/cost14 ad(14,13)*K_int/cost14 -ueq(1,14)/volume14-K_est*(somma14-pareti_esterne14(1,1))/cost14-K_int*(sum(area14))/cost14 ad(14,15)*K_int/cost14 ad(14,16)*K_int/cost14 ad(14,17)*K_int/cost14 ad(14,18)*K_int/cost14 ad(14,19)*K_int/cost14 ad(14,20)*K_int/cost14 ad(14,21)*K_int/cost14 ad(14,22)*K_int/cost14 ad(14,23)*K_int/cost14 ad(14,24)*K_int/cost14 ad(14,25)*K_int/cost14 ad(14,26)*K_int/cost14 ad(14,27)*K_int/cost14 ad(14,28)*K_int/cost14 ad(14,29)*K_int/cost14 ad(14,30)*K_int/cost14 ad(14,31)*K_int/cost14 ad(14,32)*K_int/cost14 ad(14,33)*K_int/cost14 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(15,1)*K_int/cost15 ad(15,2)*K_int/cost15 ad(15,3)*K_int/cost15 ad(15,4)*K_int/cost15 ad(15,5)*K_int/cost15 ad(15,6)*K_int/cost15 ad(15,7)*K_int/cost15 ad(15,8)*K_int/cost15 ad(15,9)*K_int/cost15 ad(15,10)*K_int/cost15 ad(15,11)*K_int/cost15 ad(15,12)*K_int/cost15 ad(15,13)*K_int/cost15 ad(15,14)*K_int/cost15 -ueq(1,15)/volume15-K_est*(somma15-pareti_esterne15(1,1))/cost15-K_int*(sum(area15))/cost15 ad(15,16)*K_int/cost15 ad(15,17)*K_int/cost15 ad(15,18)*K_int/cost15 ad(15,19)*K_int/cost15 ad(15,20)*K_int/cost15 ad(15,21)*K_int/cost15 ad(15,22)*K_int/cost15 ad(15,23)*K_int/cost15 ad(15,24)*K_int/cost15 ad(15,25)*K_int/cost15 ad(15,26)*K_int/cost15 ad(15,27)*K_int/cost15 ad(15,28)*K_int/cost15 ad(15,29)*K_int/cost15 ad(15,30)*K_int/cost15 ad(15,31)*K_int/cost15 ad(15,32)*K_int/cost15 ad(15,33)*K_int/cost15 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(16,1)*K_int/cost16 ad(16,2)*K_int/cost16 ad(16,3)*K_int/cost16 ad(16,4)*K_int/cost16 ad(16,5)*K_int/cost16 ad(16,6)*K_int/cost16 ad(16,7)*K_int/cost16 ad(16,8)*K_int/cost16 ad(16,9)*K_int/cost16 ad(16,10)*K_int/cost16 ad(16,11)*K_int/cost16 ad(16,12)*K_int/cost16 ad(16,13)*K_int/cost16 ad(16,14)*K_int/cost16 ad(16,15)*K_int/cost16 -ueq(1,16)/volume16-K_est*(somma16-pareti_esterne16(1,1))/cost16-K_int*(sum(area16))/cost16 ad(16,17)*K_int/cost16 ad(16,18)*K_int/cost16 ad(16,19)*K_int/cost16 ad(16,20)*K_int/cost16 ad(16,21)*K_int/cost16 ad(16,22)*K_int/cost16 ad(16,23)*K_int/cost16 ad(16,24)*K_int/cost16 ad(16,25)*K_int/cost16 ad(16,26)*K_int/cost16 ad(16,27)*K_int/cost16 ad(16,28)*K_int/cost16 ad(16,29)*K_int/cost16 ad(16,30)*K_int/cost16 ad(16,31)*K_int/cost16 ad(16,32)*K_int/cost16 ad(16,33)*K_int/cost16 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(17,1)*K_int/cost17 ad(17,2)*K_int/cost17 ad(17,3)*K_int/cost17 ad(17,4)*K_int/cost17 ad(17,5)*K_int/cost17 ad(17,6)*K_int/cost17 ad(17,7)*K_int/cost17 ad(17,8)*K_int/cost17 ad(17,9)*K_int/cost17 ad(17,10)*K_int/cost17 ad(17,11)*K_int/cost17 ad(17,12)*K_int/cost17 ad(17,13)*K_int/cost17 ad(17,14)*K_int/cost17 ad(17,15)*K_int/cost17 ad(17,16)*K_int/cost17 -ueq(1,17)/volume17-K_est*(somma17-pareti_esterne17(1,1))/cost17-K_int*(sum(area17))/cost17 ad(17,18)*K_int/cost17 ad(17,19)*K_int/cost17 ad(17,20)*K_int/cost17 ad(17,21)*K_int/cost17 ad(17,22)*K_int/cost17 ad(17,23)*K_int/cost17 ad(17,24)*K_int/cost17 ad(17,25)*K_int/cost17 ad(17,26)*K_int/cost17 ad(17,27)*K_int/cost17 ad(17,28)*K_int/cost17 ad(17,29)*K_int/cost17 ad(17,30)*K_int/cost17 ad(17,31)*K_int/cost17 ad(17,32)*K_int/cost17 ad(17,33)*K_int/cost17 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(18,1)*K_int/cost18 ad(18,2)*K_int/cost18 ad(18,3)*K_int/cost18 ad(18,4)*K_int/cost18 ad(18,5)*K_int/cost18 ad(18,6)*K_int/cost18 ad(18,7)*K_int/cost18 ad(18,8)*K_int/cost18 ad(18,9)*K_int/cost18 ad(18,10)*K_int/cost18 ad(18,11)*K_int/cost18 ad(18,12)*K_int/cost18 ad(18,13)*K_int/cost18 ad(18,14)*K_int/cost18 ad(18,15)*K_int/cost18 ad(18,16)*K_int/cost18 ad(18,17)*K_int/cost18 -ueq(1,18)/volume18-K_est*(somma18-pareti_esterne18(1,1))/cost18-K_int*(sum(area18))/cost18 ad(18,19)*K_int/cost18 ad(18,20)*K_int/cost18 ad(18,21)*K_int/cost18 ad(18,22)*K_int/cost18 ad(18,23)*K_int/cost18 ad(18,24)*K_int/cost18 ad(18,25)*K_int/cost18 ad(18,26)*K_int/cost18 ad(18,27)*K_int/cost18 ad(18,28)*K_int/cost18 ad(18,29)*K_int/cost18 ad(18,30)*K_int/cost18 ad(18,31)*K_int/cost18 ad(18,32)*K_int/cost18 ad(18,33)*K_int/cost18 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(19,1)*K_int/cost19 ad(19,2)*K_int/cost19 ad(19,3)*K_int/cost19 ad(19,4)*K_int/cost19 ad(19,5)*K_int/cost19 ad(19,6)*K_int/cost19 ad(19,7)*K_int/cost19 ad(19,8)*K_int/cost19 ad(19,9)*K_int/cost19 ad(19,10)*K_int/cost19 ad(19,11)*K_int/cost19 ad(19,12)*K_int/cost19 ad(19,13)*K_int/cost19 ad(19,14)*K_int/cost19 ad(19,15)*K_int/cost19 ad(19,16)*K_int/cost19 ad(19,17)*K_int/cost19 ad(19,18)*K_int/cost19 -ueq(1,19)/volume19-K_est*(somma19-pareti_esterne19(1,1))/cost19-K_int*(sum(area19))/cost19 ad(19,20)*K_int/cost19 ad(19,21)*K_int/cost19 ad(19,22)*K_int/cost19 ad(19,23)*K_int/cost19 ad(19,24)*K_int/cost19 ad(19,25)*K_int/cost19 ad(19,26)*K_int/cost19 ad(19,27)*K_int/cost19 ad(19,28)*K_int/cost19 ad(19,29)*K_int/cost19 ad(19,30)*K_int/cost19 ad(19,31)*K_int/cost19 ad(19,32)*K_int/cost19 ad(19,33)*K_int/cost19 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(20,1)*K_int/cost20 ad(20,2)*K_int/cost20 ad(20,3)*K_int/cost20 ad(20,4)*K_int/cost20 ad(20,5)*K_int/cost20 ad(20,6)*K_int/cost20 ad(20,7)*K_int/cost20 ad(20,8)*K_int/cost20 ad(20,9)*K_int/cost20 ad(20,10)*K_int/cost20 ad(20,11)*K_int/cost20 ad(20,12)*K_int/cost20 ad(20,13)*K_int/cost20 ad(20,14)*K_int/cost20 ad(20,15)*K_int/cost20 ad(20,16)*K_int/cost20 ad(20,17)*K_int/cost20 ad(20,18)*K_int/cost20 ad(20,19)*K_int/cost20 -ueq(1,20)/volume20-K_est*(somma20-pareti_esterne20(1,1))/cost20-K_int*(sum(area20))/cost20 ad(20,21)*K_int/cost20 ad(20,22)*K_int/cost20 ad(20,23)*K_int/cost20 ad(20,24)*K_int/cost20 ad(20,25)*K_int/cost20 ad(20,26)*K_int/cost20 ad(20,27)*K_int/cost20 ad(20,28)*K_int/cost20 ad(20,29)*K_int/cost20 ad(20,30)*K_int/cost20 ad(20,31)*K_int/cost20 ad(20,32)*K_int/cost20 ad(20,33)*K_int/cost20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(21,1)*K_int/cost21 ad(21,2)*K_int/cost21 ad(21,3)*K_int/cost21 ad(21,4)*K_int/cost21 ad(21,5)*K_int/cost21 ad(21,6)*K_int/cost21 ad(21,7)*K_int/cost21 ad(21,8)*K_int/cost21 ad(21,9)*K_int/cost21 ad(21,10)*K_int/cost21 ad(21,11)*K_int/cost21 ad(21,12)*K_int/cost21 ad(21,13)*K_int/cost21 ad(21,14)*K_int/cost21 ad(21,15)*K_int/cost21 ad(21,16)*K_int/cost21 ad(21,17)*K_int/cost21 ad(21,18)*K_int/cost21 ad(21,19)*K_int/cost21 ad(21,20)*K_int/cost21 -ueq(1,21)/volume21-K_est*(somma21-pareti_esterne21(1,1))/cost21-K_int*(sum(area21))/cost21 ad(21,22)*K_int/cost21 ad(21,23)*K_int/cost21 ad(21,24)*K_int/cost21 ad(21,25)*K_int/cost21 ad(21,26)*K_int/cost21 ad(21,27)*K_int/cost21 ad(21,28)*K_int/cost21 ad(21,29)*K_int/cost21 ad(21,30)*K_int/cost21 ad(21,31)*K_int/cost21 ad(21,32)*K_int/cost21 ad(21,33)*K_int/cost21 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(22,1)*K_int/cost22 ad(22,2)*K_int/cost22 ad(22,3)*K_int/cost22 ad(22,4)*K_int/cost22 ad(22,5)*K_int/cost22 ad(22,6)*K_int/cost22 ad(22,7)*K_int/cost22 ad(22,8)*K_int/cost22 ad(22,9)*K_int/cost22 ad(22,10)*K_int/cost22 ad(22,11)*K_int/cost22 ad(22,12)*K_int/cost22 ad(22,13)*K_int/cost22 ad(22,14)*K_int/cost22 ad(22,15)*K_int/cost22 ad(22,16)*K_int/cost22 ad(22,17)*K_int/cost22 ad(22,18)*K_int/cost22 ad(22,19)*K_int/cost22 ad(22,20)*K_int/cost22 ad(22,21)*K_int/cost22 -ueq(1,22)/volume22-K_est*(somma22-pareti_esterne22(1,1))/cost22-K_int*(sum(area22))/cost22 ad(22,23)*K_int/cost22 ad(22,24)*K_int/cost22 ad(22,25)*K_int/cost22 ad(22,26)*K_int/cost22 ad(22,27)*K_int/cost22 ad(22,28)*K_int/cost22 ad(22,29)*K_int/cost22 ad(22,30)*K_int/cost22 ad(22,31)*K_int/cost22 ad(22,32)*K_int/cost22 ad(22,33)*K_int/cost22 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(23,1)*K_int/cost23 ad(23,2)*K_int/cost23 ad(23,3)*K_int/cost23 ad(23,4)*K_int/cost23 ad(23,5)*K_int/cost23 ad(23,6)*K_int/cost23 ad(23,7)*K_int/cost23 ad(23,8)*K_int/cost23 ad(23,9)*K_int/cost23 ad(23,10)*K_int/cost23 ad(23,11)*K_int/cost23 ad(23,12)*K_int/cost23 ad(23,13)*K_int/cost23 ad(23,14)*K_int/cost23 ad(23,15)*K_int/cost23 ad(23,16)*K_int/cost23 ad(23,17)*K_int/cost23 ad(23,18)*K_int/cost23 ad(23,19)*K_int/cost23 ad(23,20)*K_int/cost23 ad(23,21)*K_int/cost23 ad(23,22)*K_int/cost23 -ueq(1,23)/volume23-K_est*(somma23-pareti_esterne23(1,1))/cost23-K_int*(sum(area23))/cost23 ad(23,24)*K_int/cost23 ad(23,25)*K_int/cost23 ad(23,26)*K_int/cost23 ad(23,27)*K_int/cost23 ad(23,28)*K_int/cost23 ad(23,29)*K_int/cost23 ad(23,30)*K_int/cost23 ad(23,31)*K_int/cost23 ad(23,32)*K_int/cost23 ad(23,33)*K_int/cost23 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(24,1)*K_int/cost24 ad(24,2)*K_int/cost24 ad(24,3)*K_int/cost24 ad(24,4)*K_int/cost24 ad(24,5)*K_int/cost24 ad(24,6)*K_int/cost24 ad(24,7)*K_int/cost24 ad(24,8)*K_int/cost24 ad(24,9)*K_int/cost24 ad(24,10)*K_int/cost24 ad(24,11)*K_int/cost24 ad(24,12)*K_int/cost24 ad(24,13)*K_int/cost24 ad(24,14)*K_int/cost24 ad(24,15)*K_int/cost24 ad(24,16)*K_int/cost24 ad(24,17)*K_int/cost24 ad(24,18)*K_int/cost24 ad(24,19)*K_int/cost24 ad(24,20)*K_int/cost24 ad(24,21)*K_int/cost24 ad(24,22)*K_int/cost24 ad(24,23)*K_int/cost24 -ueq(1,24)/volume24-K_est*(somma24-pareti_esterne24(1,1))/cost24-K_int*(sum(area24))/cost24 ad(24,25)*K_int/cost24 ad(24,26)*K_int/cost24 ad(24,27)*K_int/cost24 ad(24,28)*K_int/cost24 ad(24,29)*K_int/cost24 ad(24,30)*K_int/cost24 ad(24,31)*K_int/cost24 ad(24,32)*K_int/cost24 ad(24,33)*K_int/cost24 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(25,1)*K_int/cost25 ad(25,2)*K_int/cost25 ad(25,3)*K_int/cost25 ad(25,4)*K_int/cost25 ad(25,5)*K_int/cost25 ad(25,6)*K_int/cost25 ad(25,7)*K_int/cost25 ad(25,8)*K_int/cost25 ad(25,9)*K_int/cost25 ad(25,10)*K_int/cost25 ad(25,11)*K_int/cost25 ad(25,12)*K_int/cost25 ad(25,13)*K_int/cost25 ad(25,14)*K_int/cost25 ad(25,15)*K_int/cost25 ad(25,16)*K_int/cost25 ad(25,17)*K_int/cost25 ad(25,18)*K_int/cost25 ad(25,19)*K_int/cost25 ad(25,20)*K_int/cost25 ad(25,21)*K_int/cost25 ad(25,22)*K_int/cost25 ad(25,23)*K_int/cost25 ad(25,24)*K_int/cost25 -ueq(1,25)/volume25-K_est*(somma25-pareti_esterne25(1,1))/cost25-K_int*(sum(area25))/cost25 ad(25,26)*K_int/cost25 ad(25,27)*K_int/cost25 ad(25,28)*K_int/cost25 ad(25,29)*K_int/cost25 ad(25,30)*K_int/cost25 ad(25,31)*K_int/cost25 ad(25,32)*K_int/cost25 ad(25,33)*K_int/cost25 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(26,1)*K_int/cost26 ad(26,2)*K_int/cost26 ad(26,3)*K_int/cost26 ad(26,4)*K_int/cost26 ad(26,5)*K_int/cost26 ad(26,6)*K_int/cost26 ad(26,7)*K_int/cost26 ad(26,8)*K_int/cost26 ad(26,9)*K_int/cost26 ad(26,10)*K_int/cost26 ad(26,11)*K_int/cost26 ad(26,12)*K_int/cost26 ad(26,13)*K_int/cost26 ad(26,14)*K_int/cost26 ad(26,15)*K_int/cost26 ad(26,16)*K_int/cost26 ad(26,17)*K_int/cost26 ad(26,18)*K_int/cost26 ad(26,19)*K_int/cost26 ad(26,20)*K_int/cost26 ad(26,21)*K_int/cost26 ad(26,22)*K_int/cost26 ad(26,23)*K_int/cost26 ad(26,24)*K_int/cost26 ad(26,25)*K_int/cost26 -ueq(1,26)/volume26-K_est*(somma26-pareti_esterne26(1,1))/cost26-K_int*(sum(area26))/cost26 ad(26,27)*K_int/cost26 ad(26,28)*K_int/cost26 ad(26,29)*K_int/cost26 ad(26,30)*K_int/cost26 ad(26,31)*K_int/cost26 ad(26,32)*K_int/cost26 ad(26,33)*K_int/cost26 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(27,1)*K_int/cost27 ad(27,2)*K_int/cost27 ad(27,3)*K_int/cost27 ad(27,4)*K_int/cost27 ad(27,5)*K_int/cost27 ad(27,6)*K_int/cost27 ad(27,7)*K_int/cost27 ad(27,8)*K_int/cost27 ad(27,9)*K_int/cost27 ad(27,10)*K_int/cost27 ad(27,11)*K_int/cost27 ad(27,12)*K_int/cost27 ad(27,13)*K_int/cost27 ad(27,14)*K_int/cost27 ad(27,15)*K_int/cost27 ad(27,16)*K_int/cost27 ad(27,17)*K_int/cost27 ad(27,18)*K_int/cost27 ad(27,19)*K_int/cost27 ad(27,20)*K_int/cost27 ad(27,21)*K_int/cost27 ad(27,22)*K_int/cost27 ad(27,23)*K_int/cost27 ad(27,24)*K_int/cost27 ad(27,25)*K_int/cost27 ad(27,26)*K_int/cost27 -ueq(1,27)/volume27-K_est*(somma27-pareti_esterne27(1,1))/cost27-K_int*(sum(area27))/cost27 ad(27,28)*K_int/cost27 ad(27,29)*K_int/cost27 ad(27,30)*K_int/cost27 ad(27,31)*K_int/cost27 ad(27,32)*K_int/cost27 ad(27,33)*K_int/cost27 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(28,1)*K_int/cost28 ad(28,2)*K_int/cost28 ad(28,3)*K_int/cost28 ad(28,4)*K_int/cost28 ad(28,5)*K_int/cost28 ad(28,6)*K_int/cost28 ad(28,7)*K_int/cost28 ad(28,8)*K_int/cost28 ad(28,9)*K_int/cost28 ad(28,10)*K_int/cost28 ad(28,11)*K_int/cost28 ad(28,12)*K_int/cost28 ad(28,13)*K_int/cost28 ad(28,14)*K_int/cost28 ad(28,15)*K_int/cost28 ad(28,16)*K_int/cost28 ad(28,17)*K_int/cost28 ad(28,18)*K_int/cost28 ad(28,19)*K_int/cost28 ad(28,20)*K_int/cost28 ad(28,21)*K_int/cost28 ad(28,22)*K_int/cost28 ad(28,23)*K_int/cost28 ad(28,24)*K_int/cost28 ad(28,25)*K_int/cost28 ad(28,26)*K_int/cost28 ad(28,27)*K_int/cost28 -ueq(1,28)/volume28-K_est*(somma28-pareti_esterne28(1,1))/cost28-K_int*(sum(area28))/cost28 ad(28,29)*K_int/cost28 ad(28,30)*K_int/cost28 ad(28,31)*K_int/cost28 ad(28,32)*K_int/cost28 ad(28,33)*K_int/cost28 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(29,1)*K_int/cost29 ad(29,2)*K_int/cost29 ad(29,3)*K_int/cost29 ad(29,4)*K_int/cost29 ad(29,5)*K_int/cost29 ad(29,6)*K_int/cost29 ad(29,7)*K_int/cost29 ad(29,8)*K_int/cost29 ad(29,9)*K_int/cost29 ad(29,10)*K_int/cost29 ad(29,11)*K_int/cost29 ad(29,12)*K_int/cost29 ad(29,13)*K_int/cost29 ad(29,14)*K_int/cost29 ad(29,15)*K_int/cost29 ad(29,16)*K_int/cost29 ad(29,17)*K_int/cost29 ad(29,18)*K_int/cost29 ad(29,19)*K_int/cost29 ad(29,20)*K_int/cost29 ad(29,21)*K_int/cost29 ad(29,22)*K_int/cost29 ad(29,23)*K_int/cost29 ad(29,24)*K_int/cost29 ad(29,25)*K_int/cost29 ad(29,26)*K_int/cost29 ad(29,27)*K_int/cost29 ad(29,28)*K_int/cost29 -ueq(1,29)/volume29-K_est*(somma29-pareti_esterne29(1,1))/cost29-K_int*(sum(area29))/cost29 ad(29,30)*K_int/cost29 ad(29,31)*K_int/cost29 ad(29,32)*K_int/cost29 ad(29,33)*K_int/cost29 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(30,1)*K_int/cost30 ad(30,2)*K_int/cost30 ad(30,3)*K_int/cost30 ad(30,4)*K_int/cost30 ad(30,5)*K_int/cost30 ad(30,6)*K_int/cost30 ad(30,7)*K_int/cost30 ad(30,8)*K_int/cost30 ad(30,9)*K_int/cost30 ad(30,10)*K_int/cost30 ad(30,11)*K_int/cost30 ad(30,12)*K_int/cost30 ad(30,13)*K_int/cost30 ad(30,14)*K_int/cost30 ad(30,15)*K_int/cost30 ad(30,16)*K_int/cost30 ad(30,17)*K_int/cost30 ad(30,18)*K_int/cost30 ad(30,19)*K_int/cost30 ad(30,20)*K_int/cost30 ad(30,21)*K_int/cost30 ad(30,22)*K_int/cost30 ad(30,23)*K_int/cost30 ad(30,24)*K_int/cost30 ad(30,25)*K_int/cost30 ad(30,26)*K_int/cost30 ad(30,27)*K_int/cost30 ad(30,28)*K_int/cost30 ad(30,29)*K_int/cost30 -ueq(1,30)/volume30-K_est*(somma30-pareti_esterne30(1,1))/cost30-K_int*(sum(area30))/cost30 ad(30,31)*K_int/cost30 ad(30,32)*K_int/cost30 ad(30,33)*K_int/cost30 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(31,1)*K_int/cost31 ad(31,2)*K_int/cost31 ad(31,3)*K_int/cost31 ad(31,4)*K_int/cost31 ad(31,5)*K_int/cost31 ad(31,6)*K_int/cost31 ad(31,7)*K_int/cost31 ad(31,8)*K_int/cost31 ad(31,9)*K_int/cost31 ad(31,10)*K_int/cost31 ad(31,11)*K_int/cost31 ad(31,12)*K_int/cost31 ad(31,13)*K_int/cost31 ad(31,14)*K_int/cost31 ad(31,15)*K_int/cost31 ad(31,16)*K_int/cost31 ad(31,17)*K_int/cost31 ad(31,18)*K_int/cost31 ad(31,19)*K_int/cost31 ad(31,20)*K_int/cost31 ad(31,21)*K_int/cost31 ad(31,22)*K_int/cost31 ad(31,23)*K_int/cost31 ad(31,24)*K_int/cost31 ad(31,25)*K_int/cost31 ad(31,26)*K_int/cost31 ad(31,27)*K_int/cost31 ad(31,28)*K_int/cost31 ad(31,29)*K_int/cost31 ad(31,30)*K_int/cost31 -ueq(1,31)/volume31-K_est*(somma31-pareti_esterne31(1,1))/cost31-K_int*(sum(area31))/cost31 ad(31,32)*K_int/cost31 ad(31,33)*K_int/cost31 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(32,1)*K_int/cost32 ad(32,2)*K_int/cost32 ad(32,3)*K_int/cost32 ad(32,4)*K_int/cost32 ad(32,5)*K_int/cost32 ad(32,6)*K_int/cost32 ad(32,7)*K_int/cost32 ad(32,8)*K_int/cost32 ad(32,9)*K_int/cost32 ad(32,10)*K_int/cost32 ad(32,11)*K_int/cost32 ad(32,12)*K_int/cost32 ad(32,13)*K_int/cost32 ad(32,14)*K_int/cost32 ad(32,15)*K_int/cost32 ad(32,16)*K_int/cost32 ad(32,17)*K_int/cost32 ad(32,18)*K_int/cost32 ad(32,19)*K_int/cost32 ad(32,20)*K_int/cost32 ad(32,21)*K_int/cost32 ad(32,22)*K_int/cost32 ad(32,23)*K_int/cost32 ad(32,24)*K_int/cost32 ad(32,25)*K_int/cost32 ad(32,26)*K_int/cost32 ad(32,27)*K_int/cost32 ad(32,28)*K_int/cost32 ad(32,29)*K_int/cost32 ad(32,30)*K_int/cost32 ad(32,31)*K_int/cost32 -ueq(1,32)/volume32-K_est*(somma32-pareti_esterne32(1,1))/cost32-K_int*(sum(area32))/cost32 ad(32,33)*K_int/cost32 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    ad(33,1)*K_int/cost33 ad(33,2)*K_int/cost33 ad(33,3)*K_int/cost33 ad(33,4)*K_int/cost33 ad(33,5)*K_int/cost33 ad(33,6)*K_int/cost33 ad(33,7)*K_int/cost33 ad(33,8)*K_int/cost33 ad(33,9)*K_int/cost33 ad(33,10)*K_int/cost33 ad(33,11)*K_int/cost33 ad(33,12)*K_int/cost33 ad(33,13)*K_int/cost33 ad(33,14)*K_int/cost33 ad(33,15)*K_int/cost33 ad(33,16)*K_int/cost33 ad(33,17)*K_int/cost33 ad(33,18)*K_int/cost33 ad(33,19)*K_int/cost33 ad(33,20)*K_int/cost33 ad(33,21)*K_int/cost33 ad(33,22)*K_int/cost33 ad(33,23)*K_int/cost33 ad(33,24)*K_int/cost33 ad(33,25)*K_int/cost33 ad(33,26)*K_int/cost33 ad(33,27)*K_int/cost33 ad(33,28)*K_int/cost33 ad(33,29)*K_int/cost33 ad(33,30)*K_int/cost33 ad(33,31)*K_int/cost33 ad(33,32)*K_int/cost33 -ueq(1,33)/volume33-K_est*(somma33-pareti_esterne33(1,1))/cost33-K_int*(sum(area33))/cost33 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,1)*c_a/cost1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,2)*c_a/cost2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,3)*c_a/cost3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,4)*c_a/cost4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,5)*c_a/cost5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,6)*c_a/cost6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,7)*c_a/cost7 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,8)*c_a/cost8 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,9)*c_a/cost9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,10)*c_a/cost10 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,11)*c_a/cost11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,12)*c_a/cost12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,13)*c_a/cost13 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,14)*c_a/cost14 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,15)*c_a/cost15 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0  
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,16)*c_a/cost16 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,17)*c_a/cost17 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,18)*c_a/cost18 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,19)*c_a/cost19 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,20)*c_a/cost20 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,21)*c_a/cost21 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,22)*c_a/cost22 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,23)*c_a/cost23 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,24)*c_a/cost24 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,25)*c_a/cost25 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,26)*c_a/cost26 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,27)*c_a/cost27 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,28)*c_a/cost28 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,29)*c_a/cost29 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,30)*c_a/cost30 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,31)*c_a/cost31 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,32)*c_a/cost32 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -ueq(1,33)*c_a/cost33];
    
    
   
B = [-xeq(1,1)/volume1 -L/cost1 1/cost1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 -xeq(1,2)/volume2 -L/cost2 1/cost2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 -xeq(1,3)/volume3 -L/cost3 1/cost3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 -xeq(1,4)/volume4 -L/cost4 1/cost4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,5)/volume5 -L/cost5 1/cost5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,6)/volume6 -L/cost6 1/cost6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,7)/volume7 -L/cost7 1/cost7 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,8)/volume8 -L/cost8 1/cost8 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,9)/volume9 -L/cost9 1/cost9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,10)/volume10 -L/cost10 1/cost10 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,11)/volume11 -L/cost11 1/cost11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,12)/volume12 -L/cost12 1/cost12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,13)/volume13 -L/cost13 1/cost13 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,14)/volume14 -L/cost14 1/cost14 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,15)/volume15 -L/cost15 1/cost15 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,16)/volume16 -L/cost16 1/cost16 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,17)/volume17 -L/cost17 1/cost17 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,18)/volume18 -L/cost18 1/cost18 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,19)/volume19 -L/cost19 1/cost19 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,20)/volume20 -L/cost20 1/cost20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,21)/volume21 -L/cost21 1/cost21 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,22)/volume22 -L/cost22 1/cost22 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0  
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,23)/volume23 -L/cost23 1/cost23 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,24)/volume24 -L/cost24 1/cost24 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,25)/volume25 -L/cost25 1/cost25 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,26)/volume26 -L/cost26 1/cost26 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,27)/volume27 -L/cost27 1/cost27 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,28)/volume28 -L/cost28 1/cost28 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,29)/volume29 -L/cost29 1/cost29 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,30)/volume30 -L/cost30 1/cost30 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,31)/volume31 -L/cost31 1/cost31 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,32)/volume32 -L/cost32 1/cost32 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,33)/volume33 -L/cost33 1/cost33  
    -xeq(1,34)*c_a/cost1 c_a/cost1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 -xeq(1,35)*c_a/cost2 c_a/cost2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 -xeq(1,36)*c_a/cost3 c_a/cost3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 -xeq(1,37)*c_a/cost4 c_a/cost4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,38)*c_a/cost5 c_a/cost5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,39)*c_a/cost6 c_a/cost6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,40)*c_a/cost7 c_a/cost7 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,41)*c_a/cost8 c_a/cost8 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,42)*c_a/cost9 c_a/cost9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,43)*c_a/cost10 c_a/cost10 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0  
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,44)*c_a/cost11 c_a/cost11 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,45)*c_a/cost12 c_a/cost12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,46)*c_a/cost13 c_a/cost13 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,47)*c_a/cost14 c_a/cost14 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,48)*c_a/cost15 c_a/cost15 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,49)*c_a/cost16 c_a/cost16 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,50)*c_a/cost17 c_a/cost17 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,51)*c_a/cost18 c_a/cost18 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,52)*c_a/cost19 c_a/cost19 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,53)*c_a/cost20 c_a/cost20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,54)*c_a/cost21 c_a/cost21 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,55)*c_a/cost22 c_a/cost22 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,56)*c_a/cost23 c_a/cost23 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,57)*c_a/cost24 c_a/cost24 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,58)*c_a/cost25 c_a/cost25 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,59)*c_a/cost26 c_a/cost26 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,60)*c_a/cost27 c_a/cost27 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,61)*c_a/cost28 c_a/cost28 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,62)*c_a/cost29 c_a/cost29 0 0 0 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,63)*c_a/cost30 c_a/cost30 0 0 0 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,64)*c_a/cost31 c_a/cost31 0 0 0 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,65)*c_a/cost32 c_a/cost32 0 0 0 0
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -xeq(1,66)*c_a/cost33 c_a/cost33 0
    ];    
    
C= eye(66);
  
raggiungibilita= ctrb(A,B);
raggiungibilita=rank(raggiungibilita);
osservabilita = obsv (A,C);
osservabilita = rank(osservabilita);
eigenvalue = eig (A);

Q = eye (66);
%Q = 100*Q;
%Q = diag ([1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100 100]);
R = eye (99);
%R =R*10;
%R = diag([1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 1000 1 1 ]);
[G,K,e] = lqr(A,B,Q,R); % K soluzione dell'equazione di riccati, G è la matrice dei guadagni di kalman
