The three models share the same R0

syms beta_m beta_h theta_m theta_h mu_m gamma_h mu_h R_0
pars=[beta_m beta_h theta_m theta_h mu_m gamma_h mu_h];
R_0=sqrt(beta_m*beta_h*theta_h*theta_m/((theta_m+mu_m)*(gamma_h+mu_h)*(theta_h+mu_h)*(mu_m)))
R_0 = 
[T,~]=gsua_dpmat(R_0,pars,1,'R_0')
Introduce ranges in the following order:
ans = 
T = 7×2 table
 RangeNominal
1 beta_m000
2 beta_h000
3 theta_m000
4 theta_h000
5 mu_m000
6 gamma_h000
7 mu_h000
load Results7.mat
common=intersect(T7.Properties.RowNames,T.Properties.RowNames)
common = 6×1 cell array
{'beta_h' } {'beta_m' } {'gamma_h'} {'mu_m' } {'theta_h'} {'theta_m'}
T{'mu_h','Range'}=[4e-4,4e-4];
T{'mu_h','Nominal'}=4e-4;
T{common,'Range'}=T7{common,'Range'};
%T.Nominal=[];
T{common,'T7'}=T7{common,'Estlsqc'}(:,1:476);
Warning: The new variables being added to the table have fewer rows than the table. They have been extended with rows containing default values.
load Results8_complete.mat
T{common,'T8'}=T8{common,'Estlsqc'}(:,1:158);
Warning: The new variables being added to the table have fewer rows than the table. They have been extended with rows containing default values.
load Results10_complete.mat
T{common,'T10'}=T10{common,'Estlsqc'}(:,1:136);
Warning: The new variables being added to the table have fewer rows than the table. They have been extended with rows containing default values.
T.T7(end,:)=repmat(4e-4,1,476);
T.T8(end,:)=repmat(4e-4,1,158);
T.T10(end,:)=repmat(4e-4,1,136);
clf
gsua_eval(T.T7,T);
clf
gsua_eval(T.T8,T);
clf
gsua_eval(T.T10,T);