I think I've got a vague idea what happens here... The same thing also
happens when
I define the max() function myself:
sage: def my_max(x,y):
sage: if(x>y): return x
sage: else: return y
sage: fermi2(x,y,d,L) = 1 - 1/( exp( ( my_max(abs(x),abs(y))-L) /d) +
1)
sage: fermi2
(x, y, d, L)
On Mar 12, 5:47 pm, Harald Schilly wrote:
> I'm not sure but there might be a bug or problem evaluating the
> expression. Anyways, going the "pure" python way works:
>
> sage: def fermi(x,y,d,L): return 1 - 1/( exp( ( max(abs(x),abs(y))-
> L) /d) + 1)
>
> sage: plot3d(lambda x,y : fermi(x,y,0.
On Mar 12, 5:39 pm, stefan wrote:
> I have encountered a somewhat strange problem ...
Ah, I got it, the max function is evaluated during the definition of
the function, look:
sage: fermi(x,y,d,L) = 1 - 1/( exp( ( max(abs(x),abs(y))-L) /d) + 1)
sage: fermi
(x, y, d, L) |--> -1/(e^(-(L - abs(x))/
On Mar 12, 5:39 pm, stefan wrote:
> Basically, I am trying to plot this expression:
>
> #sage> fermi(x,y,d,L) = 1 - 1/( exp( ( max(abs(x),abs(y))-L) /d) + 1)
I'm not sure but there might be a bug or problem evaluating the
expression. Anyways, going the "pure" python way works:
sage: def fermi(
