pentaSolve {AdvDif4}

Solving 'Ax=b' system

Description

This software solves an 'Ax=b' pentadiagonal system using a direct method. Variables a1, a2, a3, a4, a5 are matrix A diags and b is the vector.

Usage

pentaSolve(a1,a2,a3,a4,a5,b)

Arguments

a1

A vector

a2

A vector

a3

A vector

a4

A vector

a5

A vector

b

A vector

Value

A vector with the x value

Examples

#
# Solve a 'Ax=b' easy sample
#
a1<-c(1)
a2<-c(2,2)
a3<-c(7,7,7)
a4<-c(2,2)
a5<-c(1)
b<-c(11,12,13)
pentaSolve(a1,a2,a3,a4,a5,b)