The entire trial was performed in triplicate. For each formulation (Control, SUB1, SUB2 and SUB3), results were expressed as mean and standard error of the results obtained from the three independent batches. The difference between mean of values was determined using one-way analysis of variance (ANOVA) and multiple comparisons of means were done using Bonferroni Post Hoc procedure to evaluate the statistical significance (P ≤ 0.05) among formulations. Formulation was assigned as fixed factor. Statistical analysis of the sensory test was done according to the instructions of UNE-EN ISO 4120: 2004. The statistical analysis of data was done using the STATA/IC 12.1 program