s13rout3.m

来自「Industrial Mathematics」· M 代码 · 共 14 行

% Routine 13.3
N = 6;
Id = eye(N-1);
T = ones(N-1);
T = tril(T,1);
T = triu(T,-1);
T = T - Id;
G = 2*Id/(3*N) + T/(6*N);
S = 2*N*Id - N*T;
b = ones(N-1,1)/N;
%
c= (G-S)\b;
%