Pathological floating point problems: Difference between revisions

Content deleted Content added
Line 460:
 
Task 1: Define recursive function V(n):
<lang parigp>V(n,a=2,v=0-4.)=if(n < 43,return(v[3]),!v,v=[2,-4,18.5]);V(n--,v=[v[2],v[3],111-1130/v[3]+3000/(v[3]*v[2]a)])</lang>
In order to work set precision to at least 200 digits:
<pre>\p 200: realprecision = 211 significant digits (200 digits displayed)