rom[0, 0] = h * (fun(a) + fun(b)) / 2