You would not believe how long I looked for the problem and never found it. Thanks a bunch.
This is the output I get as well. According to the published paper which I will attach, there is a second part for 3 iterations.
Problem solved!
You have used Curly brackets instead of Round brackets / parantheses in line
H[x_, t_] = (1 - t)*(g[x]) + t (f[x]);
An updated file is attached. Have fun!
Is this the result you expect?
Attachments:
