diff options
Diffstat (limited to 'src/algebra/odealg.spad.pamphlet')
| -rw-r--r-- | src/algebra/odealg.spad.pamphlet | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/src/algebra/odealg.spad.pamphlet b/src/algebra/odealg.spad.pamphlet index 653dece6..1a3d355d 100644 --- a/src/algebra/odealg.spad.pamphlet +++ b/src/algebra/odealg.spad.pamphlet @@ -152,7 +152,7 @@ SystemODESolver(F, LO): Exports == Implementation where r := r - 1 c := c - 1 firstnonzero(m, r) ~= c => error "backsolve: undetermined system" - degree(eq := m(r, c)) > 0 => error "backsolve: pivot of order > 0" + positive? degree(eq := m(r, c)) => error "backsolve: pivot of order > 0" a := leadingCoefficient(eq)::F if part? then part(r + offset) := (v(r + offset) - applyLodo(m, r, part, n)) / a |