[Date Index] [Thread Index] [Author Index]

stable precision

• To: mathgroup at smc.vnet.net
• Subject: [mg24246] stable precision
• From: pecar at masfak.masfak.ni.ac.yu (Predrag Rajkovic)
• Date: Mon, 3 Jul 2000 20:39:26 -0400 (EDT)
• Organization: The Math Forum
• Sender: owner-wri-mathgroup at wolfram.com

Hello!

How can I founded chosen precision through whole programm? 
I do not want to use every time N[]?
In FORTRAN, it can be done by REAL*8, COMPLEX*16.
Does such statement exist in Mathematica?  
For example, I want 10^-32 in next programm.
(
	f[z_]:=Sin[z]*Log[z+2]-1;
	x=-1.9;Print[x,"    ",f[x]];
	y=-1.1;Print[y,"    ",f[y]];
Label[Start];
z=y-(y-x)/(f[y]-f[x])*f[y];
	Print[z,"    ",f[z]];
	x=y;y=z;				
If[Abs[x-y]>10^ -5,  Goto[Start] ]
	)