Dr. David Harrington

5973 Reputation

21 Badges

19 years, 300 days
University of Victoria
Professor or university staff
Victoria, British Columbia, Canada

Social Networks and Content at Maplesoft.com

Maple Application Center
I am a professor of chemistry at the University of Victoria, BC, Canada, where my research areas are electrochemistry and surface science. I have been a user of Maple since about 1990.

MaplePrimes Activity

These are replies submitted by dharr

@Carl Love  To get the same output as Ronan I needed to switch jm and h in the updates line to become:

(im, istep, h, jm):= Scale2~([im, istep, h, jm], [-1$3, 1])[]

Not sure if this affects speed because I know modular was updated after the version I have (2017)

@Ronan I usually use Compiler:-Compile when I'm debugging, since there are many limitations about what Maple code can be compiled, but once I know it does compile then you can just add option autocompile after the proc() and the compiling will be automatic on first execution (and so the procedure can be used in a module or package, and the dll will be generated automatically). 

@student_1 Yes, fracdiff(v(t), t,0) gives t but fracdiff(v(t), t,0.) gives 0. The solution is to use seq(fracdiff(v(t), t,alpha/2),alpha=-2..3)

@MajorMax So this is just the plot command



@nm There are many Computer Modern fonts on my system (I use Scientific Workplace) and I don't know which are which - I just chose cmr12 and cmr16 at random, but cmbxti10 looks more like what you want. The above x may just be the standard one. 

P1 := plot(x^2):
P2 := plot(x^2,labels=["x","y"],axesfont=[cmr12,roman,16],
plots:-display(Array([ P1, P2 ]));

gives (snapshot from export to .pdf


@Earl  You/re welcome. As you probably realize, if the control point just goes  x[1] = pts[2][1]..pts[1][1] you will just get the centre part with ellipses in the quadrilateral, which you asked about in a separate question..

@mmcdara  Not at all - I expected someone could help out if that wasn't enough.  It seemed a question about why dsolve,series didn't work, and as is often the case the OP's Maple background wasn't very clear..

I'm out of time these days, but had a quick look. I didn't check the eqns in detail, so can't comment on the method, but didn't see anything obviously wrong, except that "line" is in the plottools package, which wasn't loaded. I just tried regular solve, without the parametric option, but still lost connection with the kernel, which shouldn't happen. So my two tips would be:

Change the floating point numbers like 0.5 to exact 1/2 etc, since solve is for symbolic solutions (maybe solve would do this internally, but it is better to be sure and it is anyways good practice).

Rearrange the tangent equations to remove the denominators so a/b=c changes to a=b*c. Maple's solve is very good for polynomial systems, so that may work better.

@Axel Vogt  Some would say it is zero at x = infinity (complex infinity).  I think this is fairly standard when talking about poles and zeros.

@ecterrab I agree with @nm . The journals I use request \mathrm inside math mode when typesetting equations. I guess it would be \textrm outside.

Here's a variation where one of the coefficients becomes the parameter to vary. Conceptually simpler perhaps, but what the variation leads to and the flexibility is a bit more obscure.


I just copied the output of DirectSearch to bypass that - then it seemed Pen1 should have been Pencil at the point highlighted in red. Then your ellipses work as you describe. But I don't understand what else you want to do, so I'll leave that to someone else.


Note that LegndreP and JacobiP are not in the orthopoly package (but P in that package is related).

@max125  If you want the nearest integer then ilog2(8) works, so it is interesting that log2(8) doesn't.

@max125 log10() works as a short form.

First 40 41 42 43 44 45 46 Last Page 42 of 62