很偶然地发现了一个matlab2010a的一个bug：在某种非常特殊的情况下，matlab2010a的向量内积，即*指令会产生错误的结果。

下面陈述具体现象：

1.程序源代码

(注：运行该程序时需要把den1.mat和kern1.mat这两个matlab数据文件与程序放在同一个文件夹中)

clc;

clear;

load('den1.mat');

load('kern1.mat');

x = kern;

y = den(:);

z = 0;

for i = 1:length(y)

z = z + x(1,i)*y(i,1);

end

z1 = x*y;

z2 = sum(y);

2.程序相关说明

(1)den1.mat是matlab数据文件(可从此处下载http://ishare.iask.sina.com.cn/f/22182234.html)，其中数据是1*36维行向量，如下(其中的 - 是负号)：

[1 -34.9104720617963 591.960057485552 -6494.90695093987 51826.3768981548 -320501.857479884 1598411.70513468 -6605060.33215733 23058608.5581138 -68998984.4665517 178938782.362431 -405639572.104363 809205656.532103 -1428012719.05099 2238285526.53077 -3125592609.96192 3897008179.02336 -4344352067.53865 4333253669.23331 -3867217670.62412 3085871743.18171 -2198564376.37758 1395513287.06250 -786754091.156305 392372548.725918 -172203104.526507 66062863.9131647 -21964736.3140818 6259618.77763060 -1507085.68697902 300647.808748368 -48367.8227891117 6030.54857993965 -546.832712826108 32.0845535965159 -0.914362645545420]

(2)kern1.mat也是matlab数据文件(可从此处下载http://ishare.iask.sina.com.cn/f/22182235.html)，其中数据是1*36维行向量，如下：

[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1]

(3)den和kern分别是来自于den1.mat和kern1.mat的1*36维行向量，其中kern的所有元素都为1。

(4)很显然这个程序的作用是用三种程序不同的程序方法计算den中所有元素的和，即：

z1使用向量内积计算，z2直接求元素和，z是用循环的方法求向量的内积。

3.bug现象描述

(1)在matlab7(Version 7.0.0.19920(R14))下运行该程序所得z=z1=z2=5.3111e-008的完全相同，这符合我们的数学常识。

但是令人奇怪的是，在matlab2010a中出现了奇怪的结果。

(2)在matlab210a中，z=z2=5.311088491222193e-08，z1=0.

是不是很奇怪！难道这是matlab2010a的bug！

(3)matlab210b与matlab210a出现同样问题。