//STORIA1.OX
//
// 21 aprile 2002  confindustria#2
// prende dati dal 1862 al 1924
// 
// 
// 

#include<oxstd.h>
#include<oxdraw.h>
#include "c:\dpdox\dpd.h"
#include "c:\dpdox\dpd.ox"

main()
{
      decl dpd=new DPD();     decl i; decl nprov; 
      decl t = 63 ; // dal 1962 al 1923 piu' anno doppi0
      decl ncol=24; // num. colonne. Ordine: Anno, Pie-Vaa, Lom, Trento, Ven
                                          // FVG, Udine,    Lig, EmR,    Tos
                                          // Umb, Mar,      Laz  Abr,    Mol
                                          // Cam, Pug,      Bas, Cal,    Sic
                                          // Sar, SPESEP,   Tripoli, Zara
//    Dati Fenoaltea
//    Tabelle 1, 2  e 3 Fenoaltea (1985)
//    stato, stime finali (manut., op. nuov/edifici op.nuov/altre); comunali (=), provinciali (=)
      decl fenoalt    =  loadmat("c:/dpdox/storia/dati/fenoaltea.dat", 1);  //
      fenoalt=shape(fenoalt, 9, 53)'; 
      print("\n fenoaltea", fenoalt);

//    strade
      decl strad    =  loadmat("c:/gauss/storia/destefani-strade-62-24.dat", 1);  //
      strad=shape(strad, ncol, t)'; 
//      print("\n strad ", strad);  
//      print("\n sumcol strad ", ((sumc(strad))/1000)');  
//      print("\n sumrows strad ", ((sumr(strad)-strad[][6])/1000)');  
//    ora mette a posto il friuli ecc.
      strad=strad[][0:3]~(strad[][4]-strad[][6])~(strad[][5]+strad[][6])~strad[][7:23];
//      print("\n sumcol strad ", ((sumc(strad))/1000)');  
//      print("\n sumrows strad ", ((sumr(strad))/1000)');  
//    print("\n grandtot strad", sumc(sumr(strad[][1:ncol-2])/1000000)); 
//    ricostruisce i dati per gli anni dal 1884 in poi
      decl strad2=strad[0:21][]|(strad[22][]+strad[23][]/2)|(strad[23:t-2][]/2+strad[24:t-1][]/2)|strad[t-1][]; 
//      print("\n strad2", strad2/1000);  
// print("\n grandtot strad", sumc(sumr(strad[][1:ncol-2])/1000000)); 
// print("\n grandtot strad2", sumc(sumr(strad2[][1:ncol-2])/1000000)); 

//    idrauliche
      decl idrau    =  loadmat("c:/gauss/storia/destefani-idrauliche.dat", 1);  //
      idrau=shape(idrau, ncol, t)'; 
//      print("\n idrau ", idrau);  
//      print("\n sumcol idrau ",  ((sumc(idrau[][1:ncol-1]))/1000)');  
//      print("\n sumrows idrau ", ((sumr(idrau[][1:ncol-1])-idrau[][6])/1000)');  
//    ora mette a posto il friuli ecc.
      idrau=idrau[][0:3]~(idrau[][4]-idrau[][6])~(idrau[][5]+idrau[][6])~idrau[][7:23];
//      print("\n sumcol idrau ",  ((sumc(idrau[][1:ncol-2]))/1000)');  
//      print("\n sumrows idrau ", ((sumr(idrau[][1:ncol-2]))/1000)');  
//    ricostruisce i dati per gli anni dal 1884 in poi
      decl idrau2=idrau[0:21][]|(idrau[22][]+idrau[23][]/2)|(idrau[23:t-2][]/2+idrau[24:t-1][]/2)|idrau[t-1][];       
//      print("\n idrau ", idrau/1000);  

//    alluvioni, piene e frane
      decl alluv    =  loadmat("c:/gauss/storia/destefani-alluvioni.dat", 1);  //
      alluv=shape(alluv, ncol, t)'; 
//      print("\n alluv ", alluv);  
//      print("\n sumcol alluv ",  ((sumc(alluv[][1:ncol-1]))/1000)');  
//      print("\n sumrows alluv ", ((sumr(alluv[][1:ncol-1])-alluv[][6])/1000)');  
//    ora mette a posto il friuli ecc.
      alluv=alluv[][0:3]~(alluv[][4]-alluv[][6])~(alluv[][5]+alluv[][6])~alluv[][7:23];
//      print("\n sumcol alluv ",  ((sumc(alluv[][1:ncol-2]))/1000)');  
//      print("\n sumrows alluv ", ((sumr(alluv[][1:ncol-2]))/1000)');  
//    ricostruisce i dati per gli anni dal 1884 in poi
      decl alluv2=alluv[0:21][]|(alluv[22][]+alluv[23][]/2)|(alluv[23:t-2][]/2+alluv[24:t-1][]/2)|alluv[t-1][];       
//      print("\n alluv ", alluv/1000);  


//    opere di bonifiche
      decl bonif    =  loadmat("c:/gauss/storia/destefani-bonifiche.dat", 1);  //
      bonif= shape(bonif, ncol, t)'; 
//      print("\n bonif ", bonif);  
//      print("\n sumcol bonif ",  ((sumc(bonif[][1:ncol-1]))/1000)');  
//      print("\n sumrows bonif ", ((sumr(bonif[][1:ncol-1])-bonif[][6])/1000)');  
//    ora mette a posto il friuli ecc.
      bonif=bonif[][0:3]~(bonif[][4]-bonif[][6])~(bonif[][5]+bonif[][6])~bonif[][7:23];
//      print("\n sumcol bonif ",  ((sumc(bonif[][1:ncol-2]))/1000)');  
//      print("\n sumrows bonif ", ((sumr(bonif[][1:ncol-2]))/1000)');  
//    ricostruisce i dati per gli anni dal 1884 in poi
      decl bonif2=bonif[0:21][]|(bonif[22][]+bonif[23][]/2)|(bonif[23:t-2][]/2+bonif[24:t-1][]/2)|bonif[t-1][];       
//      print("\n bonif ", bonif/1000);  

//    opere Marittime
      decl marit    =  loadmat("c:/gauss/storia/destefani-marittime.dat", 1);  //
      marit= shape(marit, ncol, t)'; 
//      print("\n marit ", marit);  
//      print("\n sumcol marit ",  ((sumc(marit[][1:ncol-1]))/1000)');  
//      print("\n sumrows marit ", ((sumr(marit[][1:ncol-1])-marit[][6])/1000)');  
//    ora mette a posto il friuli ecc.
      marit=marit[][0:3]~(marit[][4]-marit[][6])~(marit[][5]+marit[][6])~marit[][7:23];
//      print("\n sumcol marit ",  ((sumc(marit[][1:ncol-2]))/1000)');  
//      print("\n sumrows marit ", ((sumr(marit[][1:ncol-2]))/1000)');  
//    ricostruisce i dati per gli anni dal 1884 in poi
      decl marit2=marit[0:21][]|(marit[22][]+marit[23][]/2)|(marit[23:t-2][]/2+marit[24:t-1][]/2)|marit[t-1][];       
//      print("\n marit ", marit/1000);  
    
//    Edifici
      decl edifi    =  loadmat("c:/gauss/storia/destefani-edifici.dat", 1);  //
      edifi= shape(edifi, ncol, t)'; 
//      print("\n edifi ", edifi);  
//      print("\n sumcol edifi ",  ((sumc(edifi[][1:ncol-1]))/1000)');  
//      print("\n sumrows edifi ", ((sumr(edifi[][1:ncol-1])-edifi[][6])/1000)');  
//    ora mette a posto il friuli ecc.
      edifi=edifi[][0:3]~(edifi[][4]-edifi[][6])~(edifi[][5]+edifi[][6])~edifi[][7:23];
//      print("\n sumcol edifi ",  ((sumc(edifi[][1:ncol-2]))/1000)');  
//      print("\n sumrows edifi ", ((sumr(edifi[][1:ncol-2]))/1000)');  
//    ricostruisce i dati per gli anni dal 1884 in poi
      decl edifi2=edifi[0:21][]|(edifi[22][]+edifi[23][]/2)|(edifi[23:t-2][]/2+edifi[24:t-1][]/2)|edifi[t-1][];       
//    print("\n edifi ", edifi/1000);  

// ferrovie:
//   It.Sett.  It.Centr. It.Mer-ins Sommeindivise       
//   1088173501 375390501 741020808 138748661
      ncol=6;
      decl ferro    =  loadmat("c:/gauss/storia/destefani-fer.dat", 1);  //
      ferro= shape(ferro, ncol, t)'; 
//      print("\n ferro ", ferro);  
//      print("\n sumcol ferro ",  ((sumc(ferro[][1:ncol-1]))/1000)');  
      print("\n sumrows ferro ", ((sumr(ferro[][1:ncol-1]))/1000));  

// tot spese stato fenoaltea (col 5-6 tab 1) + ferro
print("\n spesestatconfer ", sumr(fenoalt[][1:2])~(0|sumr(ferro[0:51][1:ncol-1]))/1000000);
//  Attribuzione ferrovie sino al 1897
    decl mat=<1088173501, 375390501, 741020808, 138748661>;
    decl ferrorip=ferro;
    ferrorip[0:36][2:5]=ferro[0:36][1]*(mat/sumr(mat));
//  print("\n ferrorip", ferrorip);
//    ricostruisce i dati per gli anni dal 1884 in poi
      decl ferro2=ferro[0:21][]|(ferro[22][]+ferro[23][]/2)|(ferro[23:t-2][]/2+ferro[24:t-1][]/2)|ferro[t-1][];       
//    anche per ferrorip ricostruisce i dati per gli anni dal 1884 in poi
      ferrorip=ferrorip[0:21][]|(ferrorip[22][]+ferrorip[23][]/2)|(ferrorip[23:t-2][]/2+ferrorip[24:t-1][]/2)|ferrorip[t-1][];       
//      print("\n ferrorip ", ferrorip/1000);  

// Non ripartite, LLPP:
//   Trasp. capitale, Acq. pugliese, Nuove provincie
      ncol=4;
      decl nripa    =  loadmat("c:/gauss/storia/destefani-nonripartite.dat", 1);  //
      nripa= shape(nripa, ncol, t)'; 
//      print("\n nripa ", nripa);  
//      print("\n sumcol nripa ",  ((sumc(nripa[][1:ncol-1]))/1000)');  
//      print("\n sumrows nripa ", ((sumr(nripa[][1:ncol-1]))/1000)');  
//    ricostruisce i dati per gli anni dal 1884 in poi
      decl nripa2=nripa[0:21][]|(nripa[22][]+nripa[23][]/2)|(nripa[23:t-2][]/2+nripa[24:t-1][]/2)|nripa[t-1][];       
//      print("\n nripa ", nripa/1000);  

// Altri Ministeri: Prospetto 10-F pg. 1160 e sg.
//   Limitatamente a: sede minint, minmar, mingiust, ministr, minecnaz, minfin, corteconti, 
//   tesoro edifici vari, finanze edifici vari, edesami, edpostccrr, capporto, biologia marina messina, sbartolomeo cagliari,
//   ed giudiziari, ed carcerari, monumento a vittorio emanuale, ec-naz ed vari, mincolonie, univ-ist superiori, stazioni sanitarie (21 in tutto)
      ncol=22;
      decl almin    =  loadmat("c:/gauss/storia/destefani-altriministeri.dat", 1);  //
      almin= shape(almin, ncol, t)'; 
//      print("\n almin ", almin);  
//      print("\n sumcol almin ",  ((sumc(almin[][1:ncol-1]))/1000)');  
//      print("\n sumrows almin ", ((sumr(almin[][1:ncol-3]))/1000)');  
//      print("\n grandsum almin ", ((sumc(sumr(almin[][1:ncol-3]))+3792238)/1000000)');  
//    ricostruisce i dati per gli anni dal 1884 in poi
      decl almin2=almin[0:21][]|(almin[22][]+almin[23][]/2)|(almin[23:t-2][]/2+almin[24:t-1][]/2)|almin[t-1][];       
//      print("\n almin ", almin/1000);  

// Altri Ministeri 2: Prospetto 10 pg. 1128 e sg.
//   Somme pagate per le diverse specie di opere pubbliche a carico dei vari Ministeri, escluso quello
//dei LLPP Limitatamente a: edifici vari, Bonifiche-irrigazioni e derivazioni d'acqua, e impianti fissi pstelegrafonici e radiotelegrafici
      ncol=3;
      decl al2min    =  loadmat("c:/gauss/storia/destefani-altrimin2.dat", 1);  //
      al2min= shape(al2min, ncol, t)'; 
//      print("\n al2min ", al2min);  
//      print("\n sumcol al2min ",  ((sumc(al2min[][1:ncol-1]))/1000)');  
//      print("\n sumrows al2min ", ((sumr(al2min[][1:ncol-1]))/1000)');  
//    ricostruisce i dati per gli anni dal 1884 in poi
      decl al2min2=al2min[0:21][]|(al2min[22][]+al2min[23][]/2)|(al2min[23:t-2][]/2+al2min[24:t-1][]/2)|al2min[t-1][];       
//      print("\n al2min ", al2min/1000);  

// Amministrazione pubblica: spesa pubblica ecc. 
//   1a colonna: pg 42-43 Destefani, spese pubbliche
      ncol=2;
      decl pubam    =  loadmat("c:/gauss/storia/pubamm.dat", 1);  //
      pubam= shape(pubam, ncol, t)'; 
//      print("\n pubam ", pubam);  
//      print("\n sumcol pubam ",  ((sumc(pubam[][1:ncol-1]))/1000)');  
//      print("\n sumrows pubam ", ((sumr(pubam[][1:ncol-1]))/1000)');  
//    ricostruisce i dati per gli anni dal 1884 in poi
      decl pubam2=pubam[0:21][]|(pubam[22][]+pubam[23][]/2)|(pubam[23:t-2][]/2+pubam[24:t-1][]/2)|pubam[t-1][];       
//     print("\n pubam ", pubam/1000);  




ncol=24;
      decl totrip=strad2+idrau2+bonif2+marit2+edifi2;     

//    attribuzione a totrip di NRIPA: 
  //    trasp. capitale: sino '69: toscana: dal '71: lazio
      totrip[0:8][8]=totrip[0:8][8]+nripa2[0:8][1];
      totrip[9:t-1][11]=totrip[9:t-1][11]+nripa2[9:t-1][1];  
  //    acq pugliese: divisione in base dati pg. 109: 60/62: puglia,; 1/62: campania 1/62: bas.
      totrip[][15]=totrip[][15]+(60/62)*nripa2[][2];
      totrip[][14]=totrip[][14]+(1/62)*nripa2[][2];
      totrip[][16]=totrip[][16]+(1/62)*nripa2[][2];
  //    nuove provincie: 1/3 trieste, 1/3 zara, 1/3 trento
      totrip[][3]=totrip[][3]+nripa2[][3]/3;
      totrip[][21]=totrip[][21]+(1/3)*nripa2[][3];
      totrip[][22]=totrip[][22]+(1/3)*nripa2[][3];

//    attribuzione a totrip di ALMIN2: 
//    ministeri, colonne 1-7 + 9 + 18(mon. vittorio em.) + 20 (min. colonie): Lazio
      totrip[][11]=totrip[][11]+almin2[][1]+almin2[][2]+almin2[][3]+almin2[][4]+almin2[][5]+almin2[][6]+almin2[][7]+almin2[][10]+almin2[][17]+almin2[][19];
//    sicilia: messina
      totrip[][18]=totrip[][18]+almin2[][13];
//    sardegna: cagliari
      totrip[][19]=totrip[][19]+almin2[][14];

// stampa i ripartiti:
//      print("\n sum rip 1871", sumr(nripa2[9][1:3]~almin2[9][1:7]~almin2[9][10]~almin2[9][17]~almin2[9][19]~almin2[9][13:14]));
//      print("\n sum rip 1924", sumr(nripa2[t-1][1:3]~almin2[t-1][1:7]~almin2[t-1][10]~almin2[t-1][17]~almin2[t-1][19]~almin2[t-1][13:14]));
//      print("\n sum rip con strade ecc. 1871", sumr(totrip[9][1:]));
//      print("\n sum rip con strade ecc. 1924", sumr(totrip[t-1][1:]));


//    tot inv non ripartiti 0: spese non ripartite da totrip
//                        1-9:  almin (costruzioni varie, esclusi i ministeri ecc.
//                      10-12:  al2min (ed. vari, bonif. e irrig., postelegrafi)          

      decl totnorip=totrip[][20]~almin2[][8:9]~almin2[][11:12]~almin2[][15:16]~almin2[][18]~almin2[][20:21]~al2min2[][1:2];

     decl anni=strad[0:23][0]|(strad[24:t-1][0]+1);
//   stampa gli aggregati
      decl totinv=sumr(totrip[][1:19]~totrip[][21:22]~totnorip[][]~ferro2[][1:]);

      print("\n Tabella 0. Tot investimenti pubblici dell'Amministrazione centrale e la spesa pubblica. Milioni di lire a prezzi correnti. Anni 1862-1924");
      print("\n totinv con ferro ", totinv/1000000);
      print("\n 353325575 di OP nelle colonie non considerate");

//=>
      print("\n Tabella 1. Gli investimenti pubblici dell'Amministrazione centrale e la spesa pubblica. Milioni di lire a prezzi correnti. Anni 1862-1924");
      print("\n totinv-ferro, ferro, tot., spesa pubblica, ratio totinv-ferro, ratio ferro ", 
      (totinv-sumr(ferro2[][1:]))/1000000
      ~sumr(ferro2[][1:]/1000000)
      ~(totinv)/1000000~pubam2[][1]/1000000~totinv./pubam2[][1]*100
      ~sumr(fenoalt[1:][1:2]));
//);
      print("\n 353325575 di OP nelle colonie non considerate");

      DrawTMatrix(0, totinv[0:54][]', {"IP Amm. centrale dal 1862"}, 1862, 1, 1, 2, 4);
      SetDrawWindow("draw1");
      ShowDrawWindow();

      DrawTMatrix(0, totinv[55:t-1][]', {"IP Amm. centrale dal 1917"}, 1917, 1, 1, 2, 4);
      SetDrawWindow("draw2");
      ShowDrawWindow();

      DrawTMatrix(0, (totinv./pubam2[][1])'*100, {"IP / spesa pubblica"}, 1862, 1, 1, 2, 4);
      SetDrawWindow("draw3");
      ShowDrawWindow();

      DrawTMatrix(0, (sumr(totrip[][1:19]~totrip[][21:22])./(totinv-sumr(ferro2[][1:]))*100)', {"IP ripartibili / IP totali"}, 1862, 1, 1, 2, 4);
      SetDrawWindow("draw4");
      ShowDrawWindow();

// Ora ripartiamo per tipo di bene.
decl strade    =sumr(strad2[][1:]);
decl idrauliche=sumr(idrau2[][1:])+nripa2[][3];
decl bonifiche =sumr(bonif2[][1:])+al2min2[][1];
decl marittime =sumr(marit2[][1:]);
//  aggiungi: trasp capitale, ministeri vari
decl edifici   =sumr(edifi2[][1:])+nripa2[][1]+almin2[][1]+almin2[][2]+almin2[][3]+almin2[][4]+almin2[][5]+almin2[][6]+almin2[][7]+almin2[][10]+almin2[][17]+almin2[][19]+almin2[][13]+almin2[][14];
decl posteleg  =al2min2[][2];
decl totriptip =strade+idrauliche+bonifiche+marittime+edifici+posteleg;
print("\n tot divisibili in tipologia ", anni~totriptip/1000000);

//=>
      print("\n Tabella 2. Gli investimenti pubblici per tipologia di bene. Milioni di lire a prezzi correnti. Anni 1862-1924");
      print("\n Strade, Idrauliche, Bonifiche, Marittime, Edifici, Poste, Telegrafi e Telefoni ", 
      anni~strade/1000000~idrauliche/1000000~bonifiche/1000000~marittime/1000000~edifici/1000000~posteleg/1000000~totriptip/1000000);
      print("\n Tipo di IP sul totale (escluse le ferrovie)"); 

      DrawTMatrix(0, (strade./totriptip*100~bonifiche./totriptip*100~marittime./totriptip*100~(edifici./totriptip)*100~(posteleg./totriptip)*100)', 
{"strade", "bonifiche", "marittime", "edifici", "poste"}, 1862, 1, 1, 2, 4);
      SetDrawWindow("draw5");
      ShowDrawWindow();


// ripartizioni territoriali: 
//  0: Ann,  1: PiV,  2: Lom,  3: Trn,  4: Ven,  5: FVG,  6: Lig,  7: EmR,  8: Tos,  9: Umb
// 10: Mar, 11: Laz, 12: Abr, 13: Mol, 14: Cam, 15: Pug, 16: Bas, 17: Cal, 18: Sic, 19: Sar
// 20: SPS, 21: Trp, 22: Zara
   // Italia settentrionale: Piemonte, Liguria, Lombardia, Tre venezie: 1, 2, 3, 4, 5, 6
   // Italia centrale: Emilia, Toscana, Marche, Umbria, Lazio: 7, 8, 9, 10, 11
   // It mer e ins: Abr e Mol, Camp, Pug, Bas, Sic, Sard: 12, 13, 14, 15, 16, 17, 18, 19

// distribuzione territoriale: (zara e tripoli escluse), No ferrovie
      DrawTMatrix(0, 
(sumr(totrip[][1:6])./(sumr(totrip[][1:19]))*100~sumr(totrip[][7:11])./(sumr(totrip[][1:19]))*100~sumr(totrip[][12:19])./(sumr(totrip[][1:19]))*100)',
{"Nord", "Centro", "Sud e Isole"}, 1862, 1, 1, 2, 4);
      SetDrawWindow("draw6");
      ShowDrawWindow();

// distribuzione territoriale: (zara e tripoli escluse), CON ferrovie
      DrawTMatrix(0, 
((sumr(totrip[][1:6])+ferrorip[][2])./((sumr(totrip[][1:19])+sumr(+ferrorip[][2:4])))*100~(sumr(totrip[][7:11])+ferrorip[][3])./((sumr(totrip[][1:19])+sumr(+ferrorip[][2:4])))*100~(sumr(totrip[][12:19])+ferrorip[][4])./((sumr(totrip[][1:19])+sumr(+ferrorip[][2:4])))*100)',{"Nord", "Centro", "Sud e Isole"}, 1862, 1, 1, 2, 4);
      SetDrawWindow("draw7");
      ShowDrawWindow();

//==> 
      print("\n Tabella 3. IP nelle regioni. Percentuali");
      print("\n PiV, Lom, Trn, Ven, FVG, Lig, EmR, Tos, Umb, Mar, Laz, Abr, Mol, Cam, Pug, Bas, Cal, Sic, Sar");
      print("\n 62-75, 76-85, 86-95, 96-05, 06-15, 16-24"); 
      print("\n " , (sumc(totrip[0:13][1:19]./sumr(totrip[0:13][1:19]))/14)'~
                    (sumc(totrip[14:23][1:19]./sumr(totrip[14:23][1:19]))/10)'~
                    (sumc(totrip[24:33][1:19]./sumr(totrip[24:33][1:19]))/10)'~
                    (sumc(totrip[34:43][1:19]./sumr(totrip[34:43][1:19]))/10)'~
                    (sumc(totrip[44:53][1:19]./sumr(totrip[44:53][1:19]))/10)'~
                    (sumc(totrip[54:62][1:19]./sumr(totrip[54:62][1:19]))/9)'
);

//  controlliamo i dati del NITTI:
// "1862-1897-1898. Tot: 4832 milioni; nord: 1965, centro: 793, sud: 919, isole: 545.
print((sumr(sumc(totrip[0:37][1:6]))+sumr(sumc(ferrorip[0:37][2])))/1000000);
print((sumr(sumc(totrip[0:37][7:11]))+sumr(sumc(ferrorip[0:37][3])))/1000000);
print((sumr(sumc(totrip[0:37][12:19]))+sumr(sumc(ferrorip[0:37][4])))/1000000);
//      print("\n totrip", totrip);

   //  controlliamo i dati del NITTI: ora senza ferrovie
// "1862-1897-1898. Tot: 4832 milioni; nord: 1965, centro: 793, sud: 919, isole: 545.
print((sumr(sumc(totrip[0:37][1:6])))/1000000);
print((sumr(sumc(totrip[0:37][7:11])))/1000000);
print((sumr(sumc(totrip[0:37][12:19])))/1000000);
//      print("\n totrip", totrip);

   print("\n the end");
}