Given a square scalar linear system Ax = b implement verified_linear_solver(A, b) (feel free to suggest other names) which gives a rigorous bound of the solution x. Use epsilon inflation method described in section 5.5 of horacek paper, which is algorithm 10.7 in Rump's paper Verification methods: Rigorous results using floating-point arithmetic