diff options
Diffstat (limited to 'src/algebra/laplace.spad.pamphlet')
| -rw-r--r-- | src/algebra/laplace.spad.pamphlet | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/src/algebra/laplace.spad.pamphlet b/src/algebra/laplace.spad.pamphlet index 9652f91d..28e2a1dd 100644 --- a/src/algebra/laplace.spad.pamphlet +++ b/src/algebra/laplace.spad.pamphlet @@ -116,7 +116,7 @@ LaplaceTransform(R, F): Exports == Implementation where -- returns [a, c, b] iff f = a * exp(c * t + b) -- and b and c do not involve t aexp(f, t) == - is?(f, "exp"::SE) => + is?(f, 'exp) => (v := isLinear(first argument(retract(f)@K),t)) case "failed" => "failed" [1, v.nconst, v.const] @@ -162,7 +162,7 @@ LaplaceTransform(R, F): Exports == Implementation where lapkernel(f, t, tt, ss) == (k := retractIfCan(f)@Union(K, "failed")) case "failed" => "failed" empty?(arg := argument(k::K)) => "failed" - is?(op := operator k, "%diff"::SE) => + is?(op := operator k, '%diff) => not( #arg = 3) => "failed" not(is?(arg.3, t)) => "failed" fint := eval(arg.1, arg.2, tt) @@ -170,11 +170,11 @@ LaplaceTransform(R, F): Exports == Implementation where ss * locallaplace(fint, t, tt, s, ss) - eval(fint, tt = 0) not (empty?(rest arg)) => "failed" member?(t, variables(a := first(arg) / tt)) => "failed" - is?(op := operator k, "Si"::SE) => atan(a / ss) / ss - is?(op, "Ci"::SE) => log((ss**2 + a**2) / a**2) / (2 * ss) - is?(op, "Ei"::SE) => log((ss + a) / a) / ss + is?(op := operator k, 'Si) => atan(a / ss) / ss + is?(op, 'Ci) => log((ss**2 + a**2) / a**2) / (2 * ss) + is?(op, 'Ei) => log((ss + a) / a) / ss if F has SpecialFunctionCategory then - is?(op, "log"::SE) => (digamma(1) - log(a) - log(ss)) / ss + is?(op, 'log) => (digamma(1) - log(a) - log(ss)) / ss "failed" -- Below we try to apply one of the texbook rules for computing |