<div dir="ltr"><div>Interesting; even more  when using rounded numbers returns a correct solution!<br></div><div><br></div><div>To
 protect against such errors, maybe `matrix-solve` could run a 
post-check to verify that M×X=B up to some epsilon, depending on an optional argument?<br>Or maybe just 
mention to do this check systematically in the docs?<br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Apr 16, 2014 at 1:31 PM, Jens Axel Søgaard <span dir="ltr"><<a href="mailto:jensaxel@soegaard.net" target="_blank">jensaxel@soegaard.net</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="">2014-04-16 12:18 GMT+02:00 Jens Axel Søgaard <<a href="mailto:jensaxel@soegaard.net">jensaxel@soegaard.net</a>>:<br>


</div><div class="">> Hi Laurent,<br>
><br>
> I think the underlying problem is that the matrix is *very* close to<br>
> an invertible one:<br>
<br>
</div>non-invertible<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
>> (matrix-determinant<br>
>      (matrix [[ 1     0     9/10  1]<br>
>                  [ 0     1     1/10  1]<br>
>                  [ 9/10  1/10 82/100 1]<br>
>                  [ 1     1      1   0]]))<br>
> 0<br>
</div></div></blockquote></div><br></div>