Solves the symmetric indefinite linear system Ax = b.

Syntax

x = LSolveSPD(A,b)

Arguments

Name

Description

A

A square, symmetric positive definite matrix of real values.

b

The right-hand side column vector or a matrix.  If it is a matrix, each column is considered a separate vector and the system is solved separately for each, resulting in multiple solutions. The number of rows must be equal to the row/column size of the square matrix A.

Output

Name

Description

x

The solution(s) to the system(s).  It will have the same dimension as the right-hand side argument b.  Each column will be a solution corresponding to the equivalent column in b.

Example

Solve the symmetric positive definitive linear system Ax=b.

Syntax

A = [40,2,3;2,40,4;3,4,40];//symmetric positive definite

b = [1;11;21]

x = LSolveSPD(A,b)

Result

A column vector of the solution.

x =   -0.02

      0.23

      0.50

See Also:

LSolveB

LSolveT

LSolveSPDB