function pemodelan_mikrobia_BAL_MGompertz_Zwie
clear
clc

global umax Nmax q0 N0 logNbcoba Nmin Nbcoba tspan lag
umax=1.11675377;   %jam^-1
lag=2;
Nmax=10^9.14;      %CFU/ml
N0=10^1.45;       %CFU/ml
Nbcoba=[10^1.45;10^1.60;10^2.52;10^4.10;10^4.94;10^7.22;10^8.39;10^9.10;10^9.14];
Nccoba=[10^2.85;10^3.5;10^3.73;10^4.01;10^4.81;10^6.21;10^7.73;10^7.79;10^7.98];
logNbcoba=log10(Nbcoba)
tspan=[0;2;4;6;8;12;16;20;24]; 


Atebak=1;

[solution,sse,exitflag,output]=fminsearch(@funcRMSEC,Atebak)

end

function RMSE_C=funcRMSEC(A)
global umax Nmax q0 N0 logNbcoba Nmin tspan Nbcoba lnbNcoba lag
logNbcoba_hit=log(N0)+A.*exp(-exp((umax.*exp(1).*(lag-tspan)./A)+1));
Nbcoba_hit=exp(logNbcoba_hit);
log10Nbcoba_hit=log10(Nbcoba_hit)
RMSE_C=(sum((log10(Nbcoba)-log10Nbcoba_hit).^2)./length(tspan)).^0.5
plot(tspan,log10(Nbcoba),'bo',tspan,log10Nbcoba_hit,'r-')

%mencari R2
SStot=sum((logNbcoba-mean(logNbcoba)).^2);
SSres=sum((logNbcoba-log10Nbcoba_hit).^2);
R2=1-(SSres./SStot)
end