This works.

Thanks

Sorry for my late answer. I haven't been looking into this recently. But I'll let you know when I do.

Sorry for my late answer. I haven't been looking into this recently. But I'll let you know when I do.

actually you only changed the plotrange to 10^-10.

So why has maple problems calculating at k=0?

Since for deltag maple does not have problems.

For what (in particular here) do I need the transformation to another heun function of itsself?

Usually the function should be defined on this domain of definition.

Also I thought does maple perform possible analytic continuations on its own.

actually you only changed the plotrange to 10^-10.

So why has maple problems calculating at k=0?

Since for deltag maple does not have problems.

For what (in particular here) do I need the transformation to another heun function of itsself?

Usually the function should be defined on this domain of definition.

Also I thought does maple perform possible analytic continuations on its own.

what do you mean by wrapper?

what do you mean by wrapper?

Hello and thanks for your answer. haven't had time until now to look into it.

Well, using .ps gives indeed the best quality but I can only use it via plotsetup command and not via drop-down menu by rightclicking on the image. I can only find encapsulated postscript which is something else, right?

How can I setup legend and axis in plotsetup

plotsetup('png', plotoptions = "width=1000,height=1000,quality=200,boxed,noborder,portrait, legend=['curve']", plotoutput = "mypng.png")

doesnt work :-/

Hello and thanks for your answer. haven't had time until now to look into it.

Well, using .ps gives indeed the best quality but I can only use it via plotsetup command and not via drop-down menu by rightclicking on the image. I can only find encapsulated postscript which is something else, right?

How can I setup legend and axis in plotsetup

plotsetup('png', plotoptions = "width=1000,height=1000,quality=200,boxed,noborder,portrait, legend=['curve']", plotoutput = "mypng.png")

doesnt work :-/

@Carl Love jesus^^ it was just the semi-colon... why is this mandatory?

I've some last dumb question:

How exactly do I plot a pointplot of the following type

plots:-pointplot([seq(3+1/100, n = 1000 .. 5000)], l[1], axes = boxed, view = [3 .. 4, 0 .. 1])

Now this is correct for the first row...each element of l[1] shall be plotted @ 3.01

But in the same plot I want also the others

plots:-pointplot([[seq(3+1/100, n = 1000 .. 5000)], l[1]],[[seq(3+2/100, n = 1000 .. 5000)], l[2]],...,[[seq(3+100/100, n = 1000 .. 5000)], l[100]], axes = boxed, view = [3 .. 4, 0 .. 1])

where this unfortunately does not work

@Carl Love jesus^^ it was just the semi-colon... why is this mandatory?

I've some last dumb question:

How exactly do I plot a pointplot of the following type

plots:-pointplot([seq(3+1/100, n = 1000 .. 5000)], l[1], axes = boxed, view = [3 .. 4, 0 .. 1])

Now this is correct for the first row...each element of l[1] shall be plotted @ 3.01

But in the same plot I want also the others

plots:-pointplot([[seq(3+1/100, n = 1000 .. 5000)], l[1]],[[seq(3+2/100, n = 1000 .. 5000)], l[2]],...,[[seq(3+100/100, n = 1000 .. 5000)], l[100]], axes = boxed, view = [3 .. 4, 0 .. 1])

where this unfortunately does not work

thats what I tried but I keep getting:

for M from 1 to 100 do

forget(x)

 x:=proc(n)

option remember; eval(r*x(n-1)*(1-x(n-1)),r=(4*M)/(100))

end proc;

 x(0):=evalf(3/(5));

l[M]:=[seq(x(n),n=1..10)]

end:

Error, illegal use of an object as a name

thats what I tried but I keep getting:

for M from 1 to 100 do

forget(x)

 x:=proc(n)

option remember; eval(r*x(n-1)*(1-x(n-1)),r=(4*M)/(100))

end proc;

 x(0):=evalf(3/(5));

l[M]:=[seq(x(n),n=1..10)]

end:

Error, illegal use of an object as a name

I guess I need to forget the remember entries. But I cant find the correct command, since forget(x) seems to be wrong

Though Prebens Answer does work. This is simple and easy to handle^^

thx

Maybe some last question. If p ist not a constant but rather a function V(x) then this procedure does not account for the non-commutivity of Dx and V(x).