function pemodelan_mikrobia_BAL_MGompertz_pakeRMSE
clear
clc

global umax Nmax q0 N0 logNcoba Nmin tspan
umax=1.11675377;   %jam^-1
Nmax=10^9.14;      %CFU/ml
N0=10^2.85;       %CFU/ml
Nmin=(1-(1e-6))*10^2.85;       %CFU/ml
logNcoba=[1.45;1.60;2.52;4.10;4.94;7.22;8.39;9.10;9.14];
tspan=[0;2;4;6;8;12;16;20;24]; 


ytebak=[1 1 1 1];

[solution,sse,exitflag,output]=fmincon(@func_MGom,ytebak,[],[],[],[],[0 0 0 0],[])
end

function RMSE=func_MGom(y)
global  Nmax q0 N0 logNcoba tspan
W=y(1)+y(2).*exp(-exp(-y(3).*(tspan-y(4))))

RMSE=(sum((logNcoba-W).^2)./length(tspan)).^0.5
                  % huruf dalam tanda petik kedua menyatakan betuk kurva
plot(tspan,logNcoba,'bo',tspan,W,'r-')


%mencari R2
SStot=sum((logNcoba-mean(logNcoba)).^2);
SSres=sum((logNcoba-W).^2);
R2=1-(SSres./SStot)

     % huruf dalam tanda petik pertama menyatakan warna kurva
 end