highly significantp < 0.0001
The quadratic equation generated from the models ( Figure 2 ) is the following Equation (2): (2) Rupture strength = −71.72 + 19.14A + 15.99B + 3381.93C + 4.78D − 0.43AB + 136.21AC − 0.29AD − 460.12BC + 0.88BD − 79.31CD − 4.30A2 − 32560.72C2 − 0.22D2 − 3.43B2 The compressive strength response surface model was highly significant (F = 22.36, p < 0.0001, Table S3 ), accounting for 91.16% of the total variation ( R 2 = 0.9116, adjusted R 2 = 0.9116, predicted R 2 = 0.834) with a non-significant lack-of-fit ( p = 0.683), confirming the model reliability.