tomleslie

13876 Reputation

20 Badges

15 years, 175 days

MaplePrimes Activity


These are replies submitted by tomleslie

  1. Open a new blank worksheet
  2. On the toolbar, access Tools->Options.
  3. Select the 'Display' tab
  4. Ensure that the first entry is set to 'Maple Notation', and the second entry is set to '2-D Math Notation'
  5. Click 'Apply to session' at the bottom left
  6. On the toolbar, access Tools->Options again
  7. Select the 'Interface' tab
  8. Ensure that the entry 'Default format for new worksheet' is set to 'Worksheet'
  9. Click 'Apply to session' at the bottom left
  10. Now open a new worksheet in the same session. There should be a prompt '>'.
  11. At the prompt, type 'version();' without the single quotes, but with the trailing semi-colon. This ought to return details of the precise Maple version you are running. If it returns something like

    User Interface: 1133417
                Kernel: 1133417
               Library: 1133417

    (precise numbers may be different), then you ought to be good to go - so try something simple like 1+1; at the next command prompt, (remember the semicolon)

Didn't set the E-mail/redsponse notification - I have now

@Jacobbn98 

Like I said earlier, the .bak (.bkp) gets deleted when Maple 'thinks' the document is successfully saved. And if it doesn't exist then the option will be 'greyed out'

When I examine the document you have linked above ie SSH_Binomialfordeling_Rødhåret.maple , it only contains one line of text - namely

"The resource you are looking for has been removed, had its name changed, or is temporarily unavailable."

So no-one is going to be able to recover anything from this.

I haven't analysed your code closely, (and it's past my bedtime), but perhaps you ought to consider the frequent occurrence of infinite, complex numbers in the 'soln' expression towards the end of you worksheet.

Plotting complex infinities can get a little tricky:-)

From my original post

##################################################

But

BartlettWindow, BlackmanWindow, HammingWindow, HannWindow, KaiserWindow

require that the option datatype=float[8] be set in the Array() constructor, which is used as input and always return an hfarray.

###################################################

My comment was why datatyp=float[8] is necessary for five windows, and not for the other eighteen.

Also why these 'odd' five windows return an 'hfarray', rather than a vector?

Your code is computing doCal(1), doCal(2),..doCal(5), rather than doCal(0), doCal(1000),....,doCal(5000), because the final loop contains

doCal(j)

rather than

doCal(H[j])

Corrected version attached

HelpFix3.mw

See the attached.

Note that the differences in the plots for H= [100,125,140,160]  are very minor. You have to zoom in a long way just to see them.

I added 0 and 1000 to the lit of H-values, ie H- [0, 100,125,140,160, 1000], just tso that I could visulaly check that different plots were actually being generated! Differences between these plots are still pretty minor - but they do exist!

HelpFix2.mw

Before I get into this some clarification is necessary

  1. "Downwards" seems to be excluded for some reason. Now "downwards" could occur as a result of a valid left or right turn (think spirals as a trivial case). But if a valid left or right turn would result in a 'downwards' step - then what is supposed to happen? Stay in current location?
  2. I also note that you have not allowed the 'reverse' option - is this just because it lands on a previous point and that means the walk is over? Some other reason? I note with some interest that the 2nd and 3rd 'walks' (mabe others, I haven't examined them all) in the example provided by Kitonum appear to have 'reverses', whihc probably shouldn't be happening.
  1. Define 'f' explicitly as a 'function'
  2. Probably a good idea to define 'x' as evalf(Int(..)), rather than int(), becuase it *ought to speed things up.

The attached works

corr1.mw

int(Heaviside(1-x^2-y^2), x=-infinity..infinity, y=-infinity..infinity);

*ought* to be the same as

int( Heaviside(1-r^2),  r=0..infinity, theta=0..2*Pi);

which is rather obviously 2*Pi - not Pi as you state. Interestingly, Maple actually returns 2*Pi when problem is posed in this form.
 

Some very similar to those for your previous post

  1. Again you supply code which cannot be executed, since there is no call to Erroran() with an appropriate set of arguments which will illustrate your problem. Try to get used to the idea that we generally need an executable example in order to debug - preferably one whihc fails with the error you are trying to fix!
  2. Don't even think about using parallel processing to speed things up  until you have have a working solution running in normal (ie non parallel) mode. You are just adding complexity upon complexity, making debug harder. Normal development would be to spend a lot of time getting a non-parallel implementation running and optimizing this so that it runs quickly. Then parallelize
  3. For this specific case, I suggest you investigate the value of the variable 'n' - it seems to crop up in a few places, (including the upper limit of an add() statement) but I can't find anywhere that its value is defined - of course I can't be sure, because I can't meaningfully execute your code - see point (1) above
  1. Since you do not supply a 'call' to the procedure HeatPC(), no one can actually run your code - thus the only possible debug is simply visual inspection: unlikely ever to be successful. Supply the call (with the arguments) which demonstrates the error you are seeing
  2. My (strictly visual) debug shows up a few things which could have been done better/more clearly, but nothing obviously terminal -eg you should probably be using add() rather than sum()
  3. The comment in your original "To make the code easier on the eye I want to split my main function u into several smaller functions. To that end I added the function a" Bad news - the function 'a' does not exist anywhere, so you have not "split my main function u into several smaller functions". In fact for the code you posted this remark is entirely meaningless

The shoot9 packages is not part of standard Maple, so you must have downloaded/installed it from somewhere. Others have succeeded with this - see, for example

http://www.mapleprimes.com/questions/214730-Solving-A-Nonlinear-System-Of-4-Equations

so it must be possible.

In the above link, the setup/calling sequence for commands from the shoot9 package seems to be identical to those which you are using, so I can only suggest that something is wrong with your installation of this package

You define four equations, assigned to E1, E2, E3, E4 - and then attempt to dsolve(), using equations assigned to EQ1, EQ2, EQ3, EQ4. I hope you can see the problem!!!!!

Your original code, ie

with(DEtools);
DEplot3d({diff(y(x), x) = y(x)-z(x), diff(z(x), x) = z(x)-2*y(x)}, {y(x), z(x)}, x = 0 .. 3, [[y(0) = 1.638, z(0) = 2.31]], y = 0 .. 2, z = -4 .. 4, scene = [x, z(x), y(x)], linecolor = COLOR(HUE, .5));

works in Maple 18,  User Interface: 991181,  Kernel: 991181, Library: 991181

works in Maple 2015.2, User Interface: 1097895, Kernel: 1097895, Library: 1097895

works in Maple 2016.1,  User Interface: 1133417 Kernel: 1133417, Library: 1133417

OP's second command sequence does produce  black/white plot in Maple 2015.2, and a prettiily colored line in Maple 2016.1. Can only conclude that this was a 'bug' in Maple 2015, which got fixed in Maple 2016(??) - so I'm not going to wrory about this one too much

Since OP has Maple 2016.1, the interesting question is why the original posted code does not produce an output graph - because on my 64-bit Win7, Maple 2016.1 produces a nice cyan(ish) curve.

One thing to try, is to attempt the intermediate step - as in just generate a numerical solution for the ODE system, and get some output values. On my machine

sol:= dsolve( { diff(y(x), x) = y(x)-z(x),
                       diff(z(x), x) = z(x)-2*y(x),
                       y(0) = 1.638,
                       z(0) = 2.31
                    },
                    numeric
                  );
sol(1);

returns

              proc(x_rkf45)  ...  end;
              [x = 1.,
               y(x) = HFloat(1.1065942587186357),
               z(x) = HFloat(1.492495225894533)]

What happens for you - is the solution procedure generated?

Are the numerical values generated?
 

 

 

First 133 134 135 136 137 138 139 Last Page 135 of 207