Re: Zero or not?

[Bob Weigel]

On Saturday 21 July 2007 08:50, Matthias Brennwald wrote:
> Dear all
>
> I expected the below code to yield only zeros, but it does not:
>
> ----------------------------
> u = [ 0 , 1e-3 , 1e-4 ; 1e-3 , 0 , 0 ; 0 , 1e-4 , 0 ];
> d = u - diag(sum(u));
> sum(d)
> ----------------------------
>
> The output is:
>    0.0000e+00  -4.0658e-20   0.0000e+00
>
> Why is this? Is there something I can do to 'fix' this?
>

This is correct.  

See

http://docs.sun.com/app/docs/doc/800-7895

or

http://www.nabble.com/forum/Search.jtp?forum=1895&local=y&query=goldberg