% 行列と連立一次方程式求解のサンプル

n = 100;
A = hilb(n);        % Hilbert行列 (A_ij = 1/(i+j-1))
x = ones(n,1);      % サンプルの解
b = A*x;            % 問題を作る Ax = b

tic;                % 時間計測開始
xx = A\b;           % Matlab の，「Ax=b」を解くルーチンを呼び出す簡易コマンド
toc;                % 時間計測終了

norm(x-xx);         % 果たして解の精度はいかに

plot(x);            % 図示
hold on;            % 同じ図面に複数のプロットをさせるための呪文
plot(xx);