diff options
| author | dos-reis <gdr@axiomatics.org> | 2008-04-03 04:23:42 +0000 |
|---|---|---|
| committer | dos-reis <gdr@axiomatics.org> | 2008-04-03 04:23:42 +0000 |
| commit | 001e19b08ba7fb1b9e6f6bdb44a82ba3db3fc532 (patch) | |
| tree | da9e2fe5d81ff4cd7709d12e44b8c3e348b8a8e3 /src/algebra/odeef.spad.pamphlet | |
| parent | a7bab9a6c2070d05e2dbd256ce455079c8ced385 (diff) | |
| download | open-axiom-001e19b08ba7fb1b9e6f6bdb44a82ba3db3fc532.tar.gz | |
Replace `^=' with `~='.
Diffstat (limited to 'src/algebra/odeef.spad.pamphlet')
| -rw-r--r-- | src/algebra/odeef.spad.pamphlet | 16 |
1 files changed, 8 insertions, 8 deletions
diff --git a/src/algebra/odeef.spad.pamphlet b/src/algebra/odeef.spad.pamphlet index 734344a9..e33f827e 100644 --- a/src/algebra/odeef.spad.pamphlet +++ b/src/algebra/odeef.spad.pamphlet @@ -56,7 +56,7 @@ ReductionOfOrder(F, L): Exports == Impl where ithcoef(eq, i, s) == ans:F := 0 - while eq ^= 0 repeat + while eq ~= 0 repeat j := degree eq ans := ans + localbinom(j, i) * locals(s,j,i) * leadingCoefficient eq eq := reductum eq @@ -213,7 +213,7 @@ ElementaryFunctionLODESolver(R, F, L): Exports == Implementation where up := univariate(p, k := first l) l := rest l ans:P2 := 0 - while up ^= 0 repeat + while up ~= 0 repeat ans := ans + monomial(upmp(leadingCoefficient up, l), k, degree up) up := reductum up ans @@ -284,7 +284,7 @@ ElementaryFunctionLODESolver(R, F, L): Exports == Implementation where ulodo(eq, k) == op:LQ := 0 - while eq ^= 0 repeat + while eq ~= 0 repeat op := op + monomial(univariate(leadingCoefficient eq, k), degree eq) eq := reductum eq op @@ -294,7 +294,7 @@ ElementaryFunctionLODESolver(R, F, L): Exports == Implementation where op := ulodo(eq, k) empty? remove_!(k, varselect(kernels g, x)) => -- i.e. rhs is rational rc := ratDsolve(op, univariate(g, k)) - rc.particular case "failed" => -- this implies g ^= 0 + rc.particular case "failed" => -- this implies g ~= 0 doVarParams(eq, g, homosolve(eq, op, rc.basis, k, x), x) [multivariate(rc.particular::RF, k), homosolve(eq, op, rc.basis, k, x)] doVarParams(eq, g, homosolve(eq, op, ratDsolve(op, 0).basis, k, x), x) @@ -316,7 +316,7 @@ ElementaryFunctionLODESolver(R, F, L): Exports == Implementation where localmap(f, op) == ans:L := 0 - while op ^= 0 repeat + while op ~= 0 repeat ans := ans + monomial(f leadingCoefficient op, degree op) op := reductum op ans @@ -475,7 +475,7 @@ ElementaryFunctionODESolver(R, F): Exports == Implementation where p := rec.left Lx := LinearOrdinaryDifferentialOperator(F, diff x) op:Lx := 0 - while p ^= 0 repeat + while p ~= 0 repeat op := op + monomial(leadingCoefficient p, degree p) p := reductum p solve(op, rec.right, x, a, y0)$ElementaryFunctionLODESolver(R, F, Lx) @@ -495,7 +495,7 @@ ElementaryFunctionODESolver(R, F): Exports == Implementation where p := rec.left Lx := LinearOrdinaryDifferentialOperator(F, diff x) op:Lx := 0 - while p ^= 0 repeat + while p ~= 0 repeat op := op + monomial(leadingCoefficient p, degree p) p := reductum p (uuu := solve(op, rec.right, x)$ElementaryFunctionLODESolver(R, F, Lx)) @@ -504,7 +504,7 @@ ElementaryFunctionODESolver(R, F): Exports == Implementation where -- returns [M, v] s.t. the equations are D x = M x + v parseSYS(eqs, ly, x) == - (n := #eqs) ^= #ly => "failed" + (n := #eqs) ~= #ly => "failed" m:M := new(n, n, 0) v:V := new(n, 0) xx := x::F |