Omega}$ there exists a $z_{y}\in {{\mathbb R}}^{N}\backslash{\Omega}$ such that ${\overline}{B_{\theta}(z_{y})}\cap{\overline}{{\Omega}}= \{ y \}$. The boundary data $f$ and $g$ are assumed to be Lipschitz functions.

Notice that this system involves two differential operators, the usual Laplacian $$\Delta \phi = \sum\limits_{i=1}^{N} \partial_{x_{i}x_{i}}\phi$$ and the infinity Laplacian (see [@Cran]) $$\Delta_{\infty}\phi =  
 \langle D^{2}\phi \frac{\nabla \phi}{|\nabla \phi |},\frac{\nabla \phi}{| \nabla \phi |}\rangle 
 =\frac{1}{| \nabla \phi |^2}\sum\limits_{i,j=1}^{N} \partial_{x_{i}}\phi \partial_{x_{i}x_{j}}\phi
 \partial_{x_{j}}\phi.$$

This system is not variational (there is no associated energy). Therefore, to find solutions one possibility is to use monotonicity methods (Perron’s argument). Here we will look at the system in a different way and to obtain existence of solutions we find an approximation using game theory. This approach not only gives existence of solutions but it also provide us with a description that yield some light on the behaveiour of the solutions. At this point we observe that we will understand solutions to the system in the viscosity sense, this is natural since the infinity Laplacian is not variational (see Section \[sect-prelim\] for the precise definition).

The fundamental works by Doob, Feller, Hunt, Kakutani, Kolmogorov and many others show the deep connection between classical potential theory and probability theory. The main idea that is behind this relation is that harmonic functions and martingales have something in common: the mean value formulas. A well known fact is that $u$ is harmonic, that is $u$ verifies the PDE $\Delta u =0$, if and only if it verifies the mean value property $
u(x) =
\frac{1}{|B_\varepsilon (x) |} \int_{B_\varepsilon (x) } u(y) \, dy.
$ In fact, we can relax this condition by requiring that it holds asymptotically $
u(x) =
\frac{1
Omega}$ there exists a $ z_{y}\in { { \mathbb R}}^{N}\backslash{\Omega}$ such that $ { \overline}{B_{\theta}(z_{y})}\cap{\overline}{{\Omega}}= \ { y \}$. The boundary data $ f$ and $ g$ are assumed to be Lipschitz functions. 

 detect that this arrangement involves two differential operators, the usual Laplacian $ $ \Delta \phi = \sum\limits_{i=1}^{N } \partial_{x_{i}x_{i}}\phi$$ and the eternity Laplacian (see [ @Cran ]) $ $ \Delta_{\infty}\phi =  
  \langle D^{2}\phi \frac{\nabla \phi}{|\nabla \phi |},\frac{\nabla \phi}{| \nabla \phi |}\rangle 
  = \frac{1}{| \nabla \phi |^2}\sum\limits_{i, j=1}^{N } \partial_{x_{i}}\phi \partial_{x_{i}x_{j}}\phi 
  \partial_{x_{j}}\phi.$$ 

 This system is not variational (there is no associated department of energy). consequently, to find solutions one possibility is to use monotonicity methods (Perron ’s argumentation). Here we will look at the system in a unlike way and to obtain existence of solutions we find an approximation using plot theory. This approach not only yield existence of solutions but it besides provide us with a description that concede some light on the behaveiour of the solution. At this point we observe that we will understand solution to the system in the viscosity sense, this is natural since the infinity Laplacian is not variational (see Section \[sect - prelim\ ] for the accurate definition). 

 The cardinal workplace by Doob, Feller, Hunt, Kakutani, Kolmogorov and many others show the deep connection between classical potential theory and probability theory. The main mind that is behind this relation is that harmonic routine and martingales have something in common: the mean value rule. A well known fact is that $ u$ is harmonic, that is $ u$ verifies the PDE $ \Delta u = 0 $, if and only if it verifies the bastardly value property $ 
 u(x) = 
 \frac{1}{|B_\varepsilon (x) | } \int_{B_\varepsilon (x) } u(y) \, dy. 
 $ In fact, we can relax this condition by requiring that it holds asymptotically $ 
 u(x) = 
 \frac{1
Omeha}$ there exists a $z_{y}\in {{\mxthbb R}}^{N}\backslash{\Omega}$ such fhat ${\ovefline}{B_{\theta}(z_{y})}\cap{\overline}{{\Omege}}= \{ y \}$. The boundary data $f$ and $g$ are assuled to bw Likschitz functions.

Iktice tmct thjd syvvem involves twp differendial operators, tfe usual Laplacian $$\Delta \phi = \sum\limies_{i=1}^{N} \pattlal_{x_{i}x_{i}}\phi$$ and the ynfihptn Laplacian (see [@Cran]) $$\Delta_{\infty}\pgi =  
 \lengle D^{2}\phi \frac{\mabla \phi}{|\nabla \phi |},\frac{\nabpa \pji}{| \nabla \phi |}\ranglf 
 =\frac{1}{| \naboa \pru |^2}\sum\limits_{i,g=1}^{N} \partial_{q_{n}}\phi \partiam_{x_{i}x_{j}}\phi
 \partial_{x_{j}}\phi.$$

This systeo is uot variatiinql (hvere is no essociwted energy). Bnerefose, to fond solutions pne powsibility is to use mmnotonicity methodf (Perron’s axgument). Here we will ooik at the sysgwm kn z vifrerent wab and to obfain existebce of solutions we fybd an approxijation ufing game theory. This approach not only gibes existence of solutiins but it also provife us witr a description that yield some light on the behaeeiouc uf uhc rilktions. At this point we observe that we will thdtrsnand solutions to the system ij yre viscosity rense, cgia is natural since the insinitt Laplaciwn ix not variational (see Sectiin \[sect-prelii\] for the precise dzfinition).

The funcamenyal works by Doob, Fellex, Hunt, Kakutani, Kllmogorov xnd many others rhoe dhe deep connection betweeg classicel pocential gheoty and [robabilitj theory. The main idea tjat id tehind thid relation is that harmonic funrvions and marjincalvs have sjmethlng in common: tre mean value yormulas. A weul known fzct is vhat $u$ is hawmonic, that iv $u$ verifies vhe PDE $\Dqlta u =0$, uf and uvly if it verigies the mean value property $
u(x) =
\frac{1}{|B_\yareprjlon (x) |} \int_{B_\varzksioon (x) } u(y) \, dy.
$ Im fxct, wv cen rejdx this condhtiov bh requkring that lt foldx asymptotically $
u(x) =
\xrac{1
Omega}$ there exists a $z_{y}\in {{\mathbb R}}^{N}\backslash{\Omega}$ ${\overline}{B_{\theta}(z_{y})}\cap{\overline}{{\Omega}}= y \}$. boundary data $f$ be functions. Notice that system involves two operators, the usual Laplacian $$\Delta \phi \sum\limits_{i=1}^{N} \partial_{x_{i}x_{i}}\phi$$ and the infinity Laplacian (see [@Cran]) $$\Delta_{\infty}\phi = \langle D^{2}\phi \frac{\nabla \phi |},\frac{\nabla \phi}{| \nabla \phi |}\rangle =\frac{1}{| \nabla \phi |^2}\sum\limits_{i,j=1}^{N} \partial_{x_{i}}\phi \partial_{x_{i}x_{j}}\phi \partial_{x_{j}}\phi.$$ system not (there no associated energy). Therefore, to find solutions one possibility is to use monotonicity methods (Perron’s argument). we will look at the system in a way and to obtain of solutions we find an using theory. This not gives of solutions but also provide us with a description that yield some light on the behaveiour of the solutions. At point we we will solutions the in the viscosity is natural since the infinity Laplacian (see Section \[sect-prelim\] for the precise definition). The works by Feller, Hunt, Kakutani, Kolmogorov and many show the deep connection between classical potential theory probability theory. The main idea that is behind this relation is that harmonic functions and something in common: the value formulas. A known is $u$ harmonic, that $u$ verifies the PDE $\Delta u =0$, if and only if verifies the mean value property $ u(x) = \frac{1}{|B_\varepsilon (x) (x) u(y) \, dy. In fact, we can this by requiring that it $ =
Omega}$ there exists a $z_{y}\in {{\mathBb R}}^{N}\backslAsh{\OmEga}$ SucH tHat ${\oVerlIne}{B_{\theta}(z_{y})}\cap{\OVerlIne}{{\Omega}}= \{ y \}$. The boundary daTa $f$ anD $g$ ARe asSUmEd to bE LipschITz FUNctIoNs.

notIcE ThAt thiS syStem invOlves two diFfeReNtial operatoRS, tHe usual LapLacIan $$\Delta \phi = \sUm\lImits_{i=1}^{n} \pArtIAl_{x_{i}x_{I}}\phI$$ and tHe infiNIty LapLacian (see [@crAN]) $$\Delta_{\INfty}\phi =  
 \LANgLe D^{2}\pHi \frac{\nabla \phi}{|\nabLA \pHI |},\frac{\nabla \phi}{| \nAbla \phI |}\rANgLE 
 =\FraC{1}{| \naBla \phi |^2}\sum\lImIts_{i,j=1}^{n} \Partial_{X_{I}}\pHI \PArtIAl_{x_{i}x_{j}}\phi
 \partIal_{x_{j}}\phi.$$

ThiS SysTem is nOt VarIAtionaL (therE iS No aSsociated enErgy). therefore, To find SOlutionS One possIbilitY is To uSe moNOtOnIciTy MEthODs (perROn’s Argument). heRe We wilL looK AT THe sySteM in a DiffeRent way and to oBtaIn exISteNce of SolutIons We Find aN approXimatIoN using game theorY. ThiS approach Not OnLy gIvEs exiSTence oF soLutIons but It also pROviDe US WItH a description that yIeLD SoMe light oN the beHAvEiOUr of the sOlUtiOns. AT THis poInt wE ObServe thaT we wilL UnDeRstand sOlUtions To The SysTem in THe viScositY sense, thIs is nATural since the iNFinity LaplaciAN iS NOt VAriaTioNal (see SectiOn \[seCT-preLim\] fOR tHe pREcise DefinItIOn).

tHe fundamental works bY DOob, FelLer, HuNt, Kakutani, KolMogorov and MANY others sHow tHE dEEp connection beTween Classical pOTential tHeory And probaBility theORY. The main IdeA thAt iS beHINd This relation iS THat hArMonic fuNctIons and MarTinGalEs hAvE somethinG in commoN: tHe MeAn ValUe forMUlas. A welL kNowN fAct Is thaT $U$ is harMonic, That Is $U$ vERifIes the Pde $\DELTa u =0$, iF aNd Only If iT vErifiEs thE MeaN value pRoperty $
u(x) =
\FraC{1}{|b_\varEpSiLon (x) |} \int_{b_\varepsilon (x) } u(Y) \, dY.
$ In fact, we cAn RelAx this CONdition bY requiring that it holds asYMptoticAllY $
u(x) =
\frAc{1
Omega}$ there exists a $z_ {y}\in {{\ mathb b R }}^ {N }\ba cksl ash{\Omega}$ s u ch t hat ${\overline}{B_{\t heta} (z _ {y}) } \c ap{\o verline } {{ \ O meg a} }=  \{  y \} $. Th e b oundary  data $f$and  $ g$ are assum e dto be Lips chi tz functions .
Notice  t hat thissys tem i nvolve s  two d ifferenti al operat o rs, the u su al L aplacian $$\Delta \p h i = \sum\limit s_{i=1 }^ { N} \ par tia l_{x_{i}x_ {i }}\ph i $$ andt he i n fin i ty Laplacian(see [@Cran ] ) $ $\Delt a_ {\i n fty}\p hi = \la ngle D^{2}\ phi\frac{\na bla \p h i}{|\na b la \phi  |},\f rac {\n abla \p hi }{|  \ n abl a  \ phi |}\ rangle  = \f rac{1 }{|\ n a b la \ phi  |^2 }\sum \limits_{i,j= 1}^ {N}\ par tial_ {x_{i }}\p hi  \par tial_{ x_{i} x_ {j}}\phi
 \part ial_ {x_{j}}\p hi. $$ 

T hi s sys t em isnot  va riation al (the r e i sn o as sociated energy).Th e r ef ore, tofind s o lu ti o ns one p os sib ilit y is to  use mo notonici ty met h od s(Perron ’s  argum en t).  He re we will  lookat the s ystem in a different way and to ob t ai n ex i sten ceof solution s we find  ana pp rox i matio n usi ng ga m e theory. This appr oa ch not  only  gives existe nce of sol u t i ons butit a l so provide us wit h a d escription that yie ld so me light  on the b e h aveiourofthe  so lut i o ns . At this poi n t  weob serve t hat  we wil l u nde rst and  s olutionsto the s ys te min  th e vis c osity se ns e,th isis na t ural s incethein fi n ity  Laplac i an i s no tva riat ion al  (see  Sec t ion  \[sect -prelim\]  fo r  the  p re cise de finition).

T he  fundament al  wo rks by D oob, Fel ler, Hunt, Kakutani, Ko l mogorov  an d man y ot hers show  th e deep  co n nectio n betw een c la ssi c a l pot e n ti alth eory and p r o bab ility  t heor y. Themain idea that isb ehi nd this relat ion  ist h at  ha r mo n icfu n cti o n s and martingal es have so me t hi ng in comm o n:th e meanvalue f ormul a s. A we ll knownfact is t ha t $u $ isharmonic,that is$u$ verif i es th e  P DE $\ Del ta u = 0$ , i f and  onlyi f i t ver ifiesth e mean  valu eproperty  $
u(x) =
\frac{1}{|B_\ vareps ilon(x)  |} \int_ {B_ \ var epsilon ( x) }  u(y) \, d y.$ I n fac t,w e can  rel a xthi s  cond itio n  by requi r in g t h a tit holds as y m p tot icall y $ u(x) = 
\fr ac{1
Omega}$ there_exists a_$z_{y}\in {{\mathbb R}}^{N}\backslash{\Omega}$ such_that ${\overline}{B_{\theta}(z_{y})}\cap{\overline}{{\Omega}}=_\{_y \}$._The_boundary data $f$_and $g$ are_assumed to be Lipschitz_functions.

Notice that this_system_involves two differential operators, the usual Laplacian $$\Delta \phi = \sum\limits_{i=1}^{N} \partial_{x_{i}x_{i}}\phi$$ and the_infinity_Laplacian (see_[@Cran])_$$\Delta_{\infty}\phi_=  
 \langle D^{2}\phi_\frac{\nabla \phi}{|\nabla \phi |},\frac{\nabla \phi}{|_\nabla \phi_|}\rangle 
 =\frac{1}{| \nabla \phi |^2}\sum\limits_{i,j=1}^{N} \partial_{x_{i}}\phi \partial_{x_{i}x_{j}}\phi
_\partial_{x_{j}}\phi.$$

This_system is not_variational (there is no associated energy). Therefore, to find_solutions one possibility is to use_monotonicity methods (Perron’s_argument)._Here_we will look at_the system in a different way_and to obtain existence of solutions_we find an approximation using game theory._This approach not only gives existence_of solutions but it also_provide us_with a description that yield_some light on_the behaveiour_of the solutions._At this point we observe that_we will understand_solutions to the system in the_viscosity_sense, this is_natural_since_the infinity_Laplacian is not_variational_(see Section_\[sect-prelim\]_for the precise definition).

The fundamental works_by_Doob, Feller, Hunt, Kakutani, Kolmogorov and many_others show the deep_connection_between classical potential theory_and probability theory. The main_idea that is behind this relation_is that_harmonic functions_and martingales have something in common: the mean value formulas. A_well known fact is that $u$_is harmonic, that is_$u$ verifies_the_PDE $\Delta u_=0$,_if and_only if it verifies the mean value_property $
u(x)_=
\frac{1}{|B_\varepsilon (x) |} \int_{B_\varepsilon (x) }_u(y) \, dy.
$ In_fact,_we can relax this condition by_requiring that it holds asymptotically $
u(x)_=
\frac{1