# Revision history [back]

This has to do with the implementation of piecewise, I think.

I can get a numerical integral to work by doing the following:

phi(x)=(heaviside(x+1)-heaviside(x-1))*(1-abs(x))*(1-abs(x))*(1+2*abs(x))
numerical_integral(phi(x/h-1),0,pi)


Another option would be to stop the integral at the point at which the function is zero.

phi2(x)=(1-abs(x))*(1-abs(x))*(1+2*abs(x))
integral(phi2(x/h-1),x,0,2*pi/n)