别争了
for(i=0;i<n;i++)
{
l_xx=(arr[i,0]-avg_x)*(arr[i,0]-avg_x);
l_yy=(arr[i,1]-avg_y)*(arr[i,1]-avg_y);
l_xy=(arr[i,0]-avg_x)*(arr[i,1]-avg_y);
}
不谈数学上，，就程序上看上面这个就不对，，，，


这样得到的永远是数组里最后一个值的平方

for(i=0;i<n;i++)
{
l_xx+=(arr[i,0]-avg_x)*(arr[i,0]-avg_x);
l_yy+=(arr[i,1]-avg_y)*(arr[i,1]-avg_y);
l_xy+=(arr[i,0]-avg_x)*(arr[i,1]-avg_y);
}
l_xx=l_xx/n;
l_yy=l_yy/n;
l_xy=l_xy/n;

计算机数值方法,哎,好怀念!!
当时学这课的时候光睡觉了!!!

遗传算法 , 很有意思,当年上学时,我有个导师专门研究遗传算法的,我也学了一点,但工作以后,一直都没用上,