call relaxation on fine grid and compute the numerical solution
gauss_seidel_mg(nxf, nyf, dx, dy, f, un, v1)