Should sympref Change the Behavior of rectangularPulse()?

The doc page for rectangularPulse states:

• If x = a or x = b and a <> b, then the rectangular pulse function equals 1/2.

Verify this statement:

rectangularPulse(0,1,[0 1])

ans = 1×2

0.5000 0.5000

syms t real
h(t) = heaviside(t) - heaviside(t-1);
h([0 1])

ans = 

But changing the sympref for HeavisideAtOrigin yields:

sympref('HeavisideAtOrigin',sym(1));
rectangularPulse(0,1,[0 1])

ans = 1×2

1 1

h([0 1])

ans = 

So rectangularPulse appears to not be the difference between heavisides. But the sympref('HeavisideAtOrigin') does impact the edge values of rectangularPulse. Is this expected behavior and a deficiency in the documentation?

[Walter Roberson]

By the way, you can see the actual code by using

regexprep(char(evalin(symengine, 'expose(rectangularPulse)')),'\\n','\n')

but it is not straight forward.

[Paul]

After running that command, it sure looks like that under "normal" conditions, the edge values are heaviside(0). So that would be an oversight in the documentation.