myhilb.m

来自「在matlab中」· M 代码 · 共 22 行

function A=myhilb(n, m)
%MYHILB  本函数用来演示MATLAB语言的函数编写方法。
%   A=MYHILB(N, M) 将产生一个N行M列的Hilbert矩阵A；
%   A=MYHILB(N) 将产生一个NxN的方Hilbert阵A
%   MYHILB(N,M) 调用格式只显示NxM的Hilbert矩阵，但不返回任何矩阵。
%  
%See also: HILB.

%  Designed by Professor Dingyu XUE, Northeastern University, PRC
%     5 April, 1995, Last modified by DYX at 30 July, 2001
if nargout>1, error('Too many output arguments.'); end
if nargin==1, m=n; 
elseif nargin==0 | nargin>2
   error('Wrong number of input arguments.'); 
end
A1=zeros(n,m);
for i=1: n
   for j=1:m
      A1(i,j)=1/(i+j-1);
end, end
if nargout==1, A=A1; elseif nargout==0, disp(A1); end