Question: Solve linear AX=F

Hi,

I need to solve the linear system  AX=F  without using any  Linearalgebra package. 

 

A is a tridiagonal matrix ( band matrix);

X is vector of size n.

F is the right hand side of my equation AX=F.

F: given, A given, Find X???

restart;

A := n-> Matrix(n,scan=band[1,1],[[seq(b[i],i=2..n)], [seq(a[i],i=1..n)],
                [seq(c[i],i=1..n-1)]] ):
X:=n->[seq(x[i],i=1..n)]:
F:=n->[seq(d[i],i=1..n)]:

 

Thanks

Please Wait...