Algorytm liczb dużej precyzji lub BCmath Opcje

[s_w_ir]

Witam.
Zapewne wykaże się wręcz lamerską niewiedzą, ale interesuje mnie gdzie moge dostać(znaleźć) źródło funkcji(wielu nie jednej (IMG:http://forum.php.pl/style_emoticons/default/winksmiley.jpg) ) bcmath bądź dowiedzieć się jak jest skonstruowany programowo typ double, nie chodzi mi tu o mantysy itp. sprawy , tylko bardziej jak to prgramowo wysymulować.
Mój problem polega na tym że musze napisać program wyliczający liczbe e z dowolną precyzją w C. Wiadomo że żaden typ nie daje nieskończonej precyzji(mój też nie będzie), ale jakoś musze sie do tego zabrać. Myślałem żeby każdą kolejną liczbe po przecinku wciskać w tablice(osobne pole), ale żeby to zrobić musze jakoś najpierw dostać te liczby z powyżej 15 miejsca po przecinku.
Inna opcja to złożyć jakiś własny typ zmiennych, tylko czy da rade złożyć go z kilku double które razem dadzą np. 50 miejsc po przecinku?
Jak uda mi się coś wykombinować to dam znać (IMG:http://forum.php.pl/style_emoticons/default/smile.gif)
Pozdrawiam.

Ten post edytował s_w_ir 5.06.2006, 11:10:31

[s_w_ir]

Dziękuje, ale program musze napsisać i chce(naucze się pewnej rzeczy o której od dawna myślałem). Będe poszukiwał w kierunku wielomianów w takim razie.

EDit. Jakoś nie jestem sobie w stanie wyobrazić nawet jak można to programowo stworzyć, tzn. jak obliczyć 16, i dalej, miejsce po przecinku podczas gdy long double przyjmuje tylko 15 miejsc po przecinku. Do liczby Pi znalazłem wzór na obliczanie k-tej liczby po przecinku(zresztą moge obliczyć Pi lub e to nieważne, chodzi o to żeby przedstawić sposób przechowywania tych liczb), ale to nie jest to czego szukam. Może byc nawet coś bardzo obciążające procesor, ale byle można to uznać za metodę.....
....pare skeund później....
właśnie wpadłem na pomysł żeby liczyć to(1+1/n)^n i co pare(co każdą?) iterację, mnożyć to przez 10 i początek wsadzać do tablicy to co znika z lewej strony. Zresztą rozpiszę to zaraz na kartce i się okaże czy to możliwe...

Ten post edytował s_w_ir 6.06.2006, 20:34:12