1
edit
Line 664: | Line 664: | ||
function x = mldivide (A, b) | function x = mldivide (A, b) | ||
m = rows(A); | |||
[Q, R, E] = qr(A); | [Q, R, E] = qr(A); | ||
x = [A \ b, E(:, 1:m) * (R(:, 1:m) \ (Q' * b))] | x = [A \ b, E(:, 1:m) * (R(:, 1:m) \ (Q' * b))] |
edit