Re: Eigensystem
- To: mathgroup at smc.vnet.net
- Subject: [mg73317] Re: Eigensystem
- From: Jean-Marc Gulliet <jeanmarc.gulliet at gmail.com>
- Date: Mon, 12 Feb 2007 04:59:57 -0500 (EST)
- Organization: The Open University, Milton Keynes, UK
- References: <eqmcbn$bbe$1@smc.vnet.net>
Vadim N. Smelyanskiy wrote:
<snip>
> Please can you advice me how to run Eigensystem function so that its result
> has precision (say 50)
>
> that is much higher then the Machine Precision.
Vadim,
Use either exact numbers (see m1) or SetPrecision (see m3).
In[1]:=
{lambda, mu} = Eigenvalues[m1 = {{2, 1}, {3, -1}}]
Precision /@ {m1, lambda, mu}
Out[1]=
1 1
{- (1 + Sqrt[21]), - (1 - Sqrt[21])}
2 2
Out[2]=
{Infinity, Infinity, Infinity}
In[3]:=
{lambda, mu} = Eigenvalues[m2 = {{2., 1.}, {3., -1.}}]
Precision /@ {m2, lambda, mu}
Out[3]=
{2.79129, -1.79129}
Out[4]=
{MachinePrecision, MachinePrecision, MachinePrecision}
In[5]:=
{lambda, mu} = Eigenvalues[
m3 = SetPrecision[{{2., 1.}, {3., -1.}}, 50]]
Precision /@ {m3, lambda, mu}
Out[5]=
{2.7912878474779200032940235968640042444922282883840,
-1.7912878474779200032940235968640042444922282883840}
Out[6]=
{50., 50., 50.}
Regards,
Jean-Marc