Re: Coaxing N[] to work
- To: mathgroup at smc.vnet.net
- Subject: [mg73381] Re: [mg73370] Coaxing N[] to work
- From: bsyehuda at gmail.com
- Date: Thu, 15 Feb 2007 04:54:33 -0500 (EST)
FullForm[lambda] will reveal more digits, but in general, your example works
with machine precision numbers, and this depends on your machine
One possibility is to set the precision of the numbers, that is 1.0`30 (30
digits precision)
or to work with infinite precision, that is, integers and rationals, with no
decimal point
for the second case
x = {2, 3, 5};
A = {{6, 2, 1}, {2, 3, 1}, {1, 1, 1}};
For[k = 0, k < 15, ++k, lambda = x.A.x/(x.x);
y = LinearSolve[A, x];
x = y/Norm[y, Infinity];]
N[lambda, 30]
returns
0.578933385691052787623495851172
regards
yehuda
On 2/14/07, p at dirac.org <p at dirac.org> wrote:
>
> Sometimes N[,] doesn't appear to work. Like here:
>
>
> x = {2.0, 3.0, 5.0};
> A = { {6.0, 2.0, 1.0}, {2.0, 3.0, 1.0}, {1.0, 1.0, 1.0} };
> For[ k=0, k<15, ++k,
> lambda = x.A.x/(x.x);
> y = LinearSolve[A,x];
> x = y / Norm[y,Infinity];
> ]
> N[lambda, 30]
>
>
> The output is:
>
> Out[5]= 0.578933
>
> I was expecting 30 digits. Why did N[] ignore my request for 30 digits?
>
>