Use Liebeman method $($also known as Gauss Siedel method$)$ with over relaxation to determine the flow in the aquifer shown in figure $1,$ with the following properties: hydraulic conductivity $k=10\frac{m}{?}$ day, groundwater recharge $N=0m\frac{m}{?}$ day, thickness of the aquifer $D=50m.$

The boundary conditions of the aquifer are

BOUNDARY $1$: no$-$flow near the mountains;

BOUNDARY $2$: linear from $100$ m $($near mountains$)$ to $0$ m $($near sea$)$;

BOUNDARY $3$: constant seawater level of $0$ m ;

BOUNDARY $4$: no$-$flow.

Figure $1.$ Aquifer with boundary conditions

Run the model for two possible solutions:

$\backslash$table$[[$Solution I:,Solution II:$],[$Initial estimate head $=50m,-$ Initial estimate head$=50$m$],[$Convergence$-$criterion$=1$m$,-$ Convergence$-$criterion$=0\mathrm{.}01$ m$],[$Number of iterations $=13,-$ Number of iterations$=74],[$Head $1=29,17m,-$ Head $1=20,67m$