<< linalg_gausspivotal Linalg linalg_hbandL >>

Linalg >> Linalg > linalg_gausspivotalnaive

linalg_gausspivotalnaive

Computes the solution of a linear equation with Gauss and row pivoting.

Calling Sequence

x = linalg_gausspivotalnaive ( A , b , verbose )

Parameters

A :

a n-by-n matrix of doubles

b :

a n-by-1 matrix of doubles

verbose :

a 1-by-1 matrix of boolean, set to true to display intermediate messages

x :

a n-by-1 matrix of doubles

Description

Returns the solution of Ax = b with a Gauss elimination and row exchanges. This algorithm might be sensitive to inaccuracies if A near singular. There is no solution is A is singular. If A is not ill-conditionned, the algorithm is accurate. Uses a naive algorithm with all loops expanded (no vectorization).

Examples

A = [
3 17 10
2 4 -2
6 18 -12
];
b = [
67
4
6
];
x = linalg_gausspivotalnaive ( A , b , %f )
xe = [1 2 3]'

Authors

<< linalg_gausspivotal Linalg linalg_hbandL >>