6. Romberg İntegrali m: Mathematica İle Sembolik Bir Yaklaşım
2018'de E-ATA Algoritmaları'nın programlanması için bizdeki uzmanlardan yardım istemiş ve ancak bir profesörle telefonda görüşebilmiştim. Ben elimde 2003'ten kalma "E-ATA Algoritmaları Ver. 1" olduğunu ve bunların programlanması gerektiğini söylemiştim. Bana bunları nerede yayımladığımı sordu. Ben de o sırada (2002-2003) kendi kişisel web sitemde yayımladığımı söyledim. Bana inanmadı ama algoritmaların adı onu bayağı bir heyecanlandırmıştı. Bana "hadi bakalım, görelim!" diyerek aslında onu heyecanlandıran şeyin Atatürk olduğunu anladım (Bkz. "RİK3: Romberg'in Türkiye Şubesiyle Görüşmelerim!", S. 35). Bu durum şimdi yukarıda gördüğünüz gibi Atatürk'ün yapay zekayla hazırlanmış fotoğraflarıyla daha da açığa çıktı. Çünkü Atatürk'ü daha önceden bu şekilde hiç görmemiştik!
Bu programda "m" kodlamasından anlaşılacağı 4 ve 5. programdakilerin genel şeklini verdim. Fakat bu, programda "p" parametresi olarak geçer. Programın girişinde p=5 olarak aldım, dolayısıyla n=25 toleransı gerekti. Eğer sistemininiz güçlü değilse p'yi 5'ten küçük alın. Çünkü sistemim güçlü olmasına rağmen (ki sistemimi RİK 4'ün 45. sayfasının sonunda görebilirsiniz) 3. sekmedeki 2. tablodaki son çıktının alınması 17970.6 Sn=5 Saat sürdü. Bir önceki ise 3533.59 Sn=1 Saat sürdü. Bu durumda son tablonun çıktısının ne kadar zaman aldığını sormanıza gerek kalmaz. Çünkü onun için 5 saatten fazla bir süre bekledim.
Programlama sorunu mu?
İşte uzun süren bu beklemelerin nedeni programdan kaynaklanıyor. Bu nedenle size tavsiyem, eğer sistemininiz güçlü değilse programı 4. sekmedeki gibi iteratif olarak kullanmanız. Çünkü bu şekildeki hesaplar saniye bile sürmez!
2018'de konuyla ilgili resim altındaki uzmanlarla görüşmüştüm ve m=1 için Richardson ekstrapolasyonunda uzun süre yanıt vermeyen bu beklemeler için bir çözüm bulamadıklarını, bu durumun programdan kaynaklandığını söylediler. Belki ileride bir çözüm bulabiliriz. Fakat ben her programda 4. sekmeyi koyarak bu probleme kesin bir çözüm buldum!