% Logistic parameters

clear all, close all
month=[0:1:11]';
cases=[494 759 1016 1215 1619 2964 8489 22377 32618 38153 41462 42375]';
b0 = [1; 1];
global ParamsWeightLogistic
ParamsWeightLogistic = [1e4; 1e-5];
[b_log,r_log,J_log,C_log,mse_log] = nlinfit(month(2:12),cases(2:12),@logistic,b0);
b_log_true = b_log.*ParamsWeightLogistic