Laplace-ligningen

**Laplace-ligningen** er en partiell differensialligning av andre orden. Den spiller en meget viktig rolle i teoretisk fysikk og matematisk analyse. Ved bruk av laplace-operatoren &nabla;&thinsp;2 som virker på en funksjon &Phi; med *N* variable,  kan den skrives på den kompakte formen

:  \nabla^2\Phi = 0 

Generelt kan det da vises at løsningene av ligningen er harmoniske funksjoner i *N* dimensjoner.

Differensialligningen ble funnet i 1782 av Pierre-Simon Laplace i forbindelse med beregning fra Newtons gravitasjonsteori av gravitasjonspotensialet for en kontinuerlig massefordeling. I 1812 viste Siméon Denis Poisson at ligningen ikke gjelder inni en slik fordeling, men kan den utvides med et nytt ledd på høyre side som er proporsjonalt med massetettheten. Den er da gyldig overalt og kalles for Poisson-ligningen.

Laplace-ligningen opptrer ofte ved beskrivelse av stasjonære prosesser i  fysikk. For eksempel vil temperaturen i et legeme som styres av varmeledning, oppfylle ligningen når den ikke lenger forandrer seg med tiden. Også partikkeltettheten ved diffusjon vil være gitt som løsning av ligningen under stasjonære forhold. Likedan vil potensialene både innen elektrostatikken og magnetostatikken utenfor alle elektriske ladninger og strømmer være bestemt av Laplace-ligningen.

I kompleks analyse vil både den reelle og imaginære delen av en kompleks funksjon oppfylle den to-dimensjonale Laplace-ligningen. De er derfor begge harmoniske funksjoner. Dette er en direkte konsekvens av Cauchy–Riemanns ligninger.

## Opprinnelse

I sine astronomiske arbeid basert på Newtons gravitasjonslov hadde Laplace vist at mange beregninger kunne gjøres enklere ved å innføre et  gravitasjonspotensial skapt av en masse i stedet for gravitasjonskraften som den forårsaker. Betrakter man en kontinuerlig massefordeling med tetthet *&rho;*, er dette potensialet i et punkt **r** gitt ved integralet

:  \Phi(\mathbf{r}) = - \int\!d^3x' {G\rho(\mathbf{r'}) \over |\mathbf{r} - \mathbf{r'}|} 

hvor *G* er gravitasjonskonstanten. Ved bruk av kartesiske koordinater kan avstanden *R* = |**r** - **r**'| skrives som

:  R = \sqrt{(x - x')^2 + (y - y')^2 + (z - z')^2} 

slik at &part;*R*/&part;*x* = (*x* - *x'&thinsp;*)/*R*. Dermed er &part;/&part;*x*(1/*R*) = - (1/*R*3)(*x* - *x'&thinsp;*). En videre derivasjon gir så

:  \left({\partial^2\over \partial x^2}\right){1\over R} = {3\over R^5}(x - x')^2 - {1\over R^3} 

Sammen med de tilsvarende dobbeltderiverte i de to andre retningene blir dermed

:  \left({\partial^2\over \partial x^2} + {\partial^2\over \partial y^2} + {\partial^2\over \partial z^2}\right){1\over R} = {3\over R^5}\Big[(x - x')^2 + (y - y')^2 + (z - z')^2\Big] - {3\over R^3} = 0 

Da disse derivasjonene kan tas inn under integrasjonssymbolet for gravitasjonspotensialet, har man dermed

:  \left({\partial^2\over \partial x^2} + {\partial^2\over \partial y^2} + {\partial^2\over \partial z^2}\right)\Phi(x,y,z) = 0 

Dette er Laplace-ligningen i kartesiske koordinater i *N* = 3 dimensjoner.

### Elektrisk potensial

Da Coulombs lov for elektriske krefter har samme form som Newtons gravitasjonslov, er det å forvente at Laplace-ligningen også vil kunne beskrive det elektriske potensialet &phi;. Det følger fra Maxwells ligninger. Hvis man ikke har noe magnetisk felt, sier Faradays induksjonslov at  det elektriske feltet **E** må oppfylle  . Det må derfor kunne skrives som gradienten av potensialet, det vil si at  Når man i tillegg er i et område hvor det ikke er noen elektriske ladninger, sier Maxwells første ligning at . Det betyr at  som er Laplace-ligningen. Den tillater ofte å finne løsninger for det elektriske potensialet og derav feltet **E**, noe som er lettere enn å beregne det elektriske feltet direkte.

### Potensialstrømning

En ideell væske beskrives ved Euler-ligningene. Under vanlig forhold er væsken **virvelfri**. Det betyr at curl til dens hastighet **v**  er null, det vil si at . På samme måte som for det elektriske feltet, betyr det at hastighetsfeltet til væsken må kunne skrives som gradienten av et potensial &Psi;. Da har man at  og bevegelsen til væsken sies å være en **potensialstrømning**. I stedet for en beskrivelse ved de tre komponentene til  vektorfeltet **v**, kan man da klare seg med den ene potensialfunksjonen &Psi;.

Beregning av dette potensialet forenkles når væsken i tillegg er  **inkompressibel**, det vil si at den medfølgende tettheten er konstant. Kontinuitetsligningen sier da at væskens hastighet **v** tilfredsstiller .  Benytter man her kravet  om virvelfrihet, følger igjen Laplace-ligningen 

Dette er vanligvis en tredimensjonal ligning da potensialet avhenger av tre koordinater, det vil si at &Psi; = &Psi;(*x,y,z*). Men løsning av den todimensjonale ligningen kan også ha praktisk betydning. For eksempel, så følger det fra Young-Laplace-ligningen for en glatt hinne av såpevann med en høyde *H* = *H*(*x,y*) over et plant referansenivå, at den er bestemt av  hvor nå Laplace-operatoren involverer bare de to koordinatene *x* og *y*. Den glatte egenskapen til denne flaten og likedan såpebobler er en konsekvens av overflatespenningen til såpevannet.

## Generell definisjon

Laplace-ligningen i et *N*-dimensjonalt euklidsk rom med kartesiske koordinater **x** = (*x*1, *x*2, ... ,*x**N*&thinsp;) er definert som

:  \left({\partial^2\over \partial x_1^2} + {\partial^2\over \partial x_2^2} + \cdots + {\partial^2\over \partial x_N^2}\right)\Phi(x_1,x_2,\cdots,x_N) = 0 

Ved å innføre nabla-operatoren for å kunne beregne gradient og divergens, kan ligningen på mer kompakt form skrives som

:  \boldsymbol{\nabla}\cdot\boldsymbol{\nabla}\Phi(\mathbf{r}) \equiv \nabla^2\Phi(\mathbf{r}) = 0 

Her inngår på venstre side Laplace-operatoren i *N* dimensjoner. Ved å benytte andre koordinatsystem i dette rommet, vil den ta andre former som kan finnes fra denne generelle definisjonen.

## Kulesymmetrisk løsning

Den enkleste, ikke-trivielle løsningen av Laplace-ligningen finnes ved å anta at den er invariant under rotasjoner. Den er da kun en funksjon av radius

:  r = |\mathbf x| = (x_1^2+x_2^2+ \cdots +x_N^2)^{1/2} 

Potensialet kan da skrives som &Phi;(**r**) = *u*(*r*&thinsp;). Ved å benytte at  &part;*r*/&part;*x<sub>i</sub>* = *x<sub>i</sub>* /*r*,  blir dermed

:  \frac{\partial u}{\partial x_i}= u'(r)\frac{x_i}{r}, \; \; \; \frac{\partial^2 u}{\partial x_i^2} = u*(r)\frac{x_i^2}{r^2} + u'(r) \left( \frac{1}{r} - \frac{x_i^2}{r^3} \right )</math>

for hver koordinat *xi*. Dermed tar ligningen formen

:<math> \nabla^2 u = \sum_{i=1}^N \frac{\partial^2 u}{\partial x_i^2} = u*(r) + \frac{N - 1}{r}u'(r) = 0 

Herav følger at løsningen må oppfylle betingelsen

:  \frac{u*}{u'}= \frac{1 - N}{r} = \log (u')'</math>

som betyr at

: <math>u'(r) = \frac{a}{r^{N - 1}}</math>

der *a* er en konstant.  Dermed vil de kulesymmetriske løsningene være av formen

: <math> u(r) = \begin{cases} a\log r + c, \; \; \; \;   N = 2 \\ a/r^{N - 2} + c, \; \; N > 2 \end{cases}</math>

hvor *c* er en konstant som sammen med *a* kan bestemmes av grensebetingelsene.

===Standardform===
Det er vanlig å skrive disse løsningene på en standardform som benyttes ved beregning av [[Greens funksjon]]. Hvis man innfører den fulle  [[romvinkel]] i [[Kule (geometri)#Kuler i høyere dimensjoner|*N* dimensjoner]],

: <math>\Omega_{N-1} = \frac{2 \pi^{N/2}}{\Gamma(N/2)} </math>

når den uttrykkes ved [[gammafunksjon]]en &Gamma;(*z*&thinsp;),  er de kulesymmetriske løsningene

:<math>\Phi(\mathbf{r}) =
\left\{
\begin{array}{rl}
-\frac {1} {2 \pi} \ln \ , &N = 2\\ 
{1\over (N-2)\,\Omega_{N-1} |\mathbf{r}|^{N-2}}\ , &N > 2 \\
\end{array}
\right.
</math>

De gir potensialet utenfor en [[Kontinuitetsligning#Punktpartikler|punktpartikkel]] i [[origo]]. Med denne standardnormeringen er potensialet &Phi; = 1/4*&pi;&thinsp;r* i *N* = 3 dimensjoner.

[[Elektrisk potensial|Coulomb-potensialet]] oppfyller generelt Laplace-ligningen i alle områder av rommet hvor det ikke er ladninger. Da potensialet skapt av en punktladning alltid er  kulesymmetrisk, betyr det at Coulomb-potensialet utenfor en slik ladning i *N* = 2 dimensjoner varierer logaritmisk med avstanden og går ikke mot null ved store avstander. Derimot vil Coulomb-potensialet skapt av en punktladning i rom med *N* > 3 dimensjoner avta mot null raskere ved store avstander enn i tre dimensjoner.

==To dimensjoner==

Den todimensjonale Laplace-ligningen

: <math> \left({\partial^2\over \partial x^2} + {\partial^2\over \partial y^2}\right)\Phi(x,y) = 0 </math>

opptrer i [[kompleks analyse]] hvor en [[Funksjon (matematikk)|funksjon]] *f*&thinsp;(*x + iy*) = *u*(*x,y*) + *iv*(*x,y*) må oppfylle [[Cauchy-Riemanns ligninger]]

:<math>\frac{\partial u}{\partial x} = \frac {\partial v}{\partial y},\; \; \; \frac{\partial u}{\partial y} = -\frac {\partial v}{\partial x}</math>

Ved å derivere den første med hensyn på *x* og den andre med hensyn på *y*, vil en addisjon av de to resulterende ligningene gi at den reelle funksjonen *u*(*x,y*) tilfredsstiller Laplace-ligningen i to dimensjoner. Det samme gjelder for funksjonen *v*(*x,y*) som vises ved å derivere de to ligningene med hensyn på *y* i stedet.

===Homogene polynom og harmoniske løsninger===
I *N* = 2 dimensjoner er det kun løsningen ln(*x*<sup>2</sup> + *y*<sup>2</sup>)  som er rotasjonssymmetrisk. Men oppgis dette kravet, er det en uendelighet av andre løsninger. En spesiell klasse  er [[polynom|homogene polynom]]. De to enkleste er *xy* og *x*<sup>2</sup> - *y*<sup>2</sup> som er slike polynom av andre grad og oppfyller Laplace-ligningen. Uttrykkes disse i [[polarkoordinatsystem|polare koordinater]]  og , er de to løsningene proporsjonale med henholdsvis  og .

Denne sammenhengen kan man også se ved å betrakte homogene polynom av tredje grad. Av dem er der fire, nemlig , ,  og . Men bare de to kombinasjonene  og  oppfyller Laplace-ligningen. I polarkoordinater kan de skrives som henholdsvis  og  når man benytter [[Abraham de Moivre|de Moivres formler]]. Og slik fortsetter det for polynom av høyere grad. Disse spesielle løsningene sies å være *'harmoniske** da de kan uttrykkes ved [[trigonometrisk funksjon|trigonometriske funksjoner]] som beskriver [[harmonisk oscillator|harmoniske svingninger]].

===Separasjon av variable===

Mer generelle løsninger kan finnes ved å anta at løsningene kan skrives som et produkt av to funksjoner hvor hver kun avhenger av en variabel. Dette kalles for metoden med «separasjon av variable».

Med denne antagelsen &Phi;(''x,y'') = ''X''(''x'')''Y''(''y'') kan den to-dimensjonale Laplace-ligningen omformes til

: <math> {X''(x)\over X(x)} = - {Y''(y)\over Y(y)} </math>

Da ''x'' og ''y'' er uavhengige variable, kan denne ligningen bare være gyldig hvis hver side er lik med en konstant. Kalles den for ''K'', vil dermed den opprinnelige Laplace-ligningen være redusert til de to vanlige differensialligningene

: <math> X'' = KX, \; \; \; Y'' = - KY </math>

som er enklere å løse. For eksempel, hvis konstanten ''K'' er positiv slik at man kan skrive ''K = k''<sup>2</sup>, vil  og  gi den enkle løsningen  av den opprinnelige Laplace-ligningen. Mer generelle løsninger finnes ved å kombinere flere slike løsninger med verdier for ''K'' som er valgt for å oppfylle gitte grensebetingelser.

===Polarkoordinater===

[[Laplace-operator]] tar forskjell form avhengig av hvilke koordinatsystem man benytter. I to dimensjoner kan det ofte være hensiktsmessig å benytte [[polarkoordinatsystem|polarkoordinater]] . Laplace-ligningen tar da formen

: <math> \frac{1}{r} \frac{\partial}{\partial r} \left( r \frac{\partial V}{\partial r} \right) + \frac{1}{r^2} \frac{\partial^2 V}{\partial \phi^2} = \left( \frac{\partial^2}{\partial r^2} + \frac{1}{r} \frac{\partial }{\partial r} + \frac{1}{r^2} \frac{\partial^2}{\partial \phi^2}\right)V(r,\phi) = 0 </math>

Løsninger kan igjen finnes ved separasjon av variable. Man skriver da ''V'' = ''R''(''r'')''F''(''&phi;'') og setter inn i ligningen. Det vil da gi at

: <math> {1\over R} (r^2 R'' + r R') = - {F''\over F} =  \lambda </math>

hvor ''&lambda;'' er en konstant. Skal løsningene være periodiske i ''&phi;'', må ''&lambda;'' =  ''n''<sup>&thinsp;2</sup> hvor ''n'' er et helt positivt tall eller null. Vinkelavhengigheten til løsningene må derfor være

: <math> F = \begin{cases} a_0,  \; \;  n = 0   \\ a_n\cos n\phi + b_n\sin n\phi, \; \; n = 1,2,3,\ldots \end{cases}</math>

Den resterende, radielle ligningen er nå

: <math> r^2R'' + rR' - n^2R = 0 </math>

som har løsningene

: <math> R = \begin{cases} c_0 + d_0\ln r, \; \;  n = 0  \\ c_nr^n + d_nr^{-n}, \; \; n = 1,2,3,\ldots \end{cases}</math>

Her er ''a<sub>n</sub>''&thinsp;, ''b<sub>n</sub>''&thinsp;, ''c<sub>n</sub>''&thinsp; og  ''d<sub>n</sub>'' integrasjonskonstanter som må bestemmes ut fra grensebetingelsene. For eksempel, hvis man vil at løsningene skal være regulære i punktet , da må alle   slik at den mest generelle løsningen tar formen

: <math> V = a_0 + \sum_{n=1}^\infty r^n (a_n\cos n\phi + b_n\sin n\phi) </math>

hvor man her har redefinert konstantene slik at ''c<sub>n</sub>'' = 1. Hvis man lar radius ''r'' = konstant, finner man her de harmoniske løsningene cos&thinsp;''n&phi;'' og sin&thinsp;''n&phi;'' som oppsto på en mer komplisert måte som homogene polynom ved bruk av kartesiske koordinater.

==Tre dimensjoner==
Ved bruk av [[kartesisk koordinatsystem|kartesiske koordinater]] i  ''N'' = 3 dimensjoner har Laplace-ligningen formen

: <math> \left({\partial^2\over \partial x^2} + {\partial^2\over \partial y^2} + {\partial^2\over \partial z^2}\right)\Phi(x,y,z) = 0 </math>

Som for to dimensjoner kan man finne løsninger med form av [[polynom|homogene polynom]] i de tre variable ''x'', ''y'' og ''z''.  Disse har formen  og kombinasjoner av slike hvor heltallene ''p'', ''q'' og ''r&thinsp;'' har en konstant sum som er graden til polynomet. Har disse polynomene grad ℓ , er det i allminnelighet  slike løsninger av Laplace-ligningen. De danner [[sfærisk harmonisk funksjon|sfærisk harmoniske funksjoner]]  ,  også kalt for **kuleflatefunksjoner**.

Av laveste grad ℓ = 1 har man tre løsninger som ganske enkelt er de tre koordinatene ''x'', ''y'' og ''z''.  For ℓ = 2 finner man de fem løsningene ,   samt .
Uttrykkes disse i  [[kulekoordinater]] (''r, &theta;, &phi;'') slik at ,   og  ser man at denne siste kombinasjonen  tilsvarer ''Y''<sub>20</sub>. På samme måte danner ''xy'' og  de to sfæriske harmoniske funksjonene , mens  ''xz'' og ''yz'' kan kombineres til å gi  Med bruk av [[komplekst tall|komplekse variable]] viser man lett at det homogene polynomet  er direkte proporsjonal med  den spesielle kuleflatefunksjonen 

===Kulekoordinater===

Løsninger av Laplace-ligningen som involverer sfærisk harmonisk funksjoner kommer mer direkte frem ved å uttrykke [[Laplace-operator]]en i kulekoordinater. For potensialet  har den da formen

: <math> \nabla^2 V = {1\over r^2} {\partial\over\partial r}\left(r^2 {\partial V\over\partial r}\right) + {1\over r^2\sin\theta} {\partial\over\partial\theta}\left(\sin\theta{\partial V \over\partial\theta}\right) + {1\over r^2\sin^2\theta} {\partial^2 V\over\partial\phi^2} = 0 </math>

Ved en første separasjon av variable kan man skrive ''V'' = ''R''(''r'')''Y''(''&theta;,&phi;''). Det gir

: <math> {1\over R} (r^2 R'' + 2r R') = - {1\over Y}\left[{1\over \sin\theta}{\partial\over\partial\theta}\left(\sin\theta{\partial Y \over\partial\theta}\right) + {1\over \sin^2\theta} {\partial^2 Y\over\partial\phi^2}\right] =  \lambda </math>

hvor ''&lambda;'' må være en konstant. For at den  [[partiell differensialligning|partielle differensialligningen]] på høyre side skal ha regulære løsninger, må   hvor ℓ er et helt, positivt tall. Den går da over til differensialligningen som definerer [[sfærisk harmonisk funksjon|sfærisk harmoniske funksjoner]]. De betegnes med  hvor heltallet ''m'' må ligge i intervallet  og derfor tar  forskjellige verdier.

Den radielle ligningen tar nå formen

: <math> r^2R'' + 2rR' - \ell(\ell + 1) R = 0 </math>

og har løsninger av formen ''r''<sup>&thinsp;ℓ</sup> eller 1/''r''<sup>&thinsp;ℓ + 1</sup>. Disse siste divergerer når ''r'' &rarr; 0. Disse løsningene har derfor den generelle formen

: <math> V(r, \theta, \phi) = \sum_{\ell=0}^\infty \sum_{m=-\ell}^\ell \Big(A_{\ell m} r^\ell +{B_{\ell m}\over r^{\ell +1}} \Big) Y_{\ell m} (\theta, \phi)  </math>

hvor ''A''<sub>ℓ''m''</sub> og ''B''<sub>ℓ''m''</sub> er konstanter som må bestemmes ut fra grensebetingelsene. For eksempel, hvis man vil at løsningene skal være regulære i grensen ''r'' &rarr; 0, må man ha at .

===Sylinderkoordinater===

Noen løsninger av Laplace-ligningen finnes mer direkte i [[polarkoordinatsystem|sylinderkoordinater]]  (''r, &phi;, z'') der ''x'' = ''r''&thinsp;cos''&phi;'' og ''y'' = ''r''&thinsp;sin''&phi;''. Fra uttrykke for den tilsvarende [[Laplace-operator]]en tar ligningen da formen

: <math>  \nabla^2 V = \frac{1}{r} \frac{\partial}{\partial r} \left(r \frac{\partial V}{\partial r} \right) + \frac{1}{r^2} \frac{\partial^2 V}{\partial \phi^2} + \frac{\partial^2 V}{\partial z^2 } = 0 </math>

Igjen kan løsninger finnes ved en separasjon av variabler. Ved å anta formen , vil ''&phi;''-avhengigheten være gitt av løsninger til den enklere differnsialligningen  hvor på venstre side opptrer den dobbeltderiverte, her med hensyn på ''&phi;''. For at den skal ha skal periodiske løsninger, må  hvor ''m'' er et positivt heltall. Det gir

: <math> F_m(\phi) = A_m\cos m\phi + B_m\sin m\phi </math>

hvor ''A<sub>m</sub>'' og ''B<sub>m</sub>&thinsp;''  er kontanter.  Likedan vil ''z''-avhengigheten følge fra differensialligningen  hvor ''K'' er en foreløbig ubestemt [[parameter]].  Hvis man ønsker at ligningen skal ha løsninger som ikke er periodiske i ''z''-retningen, må  .  Da vil løsningene av denne differensialligningen kunne skrives som

: <math> Z(kz) = a(k)e^{-kz} + b(k) e^{kz} </math>

der integrasjonskonstantene ''a'' og ''b'' avhenger av grensebetingelsene. Hele Laplace-ligningen involverer nå bar den radielle funksjonen ''R'' som må oppfylle

: <math> r^2 R'' + r R' + (k^2r^2 - m^2)R = 0 </math>

Dette er [[Bessel-funksjon|Bessel-ligningen]] med løsninger som er gitt ved de to [[Bessel-funksjon]]ene  og . Denne siste er ikke regulær i grensen ''r'' &rarr; 0. De kan kombineres med to ubestemte konstanter ''C<sub>m</sub>'' og ''D<sub>m</sub>&thinsp;'' til å gi den generelle, radielle løsningen.

: <math> R_m(kr) = C_m(k) J_m(kr) + D_m(k)Y_m(kr) </math>

Den fulle løsningen av Laplace-ligningen i sylinderkoordinater er nå gitt som produktet av de tre delløsningene ''R''(''r''),  ''F''(''&phi;'') og  ''Z''(''z'') summert over alle mulige verdier for ''k'' og ''m''. Integrasjonskonstantene vil i allminnelighet også avhenge av disse verdiene.

Som et enkelt eksempel på en fullstendig løsning kan man tenke seg at den skal være regulær i origo ''r'' = 0. Det betyr med en gang at alle . Mulige verdier for ''k'' kan så finnes ved å forlange en viss verdi for den radielle funksjonen ''R'' for en endelig verdi av ''r''.  Hvis ikke, kan ''k'' variere fritt mellom null og uendelig og man må integrere løsningen over alle disse verdiene. Videre kan man forlange at løsningen ikke skal divergere for store verdier av koordinaten ''z'', noe som medfører at alle konstantene . Dermed vil den generelle løsningen i dette tilfellet kunne skrives som

: <math> V(r,\phi, z) = \sum_{m = 0}^\infty\int_0^\infty dk e^{-kz}J_m(kr)\big(A_m(k)\cos m\phi + B_m(k)\sin m\phi\big) </math>

==Middelverditeoremet==

Løsningene av Laplace-ligningen er [[harmonisk funksjon|harmoniske funksjoner]]. En av de viktigste egenskaper til disse, er at slike funksjoner tilfredsstiller et generelt [[teorem]] om deres midlere verdi i hvert punkt.

Dette kan nesten trivielt illustreres i det endimensjonale tilfellet for Laplace-ligningen ''d''<sup>&thinsp;2</sup>&Phi;/''dx''<sup>2</sup> = 0. Den mest generelle løsningen er  hvor ''a'' og ''b'' er konstanter. Ved direkte innsetting i løsningen ser man da at verdien til funksjon &Phi; i et vilkårlig punkt ''x'' er gitt ved verdiene i to nabopunkt som

: <math> \Phi(x) = {1\over 2}\Big(\Phi(x + h) + \Phi(x - h)\Big) </math>

hvor ''h'' her er vilkårlig stor eller liten. Dette er  den [[middelverdi]]en av de to verdiene i nabopunktene. Man ser dette også geometrisk da løsningen representerer en [[linje|rett linje]].

===Diskret versjon===

Teoremet kan lett sannsynliggjøres ved å erstatte [[Laplace-operator]]en med en diskret versjon. For en funksjon ''f''(''x'') er den [[derivasjon|deriverte]] i et punkt definert som

:<math>f'(x) = {1\over h} \big(f(x + h) - f(x)\big) = {1\over h} \big(f(x) - f(x - h)\big)</math>

i grensen der ''h'' &rarr; 0. Lar man istedet ''h'' være tilstrekkelig liten og endelig, er dette definisjonen for en «diskret derivasjon» av funksjonen. Den andrederiverte kan man da skrive som

:<math>f''(x) = {1\over h} \big(f'(x) - f'(x - h)\big) = {1\over h^2} \big(f(x + h) + f(x -h) - 2f(x)\big)</math>

Tilfredsstiller denne funksjonen den endimensjonal Laplace-ligningen ''f''&thinsp;"(''x'') = 0, vil derfor igjen ''f''(''x'') være gitt ved middelverdien i de to nabopunktene.

Ved bruk av denne diskrete derivasjonen for en funksjon ''F''(''x,y'') som tilfredsstiller den todimensjonale Laplace-ligningen, vil man da på samme måte ha

: <math> F(x,y) = {1\over 4}\big(F(x+h,y) + F(x-h,y) + F(x,y+h) + F(x,y-h)\big) </math>

Denne verdien er igjen middelverdien av de fire naboverdiene. Man kan utvide dette resultatet på samme måte til vilkårlige høyere dimensjoner, men størrelsen ''h'' må her være tilstrekkelig liten.

===Middelverditeoremet===

[[Harmonisk funksjon|Harmoniske funksjoner]] mer generelt har en meget spesiell egenskap som er uttrykt ved «middelverditeoremet». Det sier at for en generell løsning &Phi;(**x**) av Laplace i ''N'' dimensjoner hvor punktet **x** er angitt ved ''N'' koordinater , så er verdien av funksjonen i hvert punkt lik med den midlere verdi av funksjonen på en vilkårlig stor kuleflate rundt det samme punktet. På matematisk form kan dette uttrykkes som

: <math> \Phi(\mathbf{x}) =  {1\over S_{N-1}(r)}\oint_{\!S}\Phi(\mathbf{y}) dS </math>

hvor **y** angir punkt på kuleflaten ''S'' omkring **x** med vilkårlig stor radius ''r''. For en generell kule i [[Kule (geometri)#Kuler i høyere dimensjoner|''N'' dimensjoner]],  vet man at den har et areal

: <math>S_{N-1} =\frac{2 \pi^{N/2}}{\Gamma(N/2)} r^{N-1}</math>

For den mest vanlige situasjon med ''N'' = 3 dimensjoner gir denne generelle formelen det kjente resultatet  ''S''<sub>2</sub> = 4''&pi;&thinsp;r''<sup>&thinsp;2</sup>. Hvis punktet **x*' er i [[origo]] og man benytter [[kulekoordinater]] slik at flateelementet  sier teoremet at verdien der er

: <math> \Phi(0) = {1\over 4\pi}\oint_{\!S}\Phi(r,\theta,\phi)\sin\theta d\theta d\phi </math>

uavhengig av verdien til radius *r*. Dette viser med klarhet at harmoniske funksjoner må ha veldig spesielle egenskaper.

===Elektrostatisk eksempel===

En [[elektrisk ladning|elektrisk punktladning]] *Q* er plassert på *z*-aksen i avstand *z* fra origo. Det [[elektrisk potensial|elektriske potensialet]] i origo er derfor  der *ke* er [[Coulombs konstant]]. Ifølge middelverditeoremet skal dette være lik den midlere verdien av potensialet på en kuleflate med sentrum i origo. Har den radius *r*, er potensialet i hvert punkt på denne flaten  hvor

: <math> R = \sqrt{z^2 + r^2 - 2zr\cos\theta} </math>

når man benytter [[cosinussetningen]] og hvert punkt på flaten er angitt ved de sfæriske koordinatene (*&theta;,&phi;*). Den søkte middelverdien er dermed

: <math> \Phi(0) = {1\over 4\pi}\oint_{\!S} {k_e Q\over \sqrt{z^2 + r^2 - 2zr\cos\theta}}\sin\theta d\theta d\phi </math>

Integranden er uavhengig av den asimutal vinkelen *&phi;* slik at integrasjon over denne variable gir ganske enkelt 2*&pi;&thinsp;*. Integrasjonen over den polare vinkelen *&theta;* forenkles ved  å la *x* = cos*&theta;* være ny integrasjonsvariabel. Det gir

: <math>\begin{align} \Phi(0) &= {k_e Q\over 4\pi}{4\pi\over 2zr} \sqrt{z^2 + r^2 - 2zr\cos\theta}\Big|_0^\pi \\
&= {k_e Q\over 2zr}\Big[ (z + r) - (z - r)\Big] = {k_e Q\over z} \end{align}</math>
 
som er uavhengig av radius *r* og bekrefter middelverditeoremets korrekthet. For at Laplace-ligningen skal være gyldig innenfor kuleflaten,  må ladningen *Q* ligge utenfor denne. Det betyr her at .

===Maksima og minima===
Middelverditeoremet legger store krav på formen til harmoniske funksjoner. Det viktigste er at de ikke kan ha noe lokalt maksimum eller minimum innenfor et begrenset område. Slike ekstremalpunkt kan eventuelt kun finnes på randen av området. Grunnen er at hvis det for eksempel fantes et lokalt maksimum inni området, så kunne ikke dette være middelverdien av punktene på en liten kuleflate som omslutter punktet da alle verdiene på denne ville være mindre enn verdien i det gitte punktet. Tilsvarende ville gjelde for et lokalt minimum.

== Referanser ==
<references/>

== Litteratur ==
* J. Mathews and R.L. Walker, *Mathematical Methods of Physics*,  W.A. Benjamin, New York (1970). ISBN 0-8053-7002-1.
* M.L. Boas, *Mathematical Methods in the Physical Sciences'', John Wiley & Sons, New York (1983). ISBN 0-471-04409-1.

## Eksterne lenker
*