Ask Your Question

Revision history [back]

click to hide/show revision 1
initial version

Thank-you for the quick answer.

frac() gives a function which is periodic on positive reals. I just adapted your solution to get it periodic on \RR.

f = lambda x: 1 if (x - RR(x).floor()) < 1/2 else 0

I would have liked to be able to define a symbolic function though. Is it doable?

I was also trying to get a function whose plot is correct (without asking it to be pointwise) as it is the case for Piecewise().

Below are some of the things (not chronologically ordred though) I had tried (just for completness as I guess they are just full of beginer's classical mistakes).


which does not evaluate


returns ValueError: cannot convert float NaN to integer.


returns ValueError: Value not defined outside of domain.


also returns ValueError: Value not defined outside of domain.

I also tried to redefine unit_step:

def echelon_unite(x):
    if x<0:
    return hres

problem integral(echelon_unite(x),x,-10,3) returns 13

numerical integral returns a coherent result.

Other tentative with incoherent result (still with integrate and not numerical_integral)

  def Periodisation_int(f,a,b):
        x = var('x')
        h0(x) = (x-b)-(b-a)*floor((x-b)/(b-a))
        hres = compose(f,h0)
        return hres

sage: g=Periodisation_int(sin,0,1)
sage: integrate(g(x),x,0,2)
-cos(1) + 2
sage: integrate(g(x),x,0,1)
-cos(1) + 1
sage: integrate(g(x),x,1,2)
-1/2*cos(1) + 1

My guess is I was using integrate on inappropriate objects. I would still like to know how to define corresponding symbolic function (if it is possible).

Thanks again,

best regards.