Home License -- for personal use only. Not for government, academic, research, commercial, or other organizational use. 04-Nov-2026 19:04:47 beetle_dynamic_test(): MATLAB/Octave version 9.11.0.2358333 (R2021b) Update 7 Test beetle_dynamic(). beetle_parameters_test(): Test beetle_parameters(), which stores, updates, and retrieves parameter values for the beetle dynamnical system. Current internal parameter values: b = 2 cel = 0.01 cea = 0.02 cpa = 0.03 mul = 0.45 mua = 0.55 colony = 1 Updated internal parameter values: b = 10 cel = 0.02 cea = 0.04 cpa = 0.06 mul = 0.08 mua = 0.1 colony = 2 beetle_history_test(): Simulate beetle dynamics with three values of mu_a Graphics saved as "history_0.1.png" Graphics saved as "history_0.6.png" Graphics saved as "history_0.9.png" beetle_equilibrium_test(): Test beetle_equilibrium(), which computes certain equilibrium solutions of the beetle dynamical system, assuming cel = 0 and mua / ( 1 - mul ) <= b. b L P A 1.0 0 0 0 1.5 25 12 20 2.0 49 25 35 2.5 72 36 46 3.0 95 48 55 3.5 117 59 63 4.0 139 69 69 4.5 160 80 75 5.0 180 90 80 5.5 200 100 85 6.0 219 110 90 6.5 239 119 94 7.0 257 129 97 7.5 276 138 101 8.0 294 147 104 8.5 312 156 107 9.0 330 165 110 9.5 347 173 113 10.0 364 182 115 10.5 381 190 118 11.0 398 199 120 11.5 414 207 122 12.0 430 215 124 12.5 446 223 126 13.0 462 231 128 13.5 478 239 130 14.0 494 247 132 14.5 509 255 134 15.0 524 262 135 15.5 540 270 137 16.0 555 277 139 16.5 569 285 140 17.0 584 292 142 17.5 599 299 143 18.0 613 307 145 18.5 627 314 146 19.0 642 321 147 19.5 656 328 149 20.0 670 335 150 Graphics saved as "beetle_equilibrium.png" beetle_picard_test(): Test beetle_picard(), which for a range of values b seeks a nonzero fixed point of the beetle dynamic system. b F1 F2 F3 1.0 0 0 0 1.5 15 7 13 2.0 27 13 21 2.5 36 18 28 3.0 45 22 32 3.5 52 26 36 4.0 59 30 40 4.5 65 33 43 5.0 71 35 45 5.5 76 38 47 6.0 81 40 49 6.5 85 43 51 7.0 89 45 53 7.5 93 47 54 8.0 97 48 56 8.5 100 50 57 9.0 104 52 58 9.5 107 53 59 10.0 110 55 60 10.5 113 56 61 11.0 116 58 62 11.5 118 59 63 12.0 121 60 64 12.5 123 62 65 13.0 126 63 65 13.5 128 64 66 14.0 130 65 67 14.5 132 66 67 15.0 135 67 68 15.5 137 68 69 16.0 139 69 69 16.5 141 70 70 17.0 142 71 70 17.5 144 72 71 18.0 146 73 71 18.5 148 74 72 19.0 150 75 72 19.5 151 76 73 20.0 153 76 73 Graphics saved as "beetle_picard.png" beetle_bifurcation_test(): A bifurcation diagram for the beetle dynamical system based on parameter mua. Bifurcations will show up as multiple dots for a given value of mua. Graphics saved as "beetle_bifurcation" beetle_fit_function_test(): Test beetle_fit_function(), which evaluates the difference between observed data depending on unknown parameter values, and data generated by an estimate for those parameters. Test 1 Norm of discrepancy is 17.1945 Expected norm is 17.19 Test 2 Norm of discrepancy is 5.03677 Expected norm is 5.04 beetle_fit_test(): Use a nonlinear least squares solver to estimate the parameter values that produce a beetle dynamic system history that most closely matches an observed data set. Initial parameter values: b = 10 cel = 0.1 cea = 0.1 cpa = 0.1 mul = 0.5 mua = 0.5 colony = 1 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 4.99304 cel = 0.0184015 cea = 0.00766623 cpa = 0.0176303 mul = 0.165309 mua = 0.0604008 colony = 1 Residual = 25.3522 Initial parameter values: b = 10 cel = 0.1 cea = 0.1 cpa = 0.1 mul = 0.5 mua = 0.5 colony = 2 Solver stopped prematurely. lsqnonlin stopped because it exceeded the function evaluation limit, options.MaxFunctionEvaluations = 6.000000e+02. Estimated parameter values: b = 8.1884 cel = 0.021527 cea = 0.0075564 cpa = 0.0179633 mul = 0.334269 mua = 0.0486299 colony = 2 Residual = 34.8717 Initial parameter values: b = 10 cel = 0.1 cea = 0.1 cpa = 0.1 mul = 0.5 mua = 0.5 colony = 3 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 6.95395 cel = 0.0188466 cea = 0.00805866 cpa = 0.0350094 mul = 0.276799 mua = 8.19287e-08 colony = 3 Residual = 19.3744 Initial parameter values: b = 10 cel = 0.1 cea = 0.1 cpa = 0.1 mul = 0.5 mua = 0.5 colony = 4 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 6.59159 cel = 0.0172515 cea = 0.0127046 cpa = 0.0252853 mul = 0.295285 mua = 4.64834e-05 colony = 4 Residual = 44.3577 beetle_homotopy_test(): Search for the value of parameter b for which the other parameters produce a beetle dynamic system solution most closely matching the observed data for beetle colony #2. Initial parameter values: b = 0.5 cel = 0.1 cea = 0.1 cpa = 0.1 mul = 0.5 mua = 0.5 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 0.5 cel = 1.72542e-07 cea = 2.24435e-12 cpa = 0.00633495 mul = 0.264259 mua = 0.110052 colony = 2 square root of residual = 7.92705 Initial parameter values: b = 1 cel = 1.72542e-07 cea = 2.24435e-12 cpa = 0.00633495 mul = 0.264259 mua = 0.110052 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 1 cel = 7.97527e-10 cea = 6.52036e-08 cpa = 0.0122379 mul = 0.517691 mua = 0.108314 colony = 2 square root of residual = 7.36622 Initial parameter values: b = 1.5 cel = 7.97527e-10 cea = 6.52036e-08 cpa = 0.0122379 mul = 0.517691 mua = 0.108314 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 1.5 cel = 4.53993e-08 cea = 0.00320397 cpa = 0.0120439 mul = 0.559187 mua = 0.115773 colony = 2 square root of residual = 7.31046 Initial parameter values: b = 2 cel = 4.53993e-08 cea = 0.00320397 cpa = 0.0120439 mul = 0.559187 mua = 0.115773 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 2 cel = 0.00505865 cea = 1.12089e-11 cpa = 0.013203 mul = 0.574856 mua = 0.107441 colony = 2 square root of residual = 7.28652 Initial parameter values: b = 2.5 cel = 0.00505865 cea = 1.12089e-11 cpa = 0.013203 mul = 0.574856 mua = 0.107441 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 2.5 cel = 0.00702284 cea = 4.01066e-14 cpa = 0.00165477 mul = 0.5654 mua = 0.387687 colony = 2 square root of residual = 7.25447 Initial parameter values: b = 3 cel = 0.00702284 cea = 4.01066e-14 cpa = 0.00165477 mul = 0.5654 mua = 0.387687 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 3 cel = 0.00889946 cea = 3.92018e-10 cpa = 9.42996e-05 mul = 0.559916 mua = 0.449428 colony = 2 square root of residual = 7.20846 Initial parameter values: b = 3.5 cel = 0.00889946 cea = 3.92018e-10 cpa = 9.42996e-05 mul = 0.559916 mua = 0.449428 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 3.5 cel = 0.0101931 cea = 6.50842e-10 cpa = 6.18529e-11 mul = 0.506122 mua = 0.524126 colony = 2 square root of residual = 7.13541 Initial parameter values: b = 4 cel = 0.0101931 cea = 6.50842e-10 cpa = 6.18529e-11 mul = 0.506122 mua = 0.524126 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 4 cel = 0.01322 cea = 2.46319e-14 cpa = 4.15636e-11 mul = 0.122579 mua = 0.795291 colony = 2 square root of residual = 6.13724 Initial parameter values: b = 4.5 cel = 0.01322 cea = 2.46319e-14 cpa = 4.15636e-11 mul = 0.122579 mua = 0.795291 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 4.5 cel = 0.0124175 cea = 1.15602e-10 cpa = 6.37925e-11 mul = 0.213177 mua = 0.789971 colony = 2 square root of residual = 5.86628 Initial parameter values: b = 5 cel = 0.0124175 cea = 1.15602e-10 cpa = 6.37925e-11 mul = 0.213177 mua = 0.789971 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 5 cel = 0.011681 cea = 6.44836e-12 cpa = 1.36442e-12 mul = 0.288208 mua = 0.78539 colony = 2 square root of residual = 5.66422 Initial parameter values: b = 5.5 cel = 0.011681 cea = 6.44836e-12 cpa = 1.36442e-12 mul = 0.288208 mua = 0.78539 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 5.5 cel = 0.0110036 cea = 4.5592e-13 cpa = 3.98019e-14 mul = 0.350474 mua = 0.781803 colony = 2 square root of residual = 5.52182 Initial parameter values: b = 6 cel = 0.0110036 cea = 4.5592e-13 cpa = 3.98019e-14 mul = 0.350474 mua = 0.781803 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 6 cel = 0.0103891 cea = 2.42685e-10 cpa = 3.0341e-10 mul = 0.40255 mua = 0.779229 colony = 2 square root of residual = 5.42768 Initial parameter values: b = 6.5 cel = 0.0103891 cea = 2.42685e-10 cpa = 3.0341e-10 mul = 0.40255 mua = 0.779229 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 6.5 cel = 0.00984206 cea = 1.56856e-10 cpa = 1.30205e-10 mul = 0.446456 mua = 0.77761 colony = 2 square root of residual = 5.37131 Initial parameter values: b = 7 cel = 0.00984206 cea = 1.56856e-10 cpa = 1.30205e-10 mul = 0.446456 mua = 0.77761 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 7 cel = 0.00935983 cea = 9.82358e-12 cpa = 5.22028e-13 mul = 0.483831 mua = 0.776828 colony = 2 square root of residual = 5.34423 Initial parameter values: b = 7.5 cel = 0.00935983 cea = 9.82358e-12 cpa = 5.22028e-13 mul = 0.483831 mua = 0.776828 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 7.5 cel = 0.00893219 cea = 3.33138e-11 cpa = 1.67457e-12 mul = 0.515997 mua = 0.776663 colony = 2 square root of residual = 5.34002 Initial parameter values: b = 8 cel = 0.00893219 cea = 3.33138e-11 cpa = 1.67457e-12 mul = 0.515997 mua = 0.776663 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 8 cel = 0.00856072 cea = 2.6153e-10 cpa = 2.99268e-12 mul = 0.543931 mua = 0.777095 colony = 2 square root of residual = 5.35391 Initial parameter values: b = 8.5 cel = 0.00856072 cea = 2.6153e-10 cpa = 2.99268e-12 mul = 0.543931 mua = 0.777095 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 8.5 cel = 0.00823433 cea = 9.39675e-11 cpa = 1.68972e-12 mul = 0.56848 mua = 0.777938 colony = 2 square root of residual = 5.3824 Initial parameter values: b = 9 cel = 0.00823433 cea = 9.39675e-11 cpa = 1.68972e-12 mul = 0.56848 mua = 0.777938 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 9 cel = 0.00795106 cea = 1.59016e-10 cpa = 1.00723e-12 mul = 0.59029 mua = 0.779087 colony = 2 square root of residual = 5.4229 Initial parameter values: b = 9.5 cel = 0.00795106 cea = 1.59016e-10 cpa = 1.00723e-12 mul = 0.59029 mua = 0.779087 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 9.5 cel = 0.00770342 cea = 1.78381e-10 cpa = 1.42496e-10 mul = 0.60989 mua = 0.780403 colony = 2 square root of residual = 5.47346 Initial parameter values: b = 10 cel = 0.00770342 cea = 1.78381e-10 cpa = 1.42496e-10 mul = 0.60989 mua = 0.780403 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 10 cel = 0.00748654 cea = 2.46061e-10 cpa = 1.12186e-10 mul = 0.627679 mua = 0.781792 colony = 2 square root of residual = 5.53258 Initial parameter values: b = 10.5 cel = 0.00748654 cea = 2.46061e-10 cpa = 1.12186e-10 mul = 0.627679 mua = 0.781792 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 10.5 cel = 0.00729636 cea = 4.05728e-11 cpa = 2.39402e-12 mul = 0.643956 mua = 0.78319 colony = 2 square root of residual = 5.59904 Initial parameter values: b = 11 cel = 0.00729636 cea = 4.05728e-11 cpa = 2.39402e-12 mul = 0.643956 mua = 0.78319 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 11 cel = 0.00712628 cea = 1.8765e-09 cpa = 7.81036e-12 mul = 0.658981 mua = 0.784506 colony = 2 square root of residual = 5.67183 Initial parameter values: b = 11.5 cel = 0.00712628 cea = 1.8765e-09 cpa = 7.81036e-12 mul = 0.658981 mua = 0.784506 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 11.5 cel = 0.00697309 cea = 6.23469e-10 cpa = 2.23532e-14 mul = 0.672912 mua = 0.785699 colony = 2 square root of residual = 5.75008 Initial parameter values: b = 12 cel = 0.00697309 cea = 6.23469e-10 cpa = 2.23532e-14 mul = 0.672912 mua = 0.785699 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 12 cel = 0.00683852 cea = 0.000206121 cpa = 1.90164e-12 mul = 0.679156 mua = 0.790937 colony = 2 square root of residual = 5.83072 Initial parameter values: b = 12.5 cel = 0.00683852 cea = 0.000206121 cpa = 1.90164e-12 mul = 0.679156 mua = 0.790937 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 12.5 cel = 0.0067106 cea = 0.000458597 cpa = 1.38613e-10 mul = 0.683246 mua = 0.797076 colony = 2 square root of residual = 5.90833 Initial parameter values: b = 13 cel = 0.0067106 cea = 0.000458597 cpa = 1.38613e-10 mul = 0.683246 mua = 0.797076 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 13 cel = 0.00658715 cea = 0.000699319 cpa = 9.95938e-11 mul = 0.68726 mua = 0.802885 colony = 2 square root of residual = 5.98299 Initial parameter values: b = 13.5 cel = 0.00658715 cea = 0.000699319 cpa = 9.95938e-11 mul = 0.68726 mua = 0.802885 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 13.5 cel = 0.00646781 cea = 0.000927411 cpa = 1.28823e-09 mul = 0.691257 mua = 0.808316 colony = 2 square root of residual = 6.05499 Initial parameter values: b = 14 cel = 0.00646781 cea = 0.000927411 cpa = 1.28823e-09 mul = 0.691257 mua = 0.808316 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 14 cel = 0.00635267 cea = 0.00114434 cpa = 1.201e-09 mul = 0.695173 mua = 0.813479 colony = 2 square root of residual = 6.12456 Initial parameter values: b = 14.5 cel = 0.00635267 cea = 0.00114434 cpa = 1.201e-09 mul = 0.695173 mua = 0.813479 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 14.5 cel = 0.00624168 cea = 0.00135144 cpa = 1.51567e-09 mul = 0.699007 mua = 0.818392 colony = 2 square root of residual = 6.19193 Initial parameter values: b = 15 cel = 0.00624168 cea = 0.00135144 cpa = 1.51567e-09 mul = 0.699007 mua = 0.818392 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 15 cel = 0.0061348 cea = 0.00155013 cpa = 1.57652e-09 mul = 0.702749 mua = 0.823076 colony = 2 square root of residual = 6.25725 Initial parameter values: b = 15.5 cel = 0.0061348 cea = 0.00155013 cpa = 1.57652e-09 mul = 0.702749 mua = 0.823076 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 15.5 cel = 0.0060321 cea = 0.00174131 cpa = 9.00445e-11 mul = 0.70639 mua = 0.827575 colony = 2 square root of residual = 6.3207 Initial parameter values: b = 16 cel = 0.0060321 cea = 0.00174131 cpa = 9.00445e-11 mul = 0.70639 mua = 0.827575 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 16 cel = 0.00593307 cea = 0.00192483 cpa = 9.81306e-11 mul = 0.709967 mua = 0.831848 colony = 2 square root of residual = 6.38242 Initial parameter values: b = 16.5 cel = 0.00593307 cea = 0.00192483 cpa = 9.81306e-11 mul = 0.709967 mua = 0.831848 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 16.5 cel = 0.00583753 cea = 0.00210127 cpa = 1.08931e-10 mul = 0.71346 mua = 0.835941 colony = 2 square root of residual = 6.44252 Initial parameter values: b = 17 cel = 0.00583753 cea = 0.00210127 cpa = 1.08931e-10 mul = 0.71346 mua = 0.835941 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 17 cel = 0.00574552 cea = 0.00227175 cpa = 1.16358e-10 mul = 0.716866 mua = 0.839873 colony = 2 square root of residual = 6.50112 Initial parameter values: b = 17.5 cel = 0.00574552 cea = 0.00227175 cpa = 1.16358e-10 mul = 0.716866 mua = 0.839873 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 17.5 cel = 0.0056567 cea = 0.00243627 cpa = 2.76221e-13 mul = 0.720194 mua = 0.843645 colony = 2 square root of residual = 6.5583 Initial parameter values: b = 18 cel = 0.0056567 cea = 0.00243627 cpa = 2.76221e-13 mul = 0.720194 mua = 0.843645 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 18 cel = 0.00557162 cea = 0.00259746 cpa = 2.10151e-12 mul = 0.723404 mua = 0.847318 colony = 2 square root of residual = 6.61416 Initial parameter values: b = 18.5 cel = 0.00557162 cea = 0.00259746 cpa = 2.10151e-12 mul = 0.723404 mua = 0.847318 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 18.5 cel = 0.00548912 cea = 0.00275282 cpa = 2.26243e-09 mul = 0.726545 mua = 0.85085 colony = 2 square root of residual = 6.66877 Initial parameter values: b = 19 cel = 0.00548912 cea = 0.00275282 cpa = 2.26243e-09 mul = 0.726545 mua = 0.85085 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 19 cel = 0.00541046 cea = 0.00290629 cpa = 1.26025e-10 mul = 0.729586 mua = 0.854292 colony = 2 square root of residual = 6.72219 Initial parameter values: b = 19.5 cel = 0.00541046 cea = 0.00290629 cpa = 1.26025e-10 mul = 0.729586 mua = 0.854292 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 19.5 cel = 0.00533425 cea = 0.00305527 cpa = 1.09355e-10 mul = 0.732546 mua = 0.857625 colony = 2 square root of residual = 6.77448 Initial parameter values: b = 20 cel = 0.00533425 cea = 0.00305527 cpa = 1.09355e-10 mul = 0.732546 mua = 0.857625 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 20 cel = 0.00526086 cea = 0.0032013 cpa = 9.13167e-11 mul = 0.735418 mua = 0.86087 colony = 2 square root of residual = 6.8257 Initial parameter values: b = 20.5 cel = 0.00526086 cea = 0.0032013 cpa = 9.13167e-11 mul = 0.735418 mua = 0.86087 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 20.5 cel = 0.0051902 cea = 0.00334482 cpa = 1.36899e-10 mul = 0.738204 mua = 0.864038 colony = 2 square root of residual = 6.8759 Initial parameter values: b = 21 cel = 0.0051902 cea = 0.00334482 cpa = 1.36899e-10 mul = 0.738204 mua = 0.864038 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 21 cel = 0.00512401 cea = 0.00349186 cpa = 3.97944e-13 mul = 0.740856 mua = 0.867196 colony = 2 square root of residual = 6.92512 Initial parameter values: b = 21.5 cel = 0.00512401 cea = 0.00349186 cpa = 3.97944e-13 mul = 0.740856 mua = 0.867196 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 21.5 cel = 0.00505894 cea = 0.00363296 cpa = 4.03178e-13 mul = 0.743454 mua = 0.870261 colony = 2 square root of residual = 6.9734 Initial parameter values: b = 22 cel = 0.00505894 cea = 0.00363296 cpa = 4.03178e-13 mul = 0.743454 mua = 0.870261 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 22 cel = 0.00499621 cea = 0.00377254 cpa = 2.35376e-13 mul = 0.74596 mua = 0.873286 colony = 2 square root of residual = 7.02078 Initial parameter values: b = 22.5 cel = 0.00499621 cea = 0.00377254 cpa = 2.35376e-13 mul = 0.74596 mua = 0.873286 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 22.5 cel = 0.00493606 cea = 0.00391161 cpa = 2.22063e-13 mul = 0.748387 mua = 0.87626 colony = 2 square root of residual = 7.06728 Initial parameter values: b = 23 cel = 0.00493606 cea = 0.00391161 cpa = 2.22063e-13 mul = 0.748387 mua = 0.87626 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 23 cel = 0.00487821 cea = 0.00404999 cpa = 2.09786e-13 mul = 0.750732 mua = 0.879195 colony = 2 square root of residual = 7.11293 Initial parameter values: b = 23.5 cel = 0.00487821 cea = 0.00404999 cpa = 2.09786e-13 mul = 0.750732 mua = 0.879195 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 23.5 cel = 0.00482213 cea = 0.00418597 cpa = 1.90198e-13 mul = 0.753039 mua = 0.882032 colony = 2 square root of residual = 7.15777 Initial parameter values: b = 24 cel = 0.00482213 cea = 0.00418597 cpa = 1.90198e-13 mul = 0.753039 mua = 0.882032 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 24 cel = 0.00476956 cea = 0.00432802 cpa = 1.39447e-11 mul = 0.75515 mua = 0.885029 colony = 2 square root of residual = 7.20182 Initial parameter values: b = 24.5 cel = 0.00476956 cea = 0.00432802 cpa = 1.39447e-11 mul = 0.75515 mua = 0.885029 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 24.5 cel = 0.00471822 cea = 0.00446644 cpa = 2.6123e-13 mul = 0.757267 mua = 0.887869 colony = 2 square root of residual = 7.24509 Initial parameter values: b = 25 cel = 0.00471822 cea = 0.00446644 cpa = 2.6123e-13 mul = 0.757267 mua = 0.887869 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 25 cel = 0.00466971 cea = 0.00460829 cpa = 1.55073e-11 mul = 0.759286 mua = 0.890727 colony = 2 square root of residual = 7.28762 Initial parameter values: b = 25.5 cel = 0.00466971 cea = 0.00460829 cpa = 1.55073e-11 mul = 0.759286 mua = 0.890727 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 25.5 cel = 0.00462277 cea = 0.00475049 cpa = 1.56618e-11 mul = 0.761209 mua = 0.893611 colony = 2 square root of residual = 7.3294 Initial parameter values: b = 26 cel = 0.00462277 cea = 0.00475049 cpa = 1.56618e-11 mul = 0.761209 mua = 0.893611 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 26 cel = 0.00457759 cea = 0.00489333 cpa = 1.58615e-11 mul = 0.763076 mua = 0.896457 colony = 2 square root of residual = 7.37047 Initial parameter values: b = 26.5 cel = 0.00457759 cea = 0.00489333 cpa = 1.58615e-11 mul = 0.763076 mua = 0.896457 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 26.5 cel = 0.00453573 cea = 0.00504442 cpa = 3.31891e-10 mul = 0.764788 mua = 0.899434 colony = 2 square root of residual = 7.41084 Initial parameter values: b = 27 cel = 0.00453573 cea = 0.00504442 cpa = 3.31891e-10 mul = 0.764788 mua = 0.899434 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 27 cel = 0.00449398 cea = 0.00519023 cpa = 1.51962e-11 mul = 0.766525 mua = 0.902259 colony = 2 square root of residual = 7.45051 Initial parameter values: b = 27.5 cel = 0.00449398 cea = 0.00519023 cpa = 1.51962e-11 mul = 0.766525 mua = 0.902259 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 27.5 cel = 0.0044551 cea = 0.00534344 cpa = 3.13319e-10 mul = 0.768149 mua = 0.905149 colony = 2 square root of residual = 7.48951 Initial parameter values: b = 28 cel = 0.0044551 cea = 0.00534344 cpa = 3.13319e-10 mul = 0.768149 mua = 0.905149 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 28 cel = 0.0044173 cea = 0.00549718 cpa = 3.20303e-10 mul = 0.769717 mua = 0.908031 colony = 2 square root of residual = 7.52783 Initial parameter values: b = 28.5 cel = 0.0044173 cea = 0.00549718 cpa = 3.20303e-10 mul = 0.769717 mua = 0.908031 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 28.5 cel = 0.00438101 cea = 0.00565363 cpa = 3.28354e-10 mul = 0.771224 mua = 0.910911 colony = 2 square root of residual = 7.56549 Initial parameter values: b = 29 cel = 0.00438101 cea = 0.00565363 cpa = 3.28354e-10 mul = 0.771224 mua = 0.910911 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 29 cel = 0.00434618 cea = 0.00581292 cpa = 3.69521e-10 mul = 0.772663 mua = 0.913812 colony = 2 square root of residual = 7.6025 Initial parameter values: b = 29.5 cel = 0.00434618 cea = 0.00581292 cpa = 3.69521e-10 mul = 0.772663 mua = 0.913812 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 29.5 cel = 0.00431295 cea = 0.00597643 cpa = 3.82439e-10 mul = 0.774048 mua = 0.916698 colony = 2 square root of residual = 7.63886 Initial parameter values: b = 30 cel = 0.00431295 cea = 0.00597643 cpa = 3.82439e-10 mul = 0.774048 mua = 0.916698 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 30 cel = 0.00428134 cea = 0.00614494 cpa = 3.82666e-10 mul = 0.775366 mua = 0.9196 colony = 2 square root of residual = 7.67457 Initial parameter values: b = 30.5 cel = 0.00428134 cea = 0.00614494 cpa = 3.82666e-10 mul = 0.775366 mua = 0.9196 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 30.5 cel = 0.00425115 cea = 0.00631825 cpa = 3.85782e-10 mul = 0.776622 mua = 0.922508 colony = 2 square root of residual = 7.70965 Initial parameter values: b = 31 cel = 0.00425115 cea = 0.00631825 cpa = 3.85782e-10 mul = 0.776622 mua = 0.922508 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 31 cel = 0.00422334 cea = 0.00650213 cpa = 1.91148e-11 mul = 0.777788 mua = 0.925477 colony = 2 square root of residual = 7.74408 Initial parameter values: b = 31.5 cel = 0.00422334 cea = 0.00650213 cpa = 1.91148e-11 mul = 0.777788 mua = 0.925477 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 31.5 cel = 0.00419788 cea = 0.0066967 cpa = 9.3586e-13 mul = 0.778842 mua = 0.928563 colony = 2 square root of residual = 7.77788 Initial parameter values: b = 32 cel = 0.00419788 cea = 0.0066967 cpa = 9.3586e-13 mul = 0.778842 mua = 0.928563 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 32 cel = 0.00417379 cea = 0.00689871 cpa = 4.43905e-14 mul = 0.77984 mua = 0.931639 colony = 2 square root of residual = 7.81105 Initial parameter values: b = 32.5 cel = 0.00417379 cea = 0.00689871 cpa = 4.43905e-14 mul = 0.77984 mua = 0.931639 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 32.5 cel = 0.00415142 cea = 0.00710962 cpa = 2.42409e-14 mul = 0.78077 mua = 0.934755 colony = 2 square root of residual = 7.84358 Initial parameter values: b = 33 cel = 0.00415142 cea = 0.00710962 cpa = 2.42409e-14 mul = 0.78077 mua = 0.934755 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 33 cel = 0.00413003 cea = 0.00732582 cpa = 2.35896e-14 mul = 0.781668 mua = 0.93784 colony = 2 square root of residual = 7.87549 Initial parameter values: b = 33.5 cel = 0.00413003 cea = 0.00732582 cpa = 2.35896e-14 mul = 0.781668 mua = 0.93784 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 33.5 cel = 0.00410961 cea = 0.0075492 cpa = 2.34461e-14 mul = 0.782515 mua = 0.940902 colony = 2 square root of residual = 7.90678 Initial parameter values: b = 34 cel = 0.00410961 cea = 0.0075492 cpa = 2.34461e-14 mul = 0.782515 mua = 0.940902 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 34 cel = 0.00408928 cea = 0.00777205 cpa = 2.34176e-14 mul = 0.783378 mua = 0.943858 colony = 2 square root of residual = 7.93746 Initial parameter values: b = 34.5 cel = 0.00408928 cea = 0.00777205 cpa = 2.34176e-14 mul = 0.783378 mua = 0.943858 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 34.5 cel = 0.00406946 cea = 0.00799783 cpa = 2.33937e-14 mul = 0.784229 mua = 0.946743 colony = 2 square root of residual = 7.96754 Initial parameter values: b = 35 cel = 0.00406946 cea = 0.00799783 cpa = 2.33937e-14 mul = 0.784229 mua = 0.946743 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 35 cel = 0.00405127 cea = 0.00823419 cpa = 2.33732e-14 mul = 0.785027 mua = 0.949643 colony = 2 square root of residual = 7.99703 Initial parameter values: b = 35.5 cel = 0.00405127 cea = 0.00823419 cpa = 2.33732e-14 mul = 0.785027 mua = 0.949643 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 35.5 cel = 0.00403237 cea = 0.00846501 cpa = 2.33752e-14 mul = 0.785869 mua = 0.952375 colony = 2 square root of residual = 8.02596 Initial parameter values: b = 36 cel = 0.00403237 cea = 0.00846501 cpa = 2.33752e-14 mul = 0.785869 mua = 0.952375 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 36 cel = 0.00401349 cea = 0.00869494 cpa = 2.33746e-14 mul = 0.786727 mua = 0.954999 colony = 2 square root of residual = 8.05433 Initial parameter values: b = 36.5 cel = 0.00401349 cea = 0.00869494 cpa = 2.33746e-14 mul = 0.786727 mua = 0.954999 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 36.5 cel = 0.00399355 cea = 0.00891507 cpa = 2.3392e-14 mul = 0.787648 mua = 0.957428 colony = 2 square root of residual = 8.08218 Initial parameter values: b = 37 cel = 0.00399355 cea = 0.00891507 cpa = 2.3392e-14 mul = 0.787648 mua = 0.957428 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 37 cel = 0.00397455 cea = 0.00914004 cpa = 2.3384e-14 mul = 0.788554 mua = 0.95982 colony = 2 square root of residual = 8.10952 Initial parameter values: b = 37.5 cel = 0.00397455 cea = 0.00914004 cpa = 2.3384e-14 mul = 0.788554 mua = 0.95982 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 37.5 cel = 0.00395558 cea = 0.00936284 cpa = 2.33766e-14 mul = 0.789481 mua = 0.962107 colony = 2 square root of residual = 8.13637 Initial parameter values: b = 38 cel = 0.00395558 cea = 0.00936284 cpa = 2.33766e-14 mul = 0.789481 mua = 0.962107 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 38 cel = 0.00393651 cea = 0.00958227 cpa = 2.33731e-14 mul = 0.790432 mua = 0.964285 colony = 2 square root of residual = 8.16275 Initial parameter values: b = 38.5 cel = 0.00393651 cea = 0.00958227 cpa = 2.33731e-14 mul = 0.790432 mua = 0.964285 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 38.5 cel = 0.00391772 cea = 0.00980073 cpa = 2.33731e-14 mul = 0.791398 mua = 0.966381 colony = 2 square root of residual = 8.18868 Initial parameter values: b = 39 cel = 0.00391772 cea = 0.00980073 cpa = 2.33731e-14 mul = 0.791398 mua = 0.966381 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 39 cel = 0.0038981 cea = 0.0100105 cpa = 2.33731e-14 mul = 0.792396 mua = 0.968335 colony = 2 square root of residual = 8.21418 Initial parameter values: b = 39.5 cel = 0.0038981 cea = 0.0100105 cpa = 2.33731e-14 mul = 0.792396 mua = 0.968335 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 39.5 cel = 0.00387771 cea = 0.0102126 cpa = 2.33731e-14 mul = 0.793412 mua = 0.970164 colony = 2 square root of residual = 8.23927 Initial parameter values: b = 40 cel = 0.00387771 cea = 0.0102126 cpa = 2.33731e-14 mul = 0.793412 mua = 0.970164 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 40 cel = 0.00385853 cea = 0.0104177 cpa = 2.33731e-14 mul = 0.794444 mua = 0.971959 colony = 2 square root of residual = 8.26396 Initial parameter values: b = 40.5 cel = 0.00385853 cea = 0.0104177 cpa = 2.33731e-14 mul = 0.794444 mua = 0.971959 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 40.5 cel = 0.0038397 cea = 0.01062 cpa = 2.33731e-14 mul = 0.795506 mua = 0.973675 colony = 2 square root of residual = 8.28827 Initial parameter values: b = 41 cel = 0.0038397 cea = 0.01062 cpa = 2.33731e-14 mul = 0.795506 mua = 0.973675 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 41 cel = 0.00381952 cea = 0.0108102 cpa = 2.33731e-14 mul = 0.796582 mua = 0.97525 colony = 2 square root of residual = 8.31223 Initial parameter values: b = 41.5 cel = 0.00381952 cea = 0.0108102 cpa = 2.33731e-14 mul = 0.796582 mua = 0.97525 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 41.5 cel = 0.00380036 cea = 0.0110043 cpa = 2.33731e-14 mul = 0.79765 mua = 0.976812 colony = 2 square root of residual = 8.33583 Initial parameter values: b = 42 cel = 0.00380036 cea = 0.0110043 cpa = 2.33731e-14 mul = 0.79765 mua = 0.976812 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 42 cel = 0.00378152 cea = 0.0111963 cpa = 2.33731e-14 mul = 0.798729 mua = 0.978314 colony = 2 square root of residual = 8.35909 Initial parameter values: b = 42.5 cel = 0.00378152 cea = 0.0111963 cpa = 2.33731e-14 mul = 0.798729 mua = 0.978314 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 42.5 cel = 0.00376294 cea = 0.0113852 cpa = 2.33731e-14 mul = 0.799826 mua = 0.979752 colony = 2 square root of residual = 8.38202 Initial parameter values: b = 43 cel = 0.00376294 cea = 0.0113852 cpa = 2.33731e-14 mul = 0.799826 mua = 0.979752 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 43 cel = 0.00374394 cea = 0.0115608 cpa = 2.33731e-14 mul = 0.801003 mua = 0.98104 colony = 2 square root of residual = 8.40464 Initial parameter values: b = 43.5 cel = 0.00374394 cea = 0.0115608 cpa = 2.33731e-14 mul = 0.801003 mua = 0.98104 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 43.5 cel = 0.0037266 cea = 0.0117502 cpa = 2.33731e-14 mul = 0.802111 mua = 0.982406 colony = 2 square root of residual = 8.42696 Initial parameter values: b = 44 cel = 0.0037266 cea = 0.0117502 cpa = 2.33731e-14 mul = 0.802111 mua = 0.982406 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 44 cel = 0.00371022 cea = 0.0119345 cpa = 3.23623e-14 mul = 0.803162 mua = 0.983689 colony = 2 square root of residual = 8.44898 Initial parameter values: b = 44.5 cel = 0.00371022 cea = 0.0119345 cpa = 3.23623e-14 mul = 0.803162 mua = 0.983689 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 44.5 cel = 0.00369244 cea = 0.0121135 cpa = 3.25375e-14 mul = 0.804268 mua = 0.984925 colony = 2 square root of residual = 8.47071 Initial parameter values: b = 45 cel = 0.00369244 cea = 0.0121135 cpa = 3.25375e-14 mul = 0.804268 mua = 0.984925 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 45 cel = 0.0036717 cea = 0.0122826 cpa = 2.33733e-14 mul = 0.805361 mua = 0.986109 colony = 2 square root of residual = 8.49218 Initial parameter values: b = 45.5 cel = 0.0036717 cea = 0.0122826 cpa = 2.33733e-14 mul = 0.805361 mua = 0.986109 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 45.5 cel = 0.00365745 cea = 0.0124652 cpa = 3.26397e-14 mul = 0.806489 mua = 0.987275 colony = 2 square root of residual = 8.51337 Initial parameter values: b = 46 cel = 0.00365745 cea = 0.0124652 cpa = 3.26397e-14 mul = 0.806489 mua = 0.987275 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 46 cel = 0.00363821 cea = 0.0126284 cpa = 2.33731e-14 mul = 0.807659 mua = 0.98834 colony = 2 square root of residual = 8.5343 Initial parameter values: b = 46.5 cel = 0.00363821 cea = 0.0126284 cpa = 2.33731e-14 mul = 0.807659 mua = 0.98834 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 46.5 cel = 0.00363187 cea = 0.0127191 cpa = 2.92829e-14 mul = 0.808515 mua = 0.988593 colony = 2 square root of residual = 8.55528 Initial parameter values: b = 47 cel = 0.00363187 cea = 0.0127191 cpa = 2.92829e-14 mul = 0.808515 mua = 0.988593 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 47 cel = 0.00360601 cea = 0.0129763 cpa = 2.33731e-14 mul = 0.809873 mua = 0.990518 colony = 2 square root of residual = 8.57541 Initial parameter values: b = 47.5 cel = 0.00360601 cea = 0.0129763 cpa = 2.33731e-14 mul = 0.809873 mua = 0.990518 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 47.5 cel = 0.00358991 cea = 0.0131488 cpa = 2.33732e-14 mul = 0.810931 mua = 0.991581 colony = 2 square root of residual = 8.59559 Initial parameter values: b = 48 cel = 0.00358991 cea = 0.0131488 cpa = 2.33732e-14 mul = 0.810931 mua = 0.991581 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 48 cel = 0.00357391 cea = 0.0133125 cpa = 2.33741e-14 mul = 0.812037 mua = 0.992553 colony = 2 square root of residual = 8.61555 Initial parameter values: b = 48.5 cel = 0.00357391 cea = 0.0133125 cpa = 2.33741e-14 mul = 0.812037 mua = 0.992553 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 48.5 cel = 0.00355884 cea = 0.013485 cpa = 2.33739e-14 mul = 0.813099 mua = 0.993578 colony = 2 square root of residual = 8.63527 Initial parameter values: b = 49 cel = 0.00355884 cea = 0.013485 cpa = 2.33739e-14 mul = 0.813099 mua = 0.993578 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 49 cel = 0.00354471 cea = 0.0136606 cpa = 2.33826e-14 mul = 0.814169 mua = 0.9946 colony = 2 square root of residual = 8.65477 Initial parameter values: b = 49.5 cel = 0.00354471 cea = 0.0136606 cpa = 2.33826e-14 mul = 0.814169 mua = 0.9946 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 49.5 cel = 0.00355004 cea = 0.0138095 cpa = 3.56448e-07 mul = 0.814523 mua = 0.995137 colony = 2 square root of residual = 8.67451 Initial parameter values: b = 50 cel = 0.00355004 cea = 0.0138095 cpa = 3.56448e-07 mul = 0.814523 mua = 0.995137 colony = 2 Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the value of the function tolerance. Estimated parameter values: b = 50 cel = 0.00351648 cea = 0.0139993 cpa = 4.37469e-14 mul = 0.816322 mua = 0.996514 colony = 2 square root of residual = 8.69311 Graphics saved as "beetle_homotopy.png" beetle_dynamic_test(): Normal end of execution. 04-Nov-2026 19:05:06