MaplePrimes - answers and comments on Question, pdsolve:-plot changes the initial profile of the function at t=0...

Weakness

Markiyan Hirnyk — Mon, 28 Nov 2011 23:56:07 Z

Executing your code with Phi(r,0)=sqrt(r) and z:=1000 instead of z=1000, I have observed the following weakness:
> p := pds1:-plot(t = 0, numpoints = 100):
> plots[display](p);
Correct output
> p := pds1:-plot(t = 150, numpoints = 100):
> plots[display](p);
Correct output
> p := pds1:-plot(t = 0, numpoints = 100):#repeated exeqution
> plots[display](p);
Incorrect output
> p := pds1:-plot(t = 150, numpoints = 100):#repeated execution
> plots[display](p);
Correct output
After the restart; command and a repeated execution of the code the output is correct again.
See weakness.mw

Hey Markiyan,, Thanks,, but I dont

vish — Tue, 29 Nov 2011 00:14:14 Z

Hey Markiyan,,

Thanks,, but I dont quite understand.. I have not repeated the plot .. can you explain what "weakness" is please!

moreover my initial function is a complicated piecewise polynomial... is this causing problems. I thought the range I supply in pdsolve is tampering with the ranges in the piecewise polynomial.. but I have used different constant names for that ..

Any suggestion is much appreciated..

Thanks

Vish

Hi,, see below my code.. It's a bit long

vish — Tue, 29 Nov 2011 00:27:48 Z

Hi,, see below my code.. It's a bit long so you may want to look at the code around eq(30).. eq(41) .. Thanks for helping me .. If you need any more clarification on the code I would be happy to tell you..

Thanks

maple_code.mw

Opps sorry I have messed my code by writing comments...... between code..this is the right version

vish — Tue, 29 Nov 2011 00:30:55 Z

maple_code.mw

The spacestep option helps

Markiyan Hirnyk — Tue, 29 Nov 2011 02:14:24 Z

Put in the pdsolve command spacestep=0.5 and you will see the result in corrected_code3.mw

Pointlessness

Markiyan Hirnyk — Tue, 29 Nov 2011 00:19:16 Z

Can you kindly present your code in MaplePrimes by the green arrow in Menu?