在最大值

Question

说我有以下功能：

f(x) := if x<=0 then 0 else if x<=1 then 1 else -1;

或任何其他分段定义的功能。

函数定义似乎有效：

(%i9) f(-11);
f(1/2);
f(2);
(%o7) 0
(%o8) 1
(%o9) -1

然而，此处不会评估集成。以某种方式可以在Maxima中获得积分值？如果没有，可以在maxima中以数字方式完成吗？

Answer 1

load(abs_integrate)获取abs_integrate包，使integrate能够处理unit_step。您必须根据unit_step编写分段函数。 E.g：

(%i1) load (abs_integrate) $
(%i2) e : unit_step(t) - 2*unit_step(t - 1) $
(%i3) integrate (e, t, a, b);
             abs(b) - b - 2 abs(b - 1) - abs(a) + a + 2 abs(a - 1)
(%o3)        -----------------------------------------------------
                                       2

对于数值积分，quad_qags（和其他quadpack函数）可以处理unit_step和if表达式。 quad_qags不需要abs_integrate。