自己给自己顶一个

不好意思没有满足楼主所有要求，用matlab做的，软件比较大。问题有无穷多解，我就主观地加了一个最优化条件，函数如下：
function [ x ] = solution( beq )
%x是输出的矩阵，beq是Ma...Mj,Na...Nj组成的20*1列向量，这个问题有无穷解
%为了得到一个解就加了一个最优化条件f，f根据需要调整
a=zeros(20,100);
for i=1:10
    j=10*(i-1)+1;
    k=10*i;
    a(i,j:k)=1;
end
for i=11:20
    j=i-10;
    k=90+(i-10);
    a(i,j:10:k)=1;
end
f=ones(100,1);
b=zeros(100,1);
rstr=-eye(100);
x=linprog(f,rstr,b,a,beq);
x=reshape(x',10,10);
x
end

hope help~good luck

又用LINGO做了一下，它说license is too small for the given problem. 貌似不是程序问题，不知道是不是软件问题，我是破解版。楼主please试一下。个人感觉MATLAB合适一些，毕竟这是个解方程组，LINGO主要用来做优化的~
代码是这样的：
model:
sets:
row/1..10/:c,d;
links(row,row):A;
endsets
data:
!c是给出条件中的Ma-Mj,d是给出条件中的Na-Nj;
c=2132 12998 351 777 463 291 2464 8412 9814 12984;
d=3173 5927 2931 2110 10000 11160 4050 995 6 143;
enddata
[OBJ]min=@sum(row(i):@abs(@sum(row(j):A(i,j))-c(i)))+@sum(row(j):@abs(@sum(row(i):A(i,j)-d(j))));
@for(links(i,j):A(i,j)>=0);
end

好的，我晚上回去看看，加我qq吧，我们一起探讨一下251881764

lingo的那个程序，貌似算不出来吧亲