Błędy numeryczne

Podczas składania funkcji
f[x_]:=x^2
dla liczb o ustawionej przez użytkownika precyzji (x0=N[1,prec]), precyzja wyniku spada bardzo gwałtownie. Jak obejść ten błąd? Przy precyzji maszyny ten błąd się nie kumuluje. Przykład poniżej:

s1 = N[1];
s2 = N[1, 10];
TableForm[
 N /@ Precision /@ NestList[#^2 &, #, 20] & /@ {s1, s2},
 TableDirections -> Row,
 TableHeadings -> {{"Precision[x0]=MachinePrecision", 
    "Precision[x0]=10"}}]
ListPlot[%]

Ostatnio edytowany przez pawel.biernat (2010-01-30 14:40:21)