Approximate solutions to fractional subdiffusion equations