Double integral-balance method to the fractional subdiffusion equation: Approximate solutions, optimization problems to be resolved and numerical simulations