P=catenanm(5); uiopen('/home/accounts/personale/clrmrc90/aa1011/sistemi_stocastici/invariantSIM.m', true); invariantQR(P) ans = 0.0043 0.0520 0.1863 0.5431 0.2143 invariantSIM(P,100) ans = 0 0.0400 0.1800 0.5400 0.2400 invariantSIM(P,1000) ans = 0.0030 0.0800 0.2050 0.5330 0.1790 invariantSIM(P,5000) ans = 0.0040 0.0522 0.1766 0.5462 0.2210 randexp(1/2) ans = 1.8298 randexp([1/2,1/2]) ans = 0.1891 7.6539 randexp([1/2,1/3]) ans = 2.1170 1.8549 ceil(2.1) ans = 3 Narrivi=5 Narrivi = 5 Nservizi=7 Nservizi = 7 lambdaa=randexp(a*Narrivi); ??? Undefined function or variable 'a'. lambdaa=randexp(1*Narrivi); lambdab=randexp(1.5*Nservizi); arrivi=cumsum(lambdaa); servizi=cumsum(lambdab); lambdaa=randexp(1*ones(Narrivi,1)); lambdaa lambdaa = 1.7371 0.3989 2.8353 0.2799 3.1491 lambdab=randexp(1.5*ones(Nservizi,1)); arrivi=cumsum(lambdaa); servizi=cumsum(lambdab); arrivi arrivi = 1.7371 2.1360 4.9713 5.2512 8.4003 servizi servizi = 0.4767 1.0310 4.4108 5.3935 5.6445 5.6843 5.7143 arrivi=[arrivi,ones(5,1)] arrivi = 1.7371 1.0000 2.1360 1.0000 4.9713 1.0000 5.2512 1.0000 8.4003 1.0000 servizi=[servizi,-ones(7,1)] servizi = 0.4767 -1.0000 1.0310 -1.0000 4.4108 -1.0000 5.3935 -1.0000 5.6445 -1.0000 5.6843 -1.0000 5.7143 -1.0000 tempi=[arrivi;servizi] tempi = 1.7371 1.0000 2.1360 1.0000 4.9713 1.0000 5.2512 1.0000 8.4003 1.0000 0.4767 -1.0000 1.0310 -1.0000 4.4108 -1.0000 5.3935 -1.0000 5.6445 -1.0000 5.6843 -1.0000 5.7143 -1.0000 sortrows(tempi) ans = 0.4767 -1.0000 1.0310 -1.0000 1.7371 1.0000 2.1360 1.0000 4.4108 -1.0000 4.9713 1.0000 5.2512 1.0000 5.3935 -1.0000 5.6445 -1.0000 5.6843 -1.0000 5.7143 -1.0000 8.4003 1.0000 simulacodamm1 Elapsed time is 5.642955 seconds. Elapsed time is 1.008877 seconds. simulacodamm1 Elapsed time is 5.355237 seconds. Elapsed time is 0.969183 seconds. simulacodamm1 Elapsed time is 8.206122 seconds. Elapsed time is 1.097887 seconds. simulacodamm1 Elapsed time is 8.811727 seconds. Elapsed time is 1.216369 seconds. diary off