GNU Octave, version 2.1.40 (i386-redhat-linux-gnu). Copyright (C) 1996, 1997, 1998, 1999, 2000, 2001, 2002 John W. Eaton. This is free software; see the source code for copying conditions. There is ABSOLUTELY NO WARRANTY; not even for MERCHANTIBILITY or FITNESS FOR A PARTICULAR PURPOSE. For details, type `warranty'. Please contribute if you find this software useful. For more information, visit http://www.octave.org/help-wanted.html Report bugs to . octave> lsode("hill", [1,0, .1, .2], 2*pi) ans = 1.00000 0.00000 0.10000 0.20000 -0.49020 -0.14802 0.10000 0.20000 octave> lsode("hill", [1,0, .1, .2], [0,2*pi]) ans = 1.00000 0.00000 0.10000 0.20000 octave> lsode("hill", [1,0, .1, .2], [0,2*pi]) ans = 1.00000 0.00000 0.10000 0.20000 -0.49020 -0.14802 0.10000 0.20000 octave> ans(2,:) ans = -0.49020 -0.14802 0.10000 0.20000 octave> hillstrob(.1,.2) ans = -0.49020 5.13253 -0.14802 -0.49020 octave> sol = lsode("hill", [1, 0, .1, .2], linspace(0, 10, 100)); error: `linspce' undefined near line 9 column 37 error: evaluating argument list element number 3 error: evaluating assignment expression near line 9, column 5 octave> sol = lsode("hill", [1, 0, .1, .2], linspace(0, 10, 100))(:, 1); octave> sol sol = 1.00000 0.00000 0.10000 0.20000 0.99873 -0.02522 0.10000 0.20000 0.99491 -0.05025 0.10000 0.20000 0.98859 -0.07489 0.10000 0.20000 0.97980 -0.09897 0.10000 0.20000 0.96862 -0.12230 0.10000 0.20000 0.95512 -0.14472 0.10000 0.20000 0.93942 -0.16606 0.10000 0.20000 0.92162 -0.18619 0.10000 0.20000 0.90185 -0.20498 0.10000 0.20000 0.88026 -0.22233 0.10000 0.20000 0.85699 -0.23813 0.10000 0.20000 0.83220 -0.25232 0.10000 0.20000 0.80607 -0.26485 0.10000 0.20000 0.77875 -0.27570 0.10000 0.20000 0.75043 -0.28487 0.10000 0.20000 0.72126 -0.29238 0.10000 0.20000 0.69142 -0.29826 0.10000 0.20000 0.66106 -0.30259 0.10000 0.20000 0.63034 -0.30546 0.10000 0.20000 0.59940 -0.30697 0.10000 0.20000 0.56836 -0.30724 0.10000 0.20000 0.53736 -0.30642 0.10000 0.20000 0.50649 -0.30467 0.10000 0.20000 0.47584 -0.30213 0.10000 0.20000 0.44548 -0.29898 0.10000 0.20000 0.41545 -0.29539 0.10000 0.20000 0.38581 -0.29152 0.10000 0.20000 0.35657 -0.28751 0.10000 0.20000 0.32773 -0.28353 0.10000 0.20000 0.29928 -0.27969 0.10000 0.20000 0.27121 -0.27611 0.10000 0.20000 0.24349 -0.27287 0.10000 0.20000 0.21608 -0.27004 0.10000 0.20000 0.18892 -0.26765 0.10000 0.20000 0.16199 -0.26572 0.10000 0.20000 0.13523 -0.26424 0.10000 0.20000 0.10859 -0.26318 0.10000 0.20000 0.08205 -0.26250 0.10000 0.20000 0.05555 -0.26211 0.10000 0.20000 0.02909 -0.26195 0.10000 0.20000 0.00263 -0.26192 0.10000 0.20000 -0.02383 -0.26191 0.10000 0.20000 -0.05028 -0.26183 0.10000 0.20000 -0.07672 -0.26157 0.10000 0.20000 -0.10311 -0.26103 0.10000 0.20000 -0.12944 -0.26010 0.10000 0.20000 -0.15564 -0.25870 0.10000 0.20000 -0.18168 -0.25673 0.10000 0.20000 -0.20749 -0.25413 0.10000 0.20000 -0.23300 -0.25083 0.10000 0.20000 -0.25813 -0.24677 0.10000 0.20000 -0.28282 -0.24193 0.10000 0.20000 -0.30698 -0.23626 0.10000 0.20000 -0.33052 -0.22976 0.10000 0.20000 -0.35337 -0.22243 0.10000 0.20000 -0.37543 -0.21428 0.10000 0.20000 -0.39663 -0.20534 0.10000 0.20000 -0.41689 -0.19564 0.10000 0.20000 -0.43613 -0.18525 0.10000 0.20000 -0.45429 -0.17422 0.10000 0.20000 -0.47131 -0.16262 0.10000 0.20000 -0.48713 -0.15053 0.10000 0.20000 -0.50170 -0.13805 0.10000 0.20000 -0.51500 -0.12527 0.10000 0.20000 -0.52700 -0.11229 0.10000 0.20000 -0.53769 -0.09922 0.10000 0.20000 -0.54705 -0.08617 0.10000 0.20000 -0.55510 -0.07325 0.10000 0.20000 -0.56185 -0.06057 0.10000 0.20000 -0.56735 -0.04825 0.10000 0.20000 -0.57162 -0.03639 0.10000 0.20000 -0.57471 -0.02509 0.10000 0.20000 -0.57671 -0.01447 0.10000 0.20000 -0.57766 -0.00460 0.10000 0.20000 -0.57767 0.00441 0.10000 0.20000 -0.57680 0.01251 0.10000 0.20000 -0.57517 0.01962 0.10000 0.20000 -0.57287 0.02568 0.10000 0.20000 -0.57002 0.03066 0.10000 0.20000 -0.56672 0.03452 0.10000 0.20000 -0.56308 0.03725 0.10000 0.20000 -0.55923 0.03885 0.10000 0.20000 -0.55527 0.03934 0.10000 0.20000 -0.55132 0.03873 0.10000 0.20000 -0.54748 0.03708 0.10000 0.20000 -0.54386 0.03446 0.10000 0.20000 -0.54055 0.03093 0.10000 0.20000 -0.53764 0.02658 0.10000 0.20000 -0.53521 0.02152 0.10000 0.20000 -0.53331 0.01587 0.10000 0.20000 -0.53202 0.00975 0.10000 0.20000 -0.53135 0.00329 0.10000 0.20000 -0.53136 -0.00336 0.10000 0.20000 -0.53203 -0.01006 0.10000 0.20000 -0.53339 -0.01667 0.10000 0.20000 -0.53540 -0.02304 0.10000 0.20000 -0.53803 -0.02904 0.10000 0.20000 -0.54124 -0.03452 0.10000 0.20000 -0.54498 -0.03936 0.10000 0.20000 octave> sol = sol(:, 1:2) sol = 1.00000 0.00000 0.99873 -0.02522 0.99491 -0.05025 0.98859 -0.07489 0.97980 -0.09897 0.96862 -0.12230 0.95512 -0.14472 0.93942 -0.16606 0.92162 -0.18619 0.90185 -0.20498 0.88026 -0.22233 0.85699 -0.23813 0.83220 -0.25232 0.80607 -0.26485 0.77875 -0.27570 0.75043 -0.28487 0.72126 -0.29238 0.69142 -0.29826 0.66106 -0.30259 0.63034 -0.30546 0.59940 -0.30697 0.56836 -0.30724 0.53736 -0.30642 0.50649 -0.30467 0.47584 -0.30213 0.44548 -0.29898 0.41545 -0.29539 0.38581 -0.29152 0.35657 -0.28751 0.32773 -0.28353 0.29928 -0.27969 0.27121 -0.27611 0.24349 -0.27287 0.21608 -0.27004 0.18892 -0.26765 0.16199 -0.26572 0.13523 -0.26424 0.10859 -0.26318 0.08205 -0.26250 0.05555 -0.26211 0.02909 -0.26195 0.00263 -0.26192 -0.02383 -0.26191 -0.05028 -0.26183 -0.07672 -0.26157 -0.10311 -0.26103 -0.12944 -0.26010 -0.15564 -0.25870 -0.18168 -0.25673 -0.20749 -0.25413 -0.23300 -0.25083 -0.25813 -0.24677 -0.28282 -0.24193 -0.30698 -0.23626 -0.33052 -0.22976 -0.35337 -0.22243 -0.37543 -0.21428 -0.39663 -0.20534 -0.41689 -0.19564 -0.43613 -0.18525 -0.45429 -0.17422 -0.47131 -0.16262 -0.48713 -0.15053 -0.50170 -0.13805 -0.51500 -0.12527 -0.52700 -0.11229 -0.53769 -0.09922 -0.54705 -0.08617 -0.55510 -0.07325 -0.56185 -0.06057 -0.56735 -0.04825 -0.57162 -0.03639 -0.57471 -0.02509 -0.57671 -0.01447 -0.57766 -0.00460 -0.57767 0.00441 -0.57680 0.01251 -0.57517 0.01962 -0.57287 0.02568 -0.57002 0.03066 -0.56672 0.03452 -0.56308 0.03725 -0.55923 0.03885 -0.55527 0.03934 -0.55132 0.03873 -0.54748 0.03708 -0.54386 0.03446 -0.54055 0.03093 -0.53764 0.02658 -0.53521 0.02152 -0.53331 0.01587 -0.53202 0.00975 -0.53135 0.00329 -0.53136 -0.00336 -0.53203 -0.01006 -0.53339 -0.01667 -0.53540 -0.02304 -0.53803 -0.02904 -0.54124 -0.03452 -0.54498 -0.03936 octave> sol = sol(:, 1) sol = 1.0000000 0.9987254 0.9949115 0.9885872 0.9798007 0.9686184 0.9551243 0.9394187 0.9216173 0.9018490 0.8802552 0.8569868 0.8322031 0.8060693 0.7787543 0.7504283 0.7212605 0.6914169 0.6610581 0.6303365 0.5993950 0.5683646 0.5373627 0.5064922 0.4758399 0.4454761 0.4154542 0.3858109 0.3565665 0.3277266 0.2992827 0.2712143 0.2434911 0.2160751 0.1889232 0.1619894 0.1352274 0.1085931 0.0820463 0.0555529 0.0290864 0.0026289 -0.0238274 -0.0502803 -0.0767169 -0.1031137 -0.1294371 -0.1556435 -0.1816805 -0.2074875 -0.2329968 -0.2581348 -0.2828234 -0.3069811 -0.3305244 -0.3533690 -0.3754316 -0.3966308 -0.4168885 -0.4361313 -0.4542912 -0.4713074 -0.4871265 -0.5017040 -0.5150047 -0.5270034 -0.5376856 -0.5470477 -0.5550971 -0.5618529 -0.5673452 -0.5716152 -0.5747150 -0.5767070 -0.5776635 -0.5776656 -0.5768026 -0.5751713 -0.5728744 -0.5700194 -0.5667178 -0.5630832 -0.5592299 -0.5552716 -0.5513198 -0.5474822 -0.5438609 -0.5405515 -0.5376407 -0.5352058 -0.5333127 -0.5320153 -0.5313544 -0.5313568 -0.5320349 -0.5333866 -0.5353950 -0.5380291 -0.5412438 -0.5449810 octave> plot(sol) octave> plot(lsode("hill", [1, 0, .1, .2], time)(:, 1)); octave> plot(sol) octave> time=linspace(0, 100, 100); octave> plot(time, lsode("hill", [1, 0, .1, .2], time)(:, 1)); octave> plot(time, lsode("hill", [1, 0, .1, .2], time)(:, 1)); octave> time=linspace(0, 10, 100); octave> plot(time, lsode("hill", [1, 0, .2, .3], time)(:, 1)); octave> U0 = hillstrob(.1, .2); error: `hillstob' undefined near line 21 column 6 error: evaluating assignment expression near line 21, column 4 octave> U0 = hillstrob(.1, .2); octave> U0 U0 = -0.49020 5.13253 -0.14802 -0.49020 octave> poly(U0) ans = 1.00000 0.98039 1.00000 octave> hilltrace([.1,.2],[.3,.4]) ans = 1.2217 1.6195 2.4444 2.9747 octave> plot(time, lsode("hill", [1, 0, .2, .3], time)(:, 1)); octave> time=linspace(0, 20, 100); octave> plot(time, lsode("hill", [1, 0, .2, .3], time)(:, 1)); octave> time=linspace(0, 30, 100); octave> plot(time, lsode("hill", [1, 0, .2, .3], time)(:, 1)); octave> time=linspace(0, 30, 100); octave> plot(time, lsode("hill", [1, 0, .2, .3], time)(:, 1)); octave> hill_unstable octave> hill_unstable octave>