From the Dirichlet Integral to Lobachevsky's Formula: A Formalization in Lean 4

In het uitgestrekte landschap van de wiskunde is er een rustig hoekje gewijd aan de studie van hoe dingen in de loop van de tijd bij elkaar optellen, vooral wanneer die dingen heen en weer wiebelen. Dit is het domein van de reële analyse, waar wiskundigen het gedrag onderzoeken van functies die continu veranderen. Een van de beroemdste puzzels in dit veld betreft een specifieke curve die op en neer gaat als een golf, en steeds kleiner wordt naarmate deze zich uitstrekt naar oneindigheid. De vraag is eenvoudig te stellen maar lastig op te lossen: als je de oppervlakte onder deze wiebelende curve optelt vanaf het beginpunt tot het verste punt dat je je kunt voorstellen, wat totaal krijg je dan? Al meer dan een eeuw weten wiskundigen het antwoord, maar het rigoureus bewijzen ervan zonder verborgen aannames te maken, is altijd een delicate taak geweest. Dit komt omdat de curve niet snel genoeg tot rust komt voor standaard optelregels om direct toe te passen; het vertrouwt op een precieze annulering van positieve en negatieve oppervlaktes om tot een eindige som te komen. Het begrijpen van dit gedrag is niet alleen cruciaal voor de zuivere wiskunde, maar ook voor de technologie die de basis vormt van moderne communicatie, aangezien deze zelfde wiebelende patronen worden gebruikt om signalen en afbeeldingen te reconstrueren uit ruwe gegevens.

Onlangs besloten twee onderzoekers, Daniel Goldberg en Antoine Vinciguerra, deze klassieke kwestie aan te pakken met een computerprogramma dat ontworpen is om wiskundige bewijzen met absolute zekerheid te controleren. Ze schreven niet alleen de oplossing op; ze bouwden een volledig, stapsgewijs logisch argument binnen een softwaresysteem genaamd Lean 4, dat fungeert als een onvermoeibare auditor die geen enkele stap accepteert tenzij deze gerechtvaardigd is door de regels van de logica. Hun doel was om de Dirichlet-integraal te formaliseren, de naam gegeven aan de berekening van die specifieke wiebelende oppervlakte, en om te laten zien hoe deze verbonden is met een breder scala aan regels voor het integreren van periodieke functies. De uitdaging waarmee zij werden geconfronteerd, was dat de standaard manier waarop computers oppervlakteberekeningen afhandelen, bekend als de Lebesgue-integraal, deze specifieke curve niet direct kan verwerken omdat de totale omvang van haar wiebelingen oneindig is, ook al is de netto oppervlakte eindig. Om dit te omzeilen, moesten de onderzoekers een slimme omweg vinden die het oneindige probleem vermeed, terwijl ze toch tot het juiste antwoord leidden.

In plaats van te proberen de computer de oorspronkelijke wiebelende curve direct te laten accepteren, keken het team eerst naar een aangepaste versie hiervan waarbij de curve gekwadrateerd werd. Deze gekwadrateerde versie gedraagt zich veel beter; de totale oppervlakte ervan is eindig en goed beheersbaar, waardoor de computer deze met standaard methoden kan berekenen. De onderzoekers bewezen vervolgens een specifieke relatie tussen de oppervlakte onder de oorspronkelijke wiebelende curve en de oppervlakte onder deze gekwadrateerde versie. Door de oppervlakte van de gekwadrateerde curve eerst te berekenen, konden ze dat resultaat wiskundig terugbrengen naar het oorspronkelijke probleem. Deze aanpak stelde hen in staat om de moeilijkheden van conditionele convergentie, waarbij de volgorde van optellen ertoe doet, te omzeilen en tot het beroemde resultaat te komen dat de totale oppervlakte exact de helft van het getal pi is. Dit was geen gok of een simulatie; het was een rigoureus bewijs dat de limiet van de oppervlakte, naarmate de grens verder en verder uitwaaiert, convergeert naar deze specifieke waarde.

Nadat ze de hoofdpuzzel hadden opgelost, gebruikten het team hun nieuwe instrumenten om te verkennen wat er nog meer uit afgeleid kon worden. Ze lieten zien hoe deze integraal fungeert als een filter dat een vloeiende, continue golf kan omzetten in een scherpe, stap-achtige sprong, een gedrag dat fundamenteel is voor hoe digitale signalen worden verwerkt. Ze ontdekten en bewezen ook een verzameling andere identiteiten waarbij producten van deze wiebelende functies betrokken zijn, waarmee ze lieten zien hoe verschillende frequenties interageren wanneer ze met elkaar vermenigvuldigd worden. Deze resultaten zijn niet slechts abstracte nieuwsgierigheden; ze bieden het wiskundige fundament voor het begrijpen van hoe signalen gereconstrueerd kunnen worden uit monsters, een concept dat centraal staat in de Shannon-samplingtheorema die wordt gebruikt in digitale audio- en beeldverwerking. De onderzoekers demonstreerden dat men door het begrijpen van het gedrag van deze specifieke integralen, precieze formules kan afleiden voor hoe verschillende golfpatronen combineren en elkaar opheffen.

De laatste en misschien wel meest verrassende prestatie van hun werk was de formalisering van een formule ontdekt door Nikolai Lobachevsky, een wiskundige die het best bekend staat om zijn werk in de niet-Euclidische meetkunde. Lobachevsky had een regel gevonden die mogelijk maakt dat de oppervlakte onder een wiebelende curve vermenigvuldigd met een herhalend patroon berekend kan worden door slechts naar een klein deel van dat patroon te kijken. De onderzoekers bewezen dat deze regel geldt voor elke continue, herhalende functie die een specifieke soort symmetrie bezit, waarbij ze de computer gebruikten om te verifiëren dat de oneindige som van de wiebelingen kan worden teruggebracht tot een eenvoudige berekening over een kort interval. Ze bereikten dit door aan te tonen dat elke dergelijke herhalende functie nauwkeurig benaderd kan worden door een som van eenvoudige cosinusgolven, en aangezien de regel voor elke individuele golf werkt, moet deze ook voor de hele functie werken. Dit biedt een door de machine gecontroleerd bewijs van een algemene identiteit die voorheen alleen begrepen werd via menselijke intuïtie en traditionele papier-en-potloodmethoden.

Het werk van Goldberg en Vinciguerra laat zien dat zelfs eeuwenoude wiskundige waarheden kunnen profiteren van de precisie van moderne computerverificatie. Door het probleem op te delen in beheersbare stukken en om de obstakels te navigeren die standaard integratiemethoden in de war brengen, hebben zij een solide fundament gelegd voor toekomstig onderzoek in signaalverwerking en harmonische analyse. Hun formalisering bevestigt dat de Dirichlet-integraal inderdaad de limiet is van oppervlaktes over begrensde intervallen, en het vestigt een betrouwbaar kader voor de formule van Lobachevsky. Deze prestatie suggereert dat vergelijkbare rigoureuze benaderingen kunnen worden toegepast op complexere versies van deze integralen, wat potentieel kan leiden tot nieuwe inzichten in hoe we de wiskundige structuren begrijpen die de fysieke wereld beheersen. Het artikel staat als een testament voor de kracht van het combineren van diepe wiskundige inzichten met de onverbiddelijke logica van computerverificatie, waarbij een klassieke puzzel wordt omgezet in een geverifieerd feit.