07-Jan-2022 15:57:50 bdf2_test(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Test bdf2(). predator_prey_bdf2(): A pair of ordinary differential equations for a population of predators and prey are solved using the BDF2, solving the implicit function using fsolve(). The exact solution shows periodic behavior, with a fixed period and amplitude. t1 = 0 y1 = 5000 100 t2 = 0.0250 t3 = 0.0500 t1 = 0.0250 y1 = 1.0e+03 * 5.2432 0.1006 t2 = 0.0500 t3 = 0.0750 t1 = 0.0500 y1 = 1.0e+03 * 5.4982 0.1025 t2 = 0.0750 t3 = 0.1000 t1 = 0.0750 y1 = 1.0e+03 * 5.7652 0.1059 t2 = 0.1000 t3 = 0.1250 t1 = 0.1000 y1 = 1.0e+03 * 6.0445 0.1108 t2 = 0.1250 t3 = 0.1500 t1 = 0.1250 y1 = 1.0e+03 * 6.3364 0.1177 t2 = 0.1500 t3 = 0.1750 t1 = 0.1500 y1 = 1.0e+03 * 6.6410 0.1270 t2 = 0.1750 t3 = 0.2000 t1 = 0.1750 y1 = 1.0e+03 * 6.9583 0.1391 t2 = 0.2000 t3 = 0.2250 t1 = 0.2000 y1 = 1.0e+03 * 7.2882 0.1550 t2 = 0.2250 t3 = 0.2500 t1 = 0.2250 y1 = 1.0e+03 * 7.6302 0.1757 t2 = 0.2500 t3 = 0.2750 t1 = 0.2500 y1 = 1.0e+03 * 7.9833 0.2028 t2 = 0.2750 t3 = 0.3000 t1 = 0.2750 y1 = 1.0e+03 * 8.3461 0.2384 t2 = 0.3000 t3 = 0.3250 t1 = 0.3000 y1 = 1.0e+03 * 8.7160 0.2856 t2 = 0.3250 t3 = 0.3500 t1 = 0.3250 y1 = 1.0e+03 * 9.0895 0.3489 t2 = 0.3500 t3 = 0.3750 t1 = 0.3500 y1 = 1.0e+03 * 9.4608 0.4347 t2 = 0.3750 t3 = 0.4000 t1 = 0.3750 y1 = 1.0e+03 * 9.8216 0.5520 t2 = 0.4000 t3 = 0.4250 t1 = 0.4000 y1 = 1.0e+04 * 1.0160 0.0714 t2 = 0.4250 t3 = 0.4500 t1 = 0.4250 y1 = 1.0e+04 * 1.0457 0.0939 t2 = 0.4500 t3 = 0.4750 t1 = 0.4500 y1 = 1.0e+04 * 1.0690 0.1251 t2 = 0.4750 t3 = 0.5000 t1 = 0.4750 y1 = 1.0e+04 * 1.0825 0.1681 t2 = 0.5000 t3 = 0.5250 t1 = 0.5000 y1 = 1.0e+04 * 1.0823 0.2264 t2 = 0.5250 t3 = 0.5500 t1 = 0.5250 y1 = 1.0e+04 * 1.0640 0.3027 t2 = 0.5500 t3 = 0.5750 t1 = 0.5500 y1 = 1.0e+04 * 1.0241 0.3977 t2 = 0.5750 t3 = 0.6000 t1 = 0.5750 y1 = 1.0e+03 * 9.6158 5.0735 t2 = 0.6000 t3 = 0.6250 t1 = 0.6000 y1 = 1.0e+03 * 8.7884 6.2254 t2 = 0.6250 t3 = 0.6500 t1 = 0.6250 y1 = 1.0e+03 * 7.8223 7.2956 t2 = 0.6500 t3 = 0.6750 t1 = 0.6500 y1 = 1.0e+03 * 6.8037 8.1419 t2 = 0.6750 t3 = 0.7000 t1 = 0.6750 y1 = 1.0e+03 * 5.8184 8.6603 t2 = 0.7000 t3 = 0.7250 t1 = 0.7000 y1 = 1.0e+03 * 4.9301 8.8125 t2 = 0.7250 t3 = 0.7500 t1 = 0.7250 y1 = 1.0e+03 * 4.1732 8.6255 t2 = 0.7500 t3 = 0.7750 t1 = 0.7500 y1 = 1.0e+03 * 3.5553 8.1698 t2 = 0.7750 t3 = 0.8000 t1 = 0.7750 y1 = 1.0e+03 * 3.0666 7.5324 t2 = 0.8000 t3 = 0.8250 t1 = 0.8000 y1 = 1.0e+03 * 2.6887 6.7964 t2 = 0.8250 t3 = 0.8500 t1 = 0.8250 y1 = 1.0e+03 * 2.4015 6.0288 t2 = 0.8500 t3 = 0.8750 t1 = 0.8500 y1 = 1.0e+03 * 2.1864 5.2780 t2 = 0.8750 t3 = 0.9000 t1 = 0.8750 y1 = 1.0e+03 * 2.0276 4.5745 t2 = 0.9000 t3 = 0.9250 t1 = 0.9000 y1 = 1.0e+03 * 1.9130 3.9350 t2 = 0.9250 t3 = 0.9500 t1 = 0.9250 y1 = 1.0e+03 * 1.8331 3.3665 t2 = 0.9500 t3 = 0.9750 t1 = 0.9500 y1 = 1.0e+03 * 1.7806 2.8691 t2 = 0.9750 t3 = 1 t1 = 0.9750 y1 = 1.0e+03 * 1.7503 2.4391 t2 = 1 t3 = 1.0250 t1 = 1 y1 = 1.0e+03 * 1.7382 2.0706 t2 = 1.0250 t3 = 1.0500 t1 = 1.0250 y1 = 1.0e+03 * 1.7412 1.7570 t2 = 1.0500 t3 = 1.0750 t1 = 1.0500 y1 = 1.0e+03 * 1.7571 1.4912 t2 = 1.0750 t3 = 1.1000 t1 = 1.0750 y1 = 1.0e+03 * 1.7842 1.2668 t2 = 1.1000 t3 = 1.1250 t1 = 1.1000 y1 = 1.0e+03 * 1.8212 1.0778 t2 = 1.1250 t3 = 1.1500 t1 = 1.1250 y1 = 1.0e+03 * 1.8672 0.9189 t2 = 1.1500 t3 = 1.1750 t1 = 1.1500 y1 = 1.0e+03 * 1.9214 0.7853 t2 = 1.1750 t3 = 1.2000 t1 = 1.1750 y1 = 1.0e+03 * 1.9834 0.6731 t2 = 1.2000 t3 = 1.2250 t1 = 1.2000 y1 = 1.0e+03 * 2.0527 0.5788 t2 = 1.2250 t3 = 1.2500 t1 = 1.2250 y1 = 1.0e+03 * 2.1290 0.4996 t2 = 1.2500 t3 = 1.2750 t1 = 1.2500 y1 = 1.0e+03 * 2.2122 0.4329 t2 = 1.2750 t3 = 1.3000 t1 = 1.2750 y1 = 1.0e+03 * 2.3022 0.3768 t2 = 1.3000 t3 = 1.3250 t1 = 1.3000 y1 = 1.0e+03 * 2.3990 0.3296 t2 = 1.3250 t3 = 1.3500 t1 = 1.3250 y1 = 1.0e+03 * 2.5026 0.2897 t2 = 1.3500 t3 = 1.3750 t1 = 1.3500 y1 = 1.0e+03 * 2.6131 0.2560 t2 = 1.3750 t3 = 1.4000 t1 = 1.3750 y1 = 1.0e+03 * 2.7306 0.2276 t2 = 1.4000 t3 = 1.4250 t1 = 1.4000 y1 = 1.0e+03 * 2.8552 0.2036 t2 = 1.4250 t3 = 1.4500 t1 = 1.4250 y1 = 1.0e+03 * 2.9872 0.1833 t2 = 1.4500 t3 = 1.4750 t1 = 1.4500 y1 = 1.0e+03 * 3.1268 0.1661 t2 = 1.4750 t3 = 1.5000 t1 = 1.4750 y1 = 1.0e+03 * 3.2742 0.1517 t2 = 1.5000 t3 = 1.5250 t1 = 1.5000 y1 = 1.0e+03 * 3.4297 0.1396 t2 = 1.5250 t3 = 1.5500 t1 = 1.5250 y1 = 1.0e+03 * 3.5935 0.1295 t2 = 1.5500 t3 = 1.5750 t1 = 1.5500 y1 = 1.0e+03 * 3.7661 0.1211 t2 = 1.5750 t3 = 1.6000 t1 = 1.5750 y1 = 1.0e+03 * 3.9477 0.1143 t2 = 1.6000 t3 = 1.6250 t1 = 1.6000 y1 = 1.0e+03 * 4.1386 0.1089 t2 = 1.6250 t3 = 1.6500 t1 = 1.6250 y1 = 1.0e+03 * 4.3394 0.1048 t2 = 1.6500 t3 = 1.6750 t1 = 1.6500 y1 = 1.0e+03 * 4.5502 0.1019 t2 = 1.6750 t3 = 1.7000 t1 = 1.6750 y1 = 1.0e+03 * 4.7716 0.1002 t2 = 1.7000 t3 = 1.7250 t1 = 1.7000 y1 = 1.0e+03 * 5.0039 0.0996 t2 = 1.7250 t3 = 1.7500 t1 = 1.7250 y1 = 1.0e+03 * 5.2474 0.1003 t2 = 1.7500 t3 = 1.7750 t1 = 1.7500 y1 = 1.0e+03 * 5.5027 0.1022 t2 = 1.7750 t3 = 1.8000 t1 = 1.7750 y1 = 1.0e+03 * 5.7699 0.1055 t2 = 1.8000 t3 = 1.8250 t1 = 1.8000 y1 = 1.0e+03 * 6.0495 0.1105 t2 = 1.8250 t3 = 1.8500 t1 = 1.8250 y1 = 1.0e+03 * 6.3417 0.1174 t2 = 1.8500 t3 = 1.8750 t1 = 1.8500 y1 = 1.0e+03 * 6.6466 0.1267 t2 = 1.8750 t3 = 1.9000 t1 = 1.8750 y1 = 1.0e+03 * 6.9643 0.1389 t2 = 1.9000 t3 = 1.9250 t1 = 1.9000 y1 = 1.0e+03 * 7.2945 0.1548 t2 = 1.9250 t3 = 1.9500 t1 = 1.9250 y1 = 1.0e+03 * 7.6368 0.1755 t2 = 1.9500 t3 = 1.9750 t1 = 1.9500 y1 = 1.0e+03 * 7.9903 0.2026 t2 = 1.9750 t3 = 2 t1 = 1.9750 y1 = 1.0e+03 * 8.3534 0.2383 t2 = 2 t3 = 2.0250 t1 = 2 y1 = 1.0e+03 * 8.7237 0.2856 t2 = 2.0250 t3 = 2.0500 t1 = 2.0250 y1 = 1.0e+03 * 9.0974 0.3491 t2 = 2.0500 t3 = 2.0750 t1 = 2.0500 y1 = 1.0e+03 * 9.4690 0.4350 t2 = 2.0750 t3 = 2.1000 t1 = 2.0750 y1 = 1.0e+03 * 9.8300 0.5526 t2 = 2.1000 t3 = 2.1250 t1 = 2.1000 y1 = 1.0e+04 * 1.0168 0.0715 t2 = 2.1250 t3 = 2.1500 t1 = 2.1250 y1 = 1.0e+04 * 1.0466 0.0940 t2 = 2.1500 t3 = 2.1750 t1 = 2.1500 y1 = 1.0e+04 * 1.0698 0.1254 t2 = 2.1750 t3 = 2.2000 t1 = 2.1750 y1 = 1.0e+04 * 1.0832 0.1686 t2 = 2.2000 t3 = 2.2250 t1 = 2.2000 y1 = 1.0e+04 * 1.0828 0.2271 t2 = 2.2250 t3 = 2.2500 t1 = 2.2250 y1 = 1.0e+04 * 1.0642 0.3038 t2 = 2.2500 t3 = 2.2750 t1 = 2.2500 y1 = 1.0e+04 * 1.0241 0.3990 t2 = 2.2750 t3 = 2.3000 t1 = 2.2750 y1 = 1.0e+03 * 9.6119 5.0900 t2 = 2.3000 t3 = 2.3250 t1 = 2.3000 y1 = 1.0e+03 * 8.7813 6.2436 t2 = 2.3250 t3 = 2.3500 t1 = 2.3250 y1 = 1.0e+03 * 7.8125 7.3138 t2 = 2.3500 t3 = 2.3750 t1 = 2.3500 y1 = 1.0e+03 * 6.7924 8.1579 t2 = 2.3750 t3 = 2.4000 t1 = 2.3750 y1 = 1.0e+03 * 5.8067 8.6725 t2 = 2.4000 t3 = 2.4250 t1 = 2.4000 y1 = 1.0e+03 * 4.9189 8.8202 t2 = 2.4250 t3 = 2.4500 t1 = 2.4250 y1 = 1.0e+03 * 4.1631 8.6287 t2 = 2.4500 t3 = 2.4750 t1 = 2.4500 y1 = 1.0e+03 * 3.5465 8.1691 t2 = 2.4750 t3 = 2.5000 t1 = 2.4750 y1 = 1.0e+03 * 3.0591 7.5288 t2 = 2.5000 t3 = 2.5250 t1 = 2.5000 y1 = 1.0e+03 * 2.6824 6.7908 t2 = 2.5250 t3 = 2.5500 t1 = 2.5250 y1 = 1.0e+03 * 2.3963 6.0221 t2 = 2.5500 t3 = 2.5750 t1 = 2.5500 y1 = 1.0e+03 * 2.1819 5.2707 t2 = 2.5750 t3 = 2.6000 t1 = 2.5750 y1 = 1.0e+03 * 2.0239 4.5672 t2 = 2.6000 t3 = 2.6250 t1 = 2.6000 y1 = 1.0e+03 * 1.9098 3.9281 t2 = 2.6250 t3 = 2.6500 t1 = 2.6250 y1 = 1.0e+03 * 1.8303 3.3600 t2 = 2.6500 t3 = 2.6750 t1 = 2.6500 y1 = 1.0e+03 * 1.7782 2.8632 t2 = 2.6750 t3 = 2.7000 t1 = 2.6750 y1 = 1.0e+03 * 1.7482 2.4338 t2 = 2.7000 t3 = 2.7250 t1 = 2.7000 y1 = 1.0e+03 * 1.7363 2.0659 t2 = 2.7250 t3 = 2.7500 t1 = 2.7250 y1 = 1.0e+03 * 1.7395 1.7528 t2 = 2.7500 t3 = 2.7750 t1 = 2.7500 y1 = 1.0e+03 * 1.7556 1.4875 t2 = 2.7750 t3 = 2.8000 t1 = 2.7750 y1 = 1.0e+03 * 1.7828 1.2636 t2 = 2.8000 t3 = 2.8250 t1 = 2.8000 y1 = 1.0e+03 * 1.8200 1.0750 t2 = 2.8250 t3 = 2.8500 t1 = 2.8250 y1 = 1.0e+03 * 1.8660 0.9164 t2 = 2.8500 t3 = 2.8750 t1 = 2.8500 y1 = 1.0e+03 * 1.9203 0.7832 t2 = 2.8750 t3 = 2.9000 t1 = 2.8750 y1 = 1.0e+03 * 1.9823 0.6712 t2 = 2.9000 t3 = 2.9250 t1 = 2.9000 y1 = 1.0e+03 * 2.0517 0.5772 t2 = 2.9250 t3 = 2.9500 t1 = 2.9250 y1 = 1.0e+03 * 2.1280 0.4981 t2 = 2.9500 t3 = 2.9750 t1 = 2.9500 y1 = 1.0e+03 * 2.2113 0.4317 t2 = 2.9750 t3 = 3 t1 = 2.9750 y1 = 1.0e+03 * 2.3014 0.3757 t2 = 3 t3 = 3.0250 t1 = 3 y1 = 1.0e+03 * 2.3982 0.3286 t2 = 3.0250 t3 = 3.0500 t1 = 3.0250 y1 = 1.0e+03 * 2.5018 0.2888 t2 = 3.0500 t3 = 3.0750 t1 = 3.0500 y1 = 1.0e+03 * 2.6123 0.2552 t2 = 3.0750 t3 = 3.1000 t1 = 3.0750 y1 = 1.0e+03 * 2.7298 0.2269 t2 = 3.1000 t3 = 3.1250 t1 = 3.1000 y1 = 1.0e+03 * 2.8545 0.2029 t2 = 3.1250 t3 = 3.1500 t1 = 3.1250 y1 = 1.0e+03 * 2.9865 0.1827 t2 = 3.1500 t3 = 3.1750 t1 = 3.1500 y1 = 1.0e+03 * 3.1261 0.1656 t2 = 3.1750 t3 = 3.2000 t1 = 3.1750 y1 = 1.0e+03 * 3.2735 0.1512 t2 = 3.2000 t3 = 3.2250 t1 = 3.2000 y1 = 1.0e+03 * 3.4289 0.1391 t2 = 3.2250 t3 = 3.2500 t1 = 3.2250 y1 = 1.0e+03 * 3.5928 0.1290 t2 = 3.2500 t3 = 3.2750 t1 = 3.2500 y1 = 1.0e+03 * 3.7654 0.1207 t2 = 3.2750 t3 = 3.3000 t1 = 3.2750 y1 = 1.0e+03 * 3.9470 0.1139 t2 = 3.3000 t3 = 3.3250 t1 = 3.3000 y1 = 1.0e+03 * 4.1379 0.1085 t2 = 3.3250 t3 = 3.3500 t1 = 3.3250 y1 = 1.0e+03 * 4.3387 0.1044 t2 = 3.3500 t3 = 3.3750 t1 = 3.3500 y1 = 1.0e+03 * 4.5495 0.1016 t2 = 3.3750 t3 = 3.4000 t1 = 3.3750 y1 = 1.0e+03 * 4.7709 0.0998 t2 = 3.4000 t3 = 3.4250 t1 = 3.4000 y1 = 1.0e+03 * 5.0032 0.0993 t2 = 3.4250 t3 = 3.4500 t1 = 3.4250 y1 = 1.0e+03 * 5.2468 0.0999 t2 = 3.4500 t3 = 3.4750 t1 = 3.4500 y1 = 1.0e+03 * 5.5020 0.1018 t2 = 3.4750 t3 = 3.5000 t1 = 3.4750 y1 = 1.0e+03 * 5.7693 0.1051 t2 = 3.5000 t3 = 3.5250 t1 = 3.5000 y1 = 1.0e+03 * 6.0490 0.1101 t2 = 3.5250 t3 = 3.5500 t1 = 3.5250 y1 = 1.0e+03 * 6.3412 0.1170 t2 = 3.5500 t3 = 3.5750 t1 = 3.5500 y1 = 1.0e+03 * 6.6461 0.1262 t2 = 3.5750 t3 = 3.6000 t1 = 3.5750 y1 = 1.0e+03 * 6.9639 0.1383 t2 = 3.6000 t3 = 3.6250 t1 = 3.6000 y1 = 1.0e+03 * 7.2942 0.1542 t2 = 3.6250 t3 = 3.6500 t1 = 3.6250 y1 = 1.0e+03 * 7.6366 0.1748 t2 = 3.6500 t3 = 3.6750 t1 = 3.6500 y1 = 1.0e+03 * 7.9902 0.2018 t2 = 3.6750 t3 = 3.7000 t1 = 3.6750 y1 = 1.0e+03 * 8.3535 0.2373 t2 = 3.7000 t3 = 3.7250 t1 = 3.7000 y1 = 1.0e+03 * 8.7240 0.2845 t2 = 3.7250 t3 = 3.7500 t1 = 3.7250 y1 = 1.0e+03 * 9.0980 0.3477 t2 = 3.7500 t3 = 3.7750 t1 = 3.7500 y1 = 1.0e+03 * 9.4700 0.4334 t2 = 3.7750 t3 = 3.8000 t1 = 3.7750 y1 = 1.0e+03 * 9.8315 0.5506 t2 = 3.8000 t3 = 3.8250 t1 = 3.8000 y1 = 1.0e+04 * 1.0170 0.0712 t2 = 3.8250 t3 = 3.8500 t1 = 3.8250 y1 = 1.0e+04 * 1.0469 0.0937 t2 = 3.8500 t3 = 3.8750 t1 = 3.8500 y1 = 1.0e+04 * 1.0702 0.1250 t2 = 3.8750 t3 = 3.9000 t1 = 3.8750 y1 = 1.0e+04 * 1.0837 0.1681 t2 = 3.9000 t3 = 3.9250 t1 = 3.9000 y1 = 1.0e+04 * 1.0835 0.2265 t2 = 3.9250 t3 = 3.9500 t1 = 3.9250 y1 = 1.0e+04 * 1.0651 0.3031 t2 = 3.9500 t3 = 3.9750 t1 = 3.9500 y1 = 1.0e+04 * 1.0251 0.3983 t2 = 3.9750 t3 = 4 t1 = 3.9750 y1 = 1.0e+03 * 9.6226 5.0838 t2 = 4 t3 = 4.0250 t1 = 4 y1 = 1.0e+03 * 8.7921 6.2394 t2 = 4.0250 t3 = 4.0500 t1 = 4.0250 y1 = 1.0e+03 * 7.8226 7.3126 t2 = 4.0500 t3 = 4.0750 t1 = 4.0500 y1 = 1.0e+03 * 6.8012 8.1603 t2 = 4.0750 t3 = 4.1000 t1 = 4.0750 y1 = 1.0e+03 * 5.8136 8.6782 t2 = 4.1000 t3 = 4.1250 t1 = 4.1000 y1 = 1.0e+03 * 4.9240 8.8285 t2 = 4.1250 t3 = 4.1500 t1 = 4.1250 y1 = 1.0e+03 * 4.1665 8.6386 t2 = 4.1500 t3 = 4.1750 t1 = 4.1500 y1 = 1.0e+03 * 3.5485 8.1795 t2 = 4.1750 t3 = 4.2000 t1 = 4.1750 y1 = 1.0e+03 * 3.0600 7.5390 t2 = 4.2000 t3 = 4.2250 t1 = 4.2000 y1 = 1.0e+03 * 2.6826 6.8002 t2 = 4.2250 t3 = 4.2500 t1 = 4.2250 y1 = 1.0e+03 * 2.3958 6.0305 t2 = 4.2500 t3 = 4.2750 t1 = 4.2500 y1 = 1.0e+03 * 2.1811 5.2779 t2 = 4.2750 t3 = 4.3000 t1 = 4.2750 y1 = 1.0e+03 * 2.0228 4.5733 t2 = 4.3000 t3 = 4.3250 t1 = 4.3000 y1 = 1.0e+03 * 1.9085 3.9331 t2 = 4.3250 t3 = 4.3500 t1 = 4.3250 y1 = 1.0e+03 * 1.8288 3.3641 t2 = 4.3500 t3 = 4.3750 t1 = 4.3500 y1 = 1.0e+03 * 1.7766 2.8664 t2 = 4.3750 t3 = 4.4000 t1 = 4.3750 y1 = 1.0e+03 * 1.7465 2.4363 t2 = 4.4000 t3 = 4.4250 t1 = 4.4000 y1 = 1.0e+03 * 1.7345 2.0679 t2 = 4.4250 t3 = 4.4500 t1 = 4.4250 y1 = 1.0e+03 * 1.7376 1.7543 t2 = 4.4500 t3 = 4.4750 t1 = 4.4500 y1 = 1.0e+03 * 1.7536 1.4887 t2 = 4.4750 t3 = 4.5000 t1 = 4.4750 y1 = 1.0e+03 * 1.7808 1.2644 t2 = 4.5000 t3 = 4.5250 t1 = 4.5000 y1 = 1.0e+03 * 1.8178 1.0756 t2 = 4.5250 t3 = 4.5500 t1 = 4.5250 y1 = 1.0e+03 * 1.8639 0.9168 t2 = 4.5500 t3 = 4.5750 t1 = 4.5500 y1 = 1.0e+03 * 1.9181 0.7834 t2 = 4.5750 t3 = 4.6000 t1 = 4.5750 y1 = 1.0e+03 * 1.9800 0.6714 t2 = 4.6000 t3 = 4.6250 t1 = 4.6000 y1 = 1.0e+03 * 2.0492 0.5772 t2 = 4.6250 t3 = 4.6500 t1 = 4.6250 y1 = 1.0e+03 * 2.1255 0.4981 t2 = 4.6500 t3 = 4.6750 t1 = 4.6500 y1 = 1.0e+03 * 2.2087 0.4316 t2 = 4.6750 t3 = 4.7000 t1 = 4.6750 y1 = 1.0e+03 * 2.2986 0.3756 t2 = 4.7000 t3 = 4.7250 t1 = 4.7000 y1 = 1.0e+03 * 2.3953 0.3284 t2 = 4.7250 t3 = 4.7500 t1 = 4.7250 y1 = 1.0e+03 * 2.4989 0.2886 t2 = 4.7500 t3 = 4.7750 t1 = 4.7500 y1 = 1.0e+03 * 2.6092 0.2550 t2 = 4.7750 t3 = 4.8000 t1 = 4.7750 y1 = 1.0e+03 * 2.7266 0.2267 t2 = 4.8000 t3 = 4.8250 t1 = 4.8000 y1 = 1.0e+03 * 2.8512 0.2027 t2 = 4.8250 t3 = 4.8500 t1 = 4.8250 y1 = 1.0e+03 * 2.9830 0.1824 t2 = 4.8500 t3 = 4.8750 t1 = 4.8500 y1 = 1.0e+03 * 3.1225 0.1653 t2 = 4.8750 t3 = 4.9000 t1 = 4.8750 y1 = 1.0e+03 * 3.2697 0.1509 t2 = 4.9000 t3 = 4.9250 t1 = 4.9000 y1 = 1.0e+03 * 3.4250 0.1388 t2 = 4.9250 t3 = 4.9500 t1 = 4.9250 y1 = 1.0e+03 * 3.5887 0.1288 t2 = 4.9500 t3 = 4.9750 t1 = 4.9500 y1 = 1.0e+03 * 3.7611 0.1204 t2 = 4.9750 t3 = 5 Graphics saved as "predator_prey_bdf2.png" predator_prey_bdf2(): Normal end of execution. stiff_bdf2() Use bdf2() to solve the stiff ODE. t1 = 0 y1 = 0 t2 = 0.0370 t3 = 0.0741 t1 = 0.0370 y1 = 0.9614 t2 = 0.0741 t3 = 0.1111 t1 = 0.0741 y1 = 1.1246 t2 = 0.1111 t3 = 0.1481 t1 = 0.1111 y1 = 1.0767 t2 = 0.1481 t3 = 0.1852 t1 = 0.1481 y1 = 1.0211 t2 = 0.1852 t3 = 0.2222 t1 = 0.1852 y1 = 0.9917 t2 = 0.2222 t3 = 0.2593 t1 = 0.2222 y1 = 0.9783 t2 = 0.2593 t3 = 0.2963 t1 = 0.2593 y1 = 0.9698 t2 = 0.2963 t3 = 0.3333 t1 = 0.2963 y1 = 0.9612 t2 = 0.3333 t3 = 0.3704 t1 = 0.3333 y1 = 0.9509 t2 = 0.3704 t3 = 0.4074 t1 = 0.3704 y1 = 0.9390 t2 = 0.4074 t3 = 0.4444 t1 = 0.4074 y1 = 0.9257 t2 = 0.4444 t3 = 0.4815 t1 = 0.4444 y1 = 0.9111 t2 = 0.4815 t3 = 0.5185 t1 = 0.4815 y1 = 0.8952 t2 = 0.5185 t3 = 0.5556 t1 = 0.5185 y1 = 0.8781 t2 = 0.5556 t3 = 0.5926 t1 = 0.5556 y1 = 0.8598 t2 = 0.5926 t3 = 0.6296 t1 = 0.5926 y1 = 0.8403 t2 = 0.6296 t3 = 0.6667 t1 = 0.6296 y1 = 0.8197 t2 = 0.6667 t3 = 0.7037 t1 = 0.6667 y1 = 0.7979 t2 = 0.7037 t3 = 0.7407 t1 = 0.7037 y1 = 0.7751 t2 = 0.7407 t3 = 0.7778 t1 = 0.7407 y1 = 0.7512 t2 = 0.7778 t3 = 0.8148 t1 = 0.7778 y1 = 0.7262 t2 = 0.8148 t3 = 0.8519 t1 = 0.8148 y1 = 0.7003 t2 = 0.8519 t3 = 0.8889 t1 = 0.8519 y1 = 0.6734 t2 = 0.8889 t3 = 0.9259 t1 = 0.8889 y1 = 0.6456 t2 = 0.9259 t3 = 0.9630 t1 = 0.9259 y1 = 0.6168 t2 = 0.9630 t3 = 1 Graphics saved as "stiff_bdf2.png" stiff_bdf2() Normal end of execution. bdf2_test(): Normal end of execution. 07-Jan-2022 15:57:56