PARAMETER VALUES sigma beta w r 0.5000 0.9990 100.0000 0.0003 the fixed point of coefficients 1.0e+003 * 0.9971 1.3302 0.3337 0.0006 0.0000 -0.0000 ve vuaut 1.0e+004 * 2.0000 1.6759 replacement ratio for different initial values period v=16900 v=16942 v=17000 ---------------------------------------------- 1.0000 0.5101 0.8603 1.4963 2.0000 0.4881 0.8125 1.3873 3.0000 0.4675 0.7685 1.2899 4.0000 0.4482 0.7280 1.2025 5.0000 0.4301 0.6907 1.1237 6.0000 0.4131 0.6562 1.0525 7.0000 0.3971 0.6242 0.9879 8.0000 0.3820 0.5945 0.9291 9.0000 0.3677 0.5669 0.8754 10.0000 0.3542 0.5412 0.8263 11.0000 0.3415 0.5172 0.7812 12.0000 0.3294 0.4948 0.7397 13.0000 0.3180 0.4738 0.7015 14.0000 0.3071 0.4541 0.6662 15.0000 0.2968 0.4357 0.6335 16.0000 0.2870 0.4183 0.6031 17.0000 0.2777 0.4020 0.5749 18.0000 0.2689 0.3866 0.5487 19.0000 0.2604 0.3720 0.5242 20.0000 0.2524 0.3583 0.5013 21.0000 0.2447 0.3454 0.4799 22.0000 0.2374 0.3331 0.4599 23.0000 0.2303 0.3215 0.4410 24.0000 0.2236 0.3104 0.4234 25.0000 0.2172 0.3000 0.4067 PARAMETER VALUES sigma beta w r 0.5000 0.9990 100.0000 0.0003 the fixed point of coefficients 1.0e+003 * 0.9971 1.3302 0.3337 0.0006 0.0000 -0.0000 ve vuaut 1.0e+004 * 2.0000 1.6759 replacement ratio for different initial values period v=16900 v=16942 v=17000 ---------------------------------------------- ??? All matrices on a row in the bracketed expression must have the same number of rows. Error in ==> d:\books\green\hugo\hugo1a.m On line 191 ==> disp([period',rep']); eval( 'if exist(''r'',''var''), mauifunc(r), end', ''); 3.4314e-004 vuprime vuprime = 1.0e+004 * 1.6825 1.6759 1.6840 1.6800 1.6825 1.6853 1.6885 1.6917 1.6949 1.6979 1.7006 1.7029 1.7047 1.7060 1.7066 rep(i,t) ans = 0.1729 rep rep = Columns 1 through 7 0.5101 0.4881 0.4675 0.4482 0.4301 0.4131 0.3971 0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.4963 1.3873 1.2899 1.2025 1.1237 1.0525 0.9879 Columns 8 through 14 0.3820 0.3677 0.3542 0.3415 0.3294 0.3180 0.3071 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.9291 0.8754 0.8263 0.7812 0.7397 0.7015 0.6662 Columns 15 through 21 0.2968 0.2870 0.2777 0.2689 0.2604 0.2524 0.2447 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.6335 0.6031 0.5749 0.5487 0.5242 0.5013 0.4799 Columns 22 through 28 0.2374 0.2303 0.2236 0.2172 0.2111 0.2052 0.1996 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.4599 0.4410 0.4234 0.4067 0.3911 0.3763 0.3623 Columns 29 through 35 0.1942 0.1890 0.1840 0.1792 0.1746 0.1702 0.1659 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.3491 0.3367 0.3249 0.3137 0.3030 0.2929 0.2833 Columns 36 through 42 0.1618 0.1579 0.1541 0.1504 0.1468 0.1434 0.1401 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.2742 0.2655 0.2572 0.2493 0.2418 0.2346 0.2277 Columns 43 through 49 0.1369 0.1339 0.1309 0.1280 0.1252 0.1225 0.1199 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.2211 0.2148 0.2088 0.2030 0.1974 0.1921 0.1870 Columns 50 through 52 0.1174 0.1150 0.1126 0.0000 0.0000 0.0000 0.1821 0.1774 0.1729 rep rep = Columns 1 through 7 0.5101 0.4881 0.4675 0.4482 0.4301 0.4131 0.3971 0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.4963 1.3873 1.2899 1.2025 1.1237 1.0525 0.9879 Columns 8 through 14 0.3820 0.3677 0.3542 0.3415 0.3294 0.3180 0.3071 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.9291 0.8754 0.8263 0.7812 0.7397 0.7015 0.6662 Columns 15 through 21 0.2968 0.2870 0.2777 0.2689 0.2604 0.2524 0.2447 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.6335 0.6031 0.5749 0.5487 0.5242 0.5013 0.4799 Columns 22 through 28 0.2374 0.2303 0.2236 0.2172 0.2111 0.2052 0.1996 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.4599 0.4410 0.4234 0.4067 0.3911 0.3763 0.3623 Columns 29 through 35 0.1942 0.1890 0.1840 0.1792 0.1746 0.1702 0.1659 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.3491 0.3367 0.3249 0.3137 0.3030 0.2929 0.2833 Columns 36 through 42 0.1618 0.1579 0.1541 0.1504 0.1468 0.1434 0.1401 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.2742 0.2655 0.2572 0.2493 0.2418 0.2346 0.2277 Columns 43 through 49 0.1369 0.1339 0.1309 0.1280 0.1252 0.1225 0.1199 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.2211 0.2148 0.2088 0.2030 0.1974 0.1921 0.1870 Columns 50 through 52 0.1174 0.1150 0.1126 0.0000 0.0000 0.0000 0.1821 0.1774 0.1729 plot(rep) rep rep = Columns 1 through 7 0.5101 0.4881 0.4675 0.4482 0.4301 0.4131 0.3971 0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.4963 1.3873 1.2899 1.2025 1.1237 1.0525 0.9879 Columns 8 through 14 0.3820 0.3677 0.3542 0.3415 0.3294 0.3180 0.3071 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.9291 0.8754 0.8263 0.7812 0.7397 0.7015 0.6662 Columns 15 through 21 0.2968 0.2870 0.2777 0.2689 0.2604 0.2524 0.2447 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.6335 0.6031 0.5749 0.5487 0.5242 0.5013 0.4799 Columns 22 through 28 0.2374 0.2303 0.2236 0.2172 0.2111 0.2052 0.1996 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.4599 0.4410 0.4234 0.4067 0.3911 0.3763 0.3623 Columns 29 through 35 0.1942 0.1890 0.1840 0.1792 0.1746 0.1702 0.1659 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.3491 0.3367 0.3249 0.3137 0.3030 0.2929 0.2833 Columns 36 through 42 0.1618 0.1579 0.1541 0.1504 0.1468 0.1434 0.1401 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.2742 0.2655 0.2572 0.2493 0.2418 0.2346 0.2277 Columns 43 through 49 0.1369 0.1339 0.1309 0.1280 0.1252 0.1225 0.1199 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.2211 0.2148 0.2088 0.2030 0.1974 0.1921 0.1870 Columns 50 through 52 0.1174 0.1150 0.1126 0.0000 0.0000 0.0000 0.1821 0.1774 0.1729 plot(rep(2,:)) rep rep = Columns 1 through 7 0.5101 0.4881 0.4675 0.4482 0.4301 0.4131 0.3971 0 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.4963 1.3873 1.2899 1.2025 1.1237 1.0525 0.9879 Columns 8 through 14 0.3820 0.3677 0.3542 0.3415 0.3294 0.3180 0.3071 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.9291 0.8754 0.8263 0.7812 0.7397 0.7015 0.6662 Columns 15 through 21 0.2968 0.2870 0.2777 0.2689 0.2604 0.2524 0.2447 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.6335 0.6031 0.5749 0.5487 0.5242 0.5013 0.4799 Columns 22 through 28 0.2374 0.2303 0.2236 0.2172 0.2111 0.2052 0.1996 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.4599 0.4410 0.4234 0.4067 0.3911 0.3763 0.3623 Columns 29 through 35 0.1942 0.1890 0.1840 0.1792 0.1746 0.1702 0.1659 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.3491 0.3367 0.3249 0.3137 0.3030 0.2929 0.2833 Columns 36 through 42 0.1618 0.1579 0.1541 0.1504 0.1468 0.1434 0.1401 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.2742 0.2655 0.2572 0.2493 0.2418 0.2346 0.2277 Columns 43 through 49 0.1369 0.1339 0.1309 0.1280 0.1252 0.1225 0.1199 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.2211 0.2148 0.2088 0.2030 0.1974 0.1921 0.1870 Columns 50 through 52 0.1174 0.1150 0.1126 0.0000 0.0000 0.0000 0.1821 0.1774 0.1729 vautu ??? Undefined function or variable 'vautu'. who Your variables are: Cv ans fid n sigma vmin T beta flag period t vu1 ZK cinv i r test vuaut ZR coeff iter r0 ve vuprime a0 coeff1 k rep vini w aaut cons m s vmax vuaut vuaut = 1.6759e+004 vini vini = 1.0e+004 * 1.6825 1.6759 1.6840 disp(which('hugo1a')); d:\books\green\hugo\hugo1a.m clear 'd:\books\green\hugo\hugo1a.m' disp(which('hugo1a')); d:\books\green\hugo\hugo1a.m clear 'd:\books\green\hugo\hugo1a.m' plot(rep(1,:)) tic;hugo1a;toc PARAMETER VALUES sigma beta w r 0.5000 0.9990 100.0000 0.0003 the fixed point of coefficients 1.0e+003 * 0.9971 1.3302 0.3337 0.0006 0.0000 -0.0000 ve vuaut 1.0e+004 * 2.0000 1.6759 replacement ratio for different initial values period v=16900 v=16942 v=17000 ---------------------------------------------- 1.0000 0.5101 0.0000 1.4963 2.0000 0.4881 0.0000 1.3873 3.0000 0.4675 0.0000 1.2899 4.0000 0.4482 0.0000 1.2025 5.0000 0.4301 0.0000 1.1237 6.0000 0.4131 0.0000 1.0525 7.0000 0.3971 0.0000 0.9879 8.0000 0.3820 0.0000 0.9291 9.0000 0.3677 0.0000 0.8754 10.0000 0.3542 0.0000 0.8263 11.0000 0.3415 0.0000 0.7812 12.0000 0.3294 0.0000 0.7397 13.0000 0.3180 0.0000 0.7015 14.0000 0.3071 0.0000 0.6662 15.0000 0.2968 0.0000 0.6335 16.0000 0.2870 0.0000 0.6031 17.0000 0.2777 0.0000 0.5749 18.0000 0.2689 0.0000 0.5487 19.0000 0.2604 0.0000 0.5242 20.0000 0.2524 0.0000 0.5013 21.0000 0.2447 0.0000 0.4799 22.0000 0.2374 0.0000 0.4599 23.0000 0.2303 0.0000 0.4410 24.0000 0.2236 0.0000 0.4234 25.0000 0.2172 0.0000 0.4067 26.0000 0.2111 0.0000 0.3911 27.0000 0.2052 0.0000 0.3763 28.0000 0.1996 0.0000 0.3623 29.0000 0.1942 0.0000 0.3491 30.0000 0.1890 0.0000 0.3367 31.0000 0.1840 0.0000 0.3249 32.0000 0.1792 0.0000 0.3137 33.0000 0.1746 0.0000 0.3030 34.0000 0.1702 0.0000 0.2929 35.0000 0.1659 0.0000 0.2833 36.0000 0.1618 0.0000 0.2742 37.0000 0.1579 0.0000 0.2655 38.0000 0.1541 0.0000 0.2572 39.0000 0.1504 0.0000 0.2493 40.0000 0.1468 0.0000 0.2418 41.0000 0.1434 0.0000 0.2346 42.0000 0.1401 0.0000 0.2277 43.0000 0.1369 0.0000 0.2211 44.0000 0.1339 0.0000 0.2148 45.0000 0.1309 0.0000 0.2088 46.0000 0.1280 0.0000 0.2030 47.0000 0.1252 0.0000 0.1974 48.0000 0.1225 0.0000 0.1921 49.0000 0.1199 0.0000 0.1870 50.0000 0.1174 0.0000 0.1821 51.0000 0.1150 0.0000 0.1774 52.0000 0.1126 0.0000 0.1729 PARAMETER VALUES sigma beta w r 0.5000 0.9990 100.0000 0.0003 the fixed point of coefficients 1.0e+003 * 0.9971 1.3302 0.3337 0.0006 0.0000 -0.0000 ve vuaut 1.0e+004 * 2.0000 1.6759 replacement ratio for different initial values period v=16900 v=16942 v=17000 ---------------------------------------------- 1.0000 0.5101 0.0000 1.4963 2.0000 0.4881 0.0000 1.3873 3.0000 0.4675 0.0000 1.2899 4.0000 0.4482 0.0000 1.2025 5.0000 0.4301 0.0000 1.1237 6.0000 0.4131 0.0000 1.0525 7.0000 0.3971 0.0000 0.9879 8.0000 0.3820 0.0000 0.9291 9.0000 0.3677 0.0000 0.8754 10.0000 0.3542 0.0000 0.8263 11.0000 0.3415 0.0000 0.7812 12.0000 0.3294 0.0000 0.7397 13.0000 0.3180 0.0000 0.7015 14.0000 0.3071 0.0000 0.6662 15.0000 0.2968 0.0000 0.6335 16.0000 0.2870 0.0000 0.6031 17.0000 0.2777 0.0000 0.5749 18.0000 0.2689 0.0000 0.5487 19.0000 0.2604 0.0000 0.5242 20.0000 0.2524 0.0000 0.5013 21.0000 0.2447 0.0000 0.4799 22.0000 0.2374 0.0000 0.4599 23.0000 0.2303 0.0000 0.4410 24.0000 0.2236 0.0000 0.4234 25.0000 0.2172 0.0000 0.4067 26.0000 0.2111 0.0000 0.3911 27.0000 0.2052 0.0000 0.3763 28.0000 0.1996 0.0000 0.3623 29.0000 0.1942 0.0000 0.3491 30.0000 0.1890 0.0000 0.3367 31.0000 0.1840 0.0000 0.3249 32.0000 0.1792 0.0000 0.3137 33.0000 0.1746 0.0000 0.3030 34.0000 0.1702 0.0000 0.2929 35.0000 0.1659 0.0000 0.2833 36.0000 0.1618 0.0000 0.2742 37.0000 0.1579 0.0000 0.2655 38.0000 0.1541 0.0000 0.2572 39.0000 0.1504 0.0000 0.2493 40.0000 0.1468 0.0000 0.2418 41.0000 0.1434 0.0000 0.2346 42.0000 0.1401 0.0000 0.2277 43.0000 0.1369 0.0000 0.2211 44.0000 0.1339 0.0000 0.2148 45.0000 0.1309 0.0000 0.2088 46.0000 0.1280 0.0000 0.2030 47.0000 0.1252 0.0000 0.1974 48.0000 0.1225 0.0000 0.1921 49.0000 0.1199 0.0000 0.1870 50.0000 0.1174 0.0000 0.1821 51.0000 0.1150 0.0000 0.1774 52.0000 0.1126 0.0000 0.1729 PARAMETER VALUES sigma beta w r 0.5000 0.9990 100.0000 0.0003 the fixed point of coefficients 1.0e+003 * 0.9971 1.3302 0.3337 0.0006 0.0000 -0.0000 ve vuaut 1.0e+004 * 2.0000 1.6759 replacement ratio for different initial values period v=16900 v=16942 v=17000 ---------------------------------------------- 1.0000 0.5101 0.8603 1.4963 2.0000 0.4881 0.8125 1.3873 3.0000 0.4675 0.7685 1.2899 4.0000 0.4482 0.7280 1.2025 5.0000 0.4301 0.6907 1.1237 6.0000 0.4131 0.6562 1.0525 7.0000 0.3971 0.6242 0.9879 8.0000 0.3820 0.5945 0.9291 9.0000 0.3677 0.5669 0.8754 10.0000 0.3542 0.5412 0.8263 11.0000 0.3415 0.5172 0.7812 12.0000 0.3294 0.4948 0.7397 13.0000 0.3180 0.4738 0.7015 14.0000 0.3071 0.4541 0.6662 15.0000 0.2968 0.4357 0.6335 16.0000 0.2870 0.4183 0.6031 17.0000 0.2777 0.4020 0.5749 18.0000 0.2689 0.3866 0.5487 19.0000 0.2604 0.3720 0.5242 20.0000 0.2524 0.3583 0.5013 21.0000 0.2447 0.3454 0.4799 22.0000 0.2374 0.3331 0.4599 23.0000 0.2303 0.3215 0.4410 24.0000 0.2236 0.3104 0.4234 25.0000 0.2172 0.3000 0.4067 26.0000 0.2111 0.2900 0.3911 27.0000 0.2052 0.2806 0.3763 28.0000 0.1996 0.2716 0.3623 29.0000 0.1942 0.2630 0.3491 30.0000 0.1890 0.2548 0.3367 31.0000 0.1840 0.2470 0.3249 32.0000 0.1792 0.2396 0.3137 33.0000 0.1746 0.2325 0.3030 34.0000 0.1702 0.2257 0.2929 35.0000 0.1659 0.2192 0.2833 36.0000 0.1618 0.2130 0.2742 37.0000 0.1579 0.2070 0.2655 38.0000 0.1541 0.2013 0.2572 39.0000 0.1504 0.1958 0.2493 40.0000 0.1468 0.1906 0.2418 41.0000 0.1434 0.1855 0.2346 42.0000 0.1401 0.1807 0.2277 43.0000 0.1369 0.1760 0.2211 44.0000 0.1339 0.1715 0.2148 45.0000 0.1309 0.1672 0.2088 46.0000 0.1280 0.1631 0.2030 47.0000 0.1252 0.1591 0.1974 48.0000 0.1225 0.1552 0.1921 49.0000 0.1199 0.1515 0.1870 50.0000 0.1174 0.1479 0.1821 51.0000 0.1150 0.1445 0.1774 52.0000 0.1126 0.1411 0.1729