//STORIA1.OX // // 28 aprile 2002 confindustria#2 // prende dati dal 1928 al 1995 da P. Montano // // // #include #include #include "c:\dpdox\dpd.h" #include "c:\dpdox\dpd.ox" main() { decl dpd=new DPD(); decl i; decl nprov; decl t = 68 ; // dal 1928 al 1995 //Stradali Ferroviarie Edilizia Idrauliche Altre TOTALE // // pie decl pie = loadmat("c:/dpdox/storia/dati/pie.dat", 1); // pie=shape(pie, 7, t)'; // vaa decl vaa = loadmat("c:/dpdox/storia/dati/vaa.dat", 1); // vaa=shape(vaa, 7, t)'; // print("\n vaa", vaa); decl piv=pie+vaa; // lom decl lom = loadmat("c:/dpdox/storia/dati/lom.dat", 1); // lom=shape(lom, 7, t)'; // print("\n lom", lom); // taa decl taa = loadmat("c:/dpdox/storia/dati/taa.dat", 1); // taa=shape(taa, 7, t)'; // print("\n taa", taa); // ven decl ven = loadmat("c:/dpdox/storia/dati/ven.dat", 1); // ven=shape(ven, 7, t)'; // print("\n ven", ven); // fvg decl fvg = loadmat("c:/dpdox/storia/dati/fvg.dat", 1); // fvg=shape(fvg, 7, t)'; // print("\n fvg", fvg); // lig decl lig = loadmat("c:/dpdox/storia/dati/lig.dat", 1); // lig=shape(lig, 7, t)'; // print("\n lig", lig); // emr decl emr = loadmat("c:/dpdox/storia/dati/emr.dat", 1); // emr=shape(emr, 7, t)'; // print("\n emr", emr); // tos decl tos = loadmat("c:/dpdox/storia/dati/tos.dat", 1); // tos=shape(tos, 7, t)'; // print("\n tos", tos); // mar decl mar = loadmat("c:/dpdox/storia/dati/mar.dat", 1); // mar=shape(mar, 7, t)'; // print("\n mar", mar); // umb decl umb = loadmat("c:/dpdox/storia/dati/umb.dat", 1); // umb=shape(umb, 7, t)'; // print("\n umb", umb); // laz decl laz = loadmat("c:/dpdox/storia/dati/laz.dat", 1); // laz=shape(laz, 7, t)'; // print("\n laz", laz); // cam decl cam = loadmat("c:/dpdox/storia/dati/cam.dat", 1); // cam=shape(cam, 7, t)'; // print("\n cam", cam); // abm decl abm = loadmat("c:/dpdox/storia/dati/abm.dat", 1); // abm=shape(abm, 7, t)'; // print("\n abm", abm); // pug decl pug = loadmat("c:/dpdox/storia/dati/pug.dat", 1); // pug=shape(pug, 7, t)'; // print("\n pug", pug); // bas decl bas = loadmat("c:/dpdox/storia/dati/bas.dat", 1); // bas=shape(bas, 7, t)'; // print("\n bas", bas); // cal decl cal = loadmat("c:/dpdox/storia/dati/cal.dat", 1); // cal=shape(cal, 7, t)'; // print("\n cal", cal); // sic decl sic = loadmat("c:/dpdox/storia/dati/sic.dat", 1); // sic=shape(sic, 7, t)'; // print("\n sic", sic); // sar decl sar = loadmat("c:/dpdox/storia/dati/sar.dat", 1); // sar=shape(sar, 7, t)'; // print("\n sar", sar); // ita decl ita = loadmat("c:/dpdox/storia/dati/ita.dat", 1); // ita=shape(ita, 7, t)'; print("\n ita", ita); // ip prezzi correnti, tab 4 rossi, 1890-1998 decl ip = loadmat("c:/rossi/tab4.dat", 1); // ip=shape(ip, 9, 108)'; ip=ip[][6]; // print("\n ip", ip[38:105][]); decl ita2=piv+taa+lom+ven+fvg+lig+emr+tos+mar+umb+laz+abm+cam+pug+bas+cal+sic+sar; //print("\n ita2 " , ita2); // tabella con le frazioni per regione: decl ipreg = piv[][6].*ip[38:105]~lom[][6].*ip[38:105]~taa[][6].*ip[38:105]~ ven[][6].*ip[38:105]~fvg[][6].*ip[38:105]~lig[][6].*ip[38:105]~ emr[][6].*ip[38:105]~tos[][6].*ip[38:105]~umb[][6].*ip[38:105]~ mar[][6].*ip[38:105]~laz[][6].*ip[38:105]~abm[][6].*ip[38:105]~ cam[][6].*ip[38:105]~pug[][6].*ip[38:105]~bas[][6].*ip[38:105]~ cal[][6].*ip[38:105]~sic[][6].*ip[38:105]~sar[][6].*ip[38:105]; // print("\ipreg ", (ipreg)); print("\sumr ipreg ", sumr(ipreg)); // Grafico: tipo di bene: frazione del totale. dal '28 al '52 DrawTMatrix(0, (ita[0:24][1:5])'*100, {"Strade", "Ferrovie", "Edifici pubblici", "Idrauliche", "Altro"}, 1928, 1, 1, 2, 4); SetDrawWindow("draw1"); ShowDrawWindow(); //Tabella: frazione del totale per ciascuna regione print("\n Tabella ?. IP nelle regioni. Percentuali"); print("\n PiV, Lom, Trn, Ven, FVG, Lig, EmR, Tos, Umb, Mar, Laz, Abm, Cam, Pug, Bas, Cal, Sic, Sar"); print("\n 28-32, 33-37, 38-42, 43-47, 48-52"); print("\n " , (sumc(ipreg[0:4][0:17]./sumr(ipreg[0:4][0:17]))/5)'~ (sumc(ipreg[5:9][0:17]./sumr(ipreg[5:9][0:17]))/5)'~ (sumc(ipreg[10:14][0:17]./sumr(ipreg[10:14][0:17]))/5)'~ (sumc(ipreg[15:19][0:17]./sumr(ipreg[15:19][0:17]))/5)'~ (sumc(ipreg[20:24][0:17]./sumr(ipreg[20:24][0:17]))/5)' ); // Grafico: distribuzione per macroregione. dal '28 al '52 // Italia settentrionale: Piemonte, Liguria, Lombardia, Tre venezie: 0, 1, 2, 3, 4, 5, // Italia centrale: Emilia, Toscana, Marche, Umbria, Lazio: 6, 7, 8, 9, 10, // It mer e ins: Abr e Mol, Camp, Pug, Bas, Sic, Sard: 11, 12, 13, 14, 15, 16, 17 DrawTMatrix(0, ((sumr(ipreg[0:24][0:5])~sumr(ipreg[0:24][6:10])~sumr(ipreg[0:24][11:17]))./sumr(ipreg[0:24][0:17]))'*100 , {"Nord", "Centro", "Sud e Isole"}, 1928, 1, 1, 2, 4); SetDrawWindow("draw2"); ShowDrawWindow(); // 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)); // DrawTMatrix(0, totinv[0:54][]', {"IP Amm. centrale dal 1862"}, 1862, 1, 1, 2, 4); // SetDrawWindow("draw1"); // ShowDrawWindow(); print("\n the end"); }