## Documentation |

Solve almost block-diagonal linear system

`x = slvblk(blokmat,b) x = slvblk(blockmat,b,w) `

`x = slvblk(blokmat,b) ` returns
the solution (if any) of the linear system `Ax = b`,
with the matrix `A` stored in `blokmat` in
the spline almost block-diagonal form. At present, only the command `spcol` provides
such a description, of the matrix whose typical entry is the value
of some derivative (including the 0th derivative, i.e., the value)
of a B-spline at some site. If the linear system is overdetermined
(i.e., has more equations than unknowns but is of full rank), then
the least-squares solution is returned.

The right side `b` may contain several columns,
and is expected to contain as many rows as there are rows in the matrix
described by `blokmat`.

`x = slvblk(blockmat,b,w) `
returns the vector `x` that minimizes the *weighted* sum
Σ_{j}*w*(*j*)((*Ax* – *b*)(*j*))^{2}.

`sp=spmak(knots,slvblk(spcol(knots,k,x,1),y.'))` provides
in `sp` the B-form of the spline s of order `k` with
knot sequence `knots` that matches the given data `(x,y)`,
i.e., for which s`(x)` equals `y`.

Was this topic helpful?