Carl Love

Carl Love

28100 Reputation

25 Badges

13 years, 107 days
Himself
Wayland, Massachusetts, United States
My name was formerly Carl Devore.

MaplePrimes Activity


These are replies submitted by Carl Love

@Preben Alsholm My square brackets after D were unintentional. I'm also surprised that they worked. Thanks for spotting that. I'll fix the original.

@Axel Vogt Continuity of the integrand is obvious over a nonnegative interval. But Maple refuses to do the definite integral even when specific nonnegative limits of integration are supplied. Here's a simplified example:

Maple does the indefinite integral:

int(t^(5/2)*exp(t^(5/2)), t);

Maple does not do the definite integral:

int(t^(5/2)*exp(t^(5/2)), t= 0..1);

@dohashi I agree with itsme that Grid:-Map should just go ahead and copy the full memory state of the main kernel.

Regarding that fourth paragraph of ?Grid,Map : Could you show an example where Map is not "called outside the grid, on a non-compute node, not as part of a parallel computation"? I.e., an example where it is called as it "is designed to be called on all the nodes of a parallel computation." Would the last example at ?Grid,Map , where Launch is used, be an example of what the fourth paragraph refers to as "robust" use?

@itsme: Does your wrapped and Mapped procedure use global variables? I wonder if that is the issue. If it is, then maybe you could put the procedure in a module.

@mehdi jafari Numeric solution of PDEs without a time component is not currently implemented.

@Edinburgh Please explain what you are trying to plot. Is it simply the bezier curve with control points P[0]...P[5]? Why are you using complex numbers?

@Ahm3d Are you sure that the range of the integrals is supposed to be n..(n+1)*h? That seems strange to me. Are you sure that it is not n*h..(n+1)*h?

@Ahm3d Sorry about that. I had eval(F, n*h). That was supposed to be eval(F, t= n*h). I corrected the Answer. Please try again.

Since there is no actual error in the integration after you do the restart, what you say about the numerator is irrelevant.

You might be better off using a procedure-based numerical integration for the Optimization. I'll try it both ways in a few hours. Can I assume that is positive integer and h is positive?

I think that everyone is having that problem as a result of the MaplePrimes upgrade on Wednesday April 16.

@geischtli The main tool to use is plottools:-scale. This operates on a frame after it has been created by plot3d. In the code below, I take the fourth and final frame from your animation, l(4), and put it on the same scale (same axes) as the first frame. Then I display these in a simple two-frame animation:

x1:= -2..1:              y1:= -1.5..1.5:     #first frame axes
x4:= -0.82..-0.7:   y4:= -0.2..-0.08: #last frame axes
L(4):= plottools:-scale(
     l(4), #first frame from your code
     (op(2,x1)-op(1,x1))/(op(2,x4)-op(1,x4)), #scale factor for x
     (op(2,y1)-op(1,y1))/(op(2,y4)-op(1,y4)), #scale factor for y
     1, #don't change z
     [(op(1,x4)+op(2,x4))/2, (op(1,y4)+op(2,y4))/2, 0] #central point of last frame
):

plots:-display([l(1), L(4)], insequence = true);

Unfortunately, the tickmarks are the same on each frame. I don't think that this can be changed. The ordinary tickmarks option will not correct this: The tickmarks of the last frame are applied to all frames once it's animated.

@sarra It's a periodic function, so does the actual interval matter as long as the width is correct? Here's a plot of Markiyan's Fourier Series for the convolution:

Pi/2+Sum(2*((-1)^i-1)*cos(i*x)/i^2/Pi, i= 1..infinity);

plot(Re(value(%)), x= -2*Pi..2*Pi);

(Can't upload plots now. But you can generate it from the code above. It looks like the correct convolution.)

 

Are you trying to plot a surface or a curve? complexplot3d is a command that plots a surface. Look at ?complexplot3d . There are four calling sequences. Which one were you trying to emulate?

@itsme Looking at the first few lines of showstat(Grid:-Map) and executing the code below shows that it behaves differently depending on whether the first argument is a named procedure or an anonymous procedure. This is because of the "last name evaluation" property of named procedures. (The issue is not really the use of the arrow, although the most common anonymous procedures are arrow expressions.)

debug(Grid:-Map):
Grid:-Map(x-> x^2, [1,2,3]);
f:= x-> x^2:
Grid:-Map(f, [1,2,3]);

Certainly this should be considered a bug since the help page ?Grid,Map shows examples with anonymous procedures.

There is a minor anomaly with the pizza problem. I wonder if MapleTA can resolve it. I repeat the problem here for convenience.

A group of 11 male and 13 female students is planning to go out for pizza. If 82% of the male students go and 62% of the female students go, find the probability that a random student who goes out for pizza is female.

The "expected" answer, judging from the multiple-choice answers, is

(.62*13)/(.62*13+.82*11) = 0.4719.

However, the realistic answer uses a whole number of males (.82*11 = 9) and females (.62*13 = 8) and is thus

8/(8+9) = 0.4706,

which is not one of the choices. Can MapleTA automatically adjust the answers to reflect this? It would just involve using Maple's round function. This anomaly also illustrates the inherent problem of requiring more significant digits in the answer (4 sig. digits) than were provided in the problem data (2 sig. digits).

I also wonder if MapleTA was used to include the drawing of the pizza to the right of the problem.

You may have a problem displaying 3d plots. Try this:

plot3d(0, x= 0..1, y= 0..1, axes= boxed);

@Preben Alsholm 

This avoids the step of having to substitute z = exp(2*I*x).

S:= -ln(2)+Sum((-1)^(n+1)/n*cos(2*n*x), n= 1..infinity);

expand(simplify(evalc(Re(value(S))))) assuming cos(x) > 0;

 

First 550 551 552 553 554 555 556 Last Page 552 of 709