function pemodelan_mikrobia_coli_2_MSchnute_reparameter
clear
clc

global umax Nmax q0 N0 logNcoba Nmin tspan Nccoba lnNccoba lag
umax=0.695709639;   %jam^-1
lag=2;
Nmax=10^7.98;      %CFU/ml
N0=10^2.85;       %CFU/ml
Nmin=(1-(1e-6))*10^2.85;       %CFU/ml
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];
lnNccoba=log(Nccoba)
tspan=[0;2;4;6;8;12;16;20;24]; 


Atebak=[0.1;0.1];

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

end

function RMSE_C=funcRMSEC(g)
global umax Nmax q0 N0 logNcoba Nmin tspan Nccoba lncNcoba lag
logNccoba_hit=log(N0)+((umax.*((1-g(2))./g(1))).*((1-(g(2).*exp((g(1).*lag)+1-g(2)-(g(1).*tspan))))./(1-g(2))).^(1./g(2)));
Nccoba_hit=exp(logNccoba_hit);
log10Nccoba_hit=log10(Nccoba_hit)

%mencari R2
SStot=sum((log10(Nccoba)-mean(log10(Nccoba))).^2);
SSres=sum((log10(Nccoba)-log10Nccoba_hit).^2);
R2=1-(SSres./SStot)

RMSE_C=(sum((log10(Nccoba)-log10Nccoba_hit).^2)./length(tspan)).^0.5
plot(tspan,log10(Nccoba),'bo',tspan,log10Nccoba_hit,'r-')
end