function pemodelan_mikrobia_coli_2_MGompertz_3_pakeRMSE_umax2
clear
clc

global umax Nmax q0 N0 logNcoba Nmin tspan
umax=0.695709639;   %jam^-1
Nmax=10^7.98;      %CFU/ml
N0=10^2.85;       %CFU/ml
Nmin=(1-(1e-6))*10^2.85;       %CFU/ml
logNcoba=[2.85;3.5;3.73;4.01;4.81;6.21;7.73;7.79;7.98];
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 umax Nmax q0 N0 logNcoba tspan
W=y(1)+y(2).*exp(-exp(-y(3).*(tspan-y(4))))

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

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

     % huruf dalam tanda petik pertama menyatakan warna kurva
 end