Hello , I have been trying to plot the solution to a boundary layer problem using Maple via the following programme: >BVP:={0.1*diff(y(x),[x$2]) + (x^0.5)*diff(y(x),x) - y(x) = 0, y(0)=0,y(1)=exp(2)}; >J:=dsolve(BVP,y(x), type=numeric,method=bvp[midrich]); >plots[odeplot](J,[x,y(x)], 0..1,color=blue,style=line,symbol=circle, symbolsize=12); However,before I can write the instructions to plot, I get an error message saying: "Error,(in dsolve/numeric,bvp) unable to achieve requested accuracy of 0.1e-5 with maximum 128 point mesh (was able to get 0.17e-4), consider increasing "maxmesh" or using larger "abserr" ."

I'm running Maple 10.02, Mac OS X.4.3, and Java 1.5, and when I give Maple an inverse laplace transform with a number, it returns the correct answer. But when I give it the same thing but replace the number with an arbitrary constant (i.e., I replaced 4 with k) it doesn't seem to be able to figure it out. Any suggestions on this?

>

I just push the last version of my sheet "proba4maple" dedicated to discrete probabilistic computations on my SourceForge hive.

The NIST Engineering Statistics handbook. It's a comprehensive collection of information on the effective application of statistical methods. Submitted by Mariner.

As a newbie to maple, I was wondering how do I find the maximum of a equation under contraint. i can find some but for this one I can't ... and it's getting to me f=1000*x+7500*y under constraint: g=10000000-(20000*x^0.8)*y^1.2 and for f=(20000*x^0.8)*y^1.2 under constraint: g=125000-1000*x-7500*y Who do I find the maximum using MAple 10 !!!????

Let's all congratulate Carl Friedrich Gauss for winning our prestigious Greatest Mathematician contest. Our next poll brings us back down to earth and requests your input and comments on what facilities you work with most in Maple. T4.

Maplesoft just announced a patch available for download here.

I had been following a thread on comp.soft-sys.math.maple titled: Maple missing solutions?

http://groups.google.com/group/comp.soft-sys.math.maple/browse_thread/thread/0dc424e68954103e/9511b79909ecf139?hl=en#9511b79909ecf139

In, short the OP was missing solutions to a set of equations involving sin and cos.

A respondent suggested converting to exp, which seems to have worked very well in that case.

For no particular reason I decided to try this with another equation and got some

results which I cannot explain. As follows:

> S1:=[solve(sin(x)=sin(1/x))];

A student of mine recently asked what algorithm Maple uses to calculate Eigenvalues. So, I tried diving into Maple Procedures. For example: >showstat(`Eigenvalues`); LinearAlgebra:-Eigenvalues := proc(A) .. [error catchers ommitted] .. 17 LinearAlgebra:-LA_Main:-Eigenvalues(`if` ... end proc How should I interpret that last line? Can I get Maple to show me what that procedure does? So far I've no luck cracking it. Also, is this the best/easiest way of going about learning what types of algorithms Maple uses. I understand that some (/many/all) of them may be proprietary and the persons at Waterloo may not want the general public knowing them, but even an answer as simple as "Accelerated Newton-Jacobi algorithm" would be good.

Hi, I have plotted the following chart using maple 10 on windows xp > plot(((x^2)+x-4)/(x-2),x=-4..6,y=-2..10); X being -4, 2 Which command needs to be executed to find the local extremity? Thanks in advance for any help… Robert

I've posted a worksheet and some code for simplification with side relations. What makes this code interesting is that it properly handles rational expressions, using an algorithm is from my M.Sc. thesis. You can download the code from the Maple Application Center here, or from my personal webpage here. The code requires Maple 10.

Another hour, another possible problem. Consider the following Maple statements:

> f := table([index=m]);
f := table([index = m])
> g := copy(eval(f));
g := table([index = m])
> g[index] := one;
g[index] := one
> h := g;
h := g
> f[index], h[index];
m, one
> g := copy(eval(h));
g := table([index = one])
> g[index] := two;
g[index] := two
> l := g;
l := g
> f[index], h[index], l[index];
m, two, two

So, after the first copy, the second table h has value one for index, while f still has value m for index. Yet, after the second copy, not only the third table l has value two for index, but also the second table h now has value two for index, while I have explicitly used copy to make a copy of h. Does this indicate a problem with the copy function?

Hello,

a colleague just noticed the following (erroneous) behavior in Maple 10 :

> a := Psi(2, 2+2*I);
a := Psi(2, 2 + 2 I)
> printf("%+.6e\n", Re(a));
Error, (in fprintf) number expected for floating point format
> evalb(Im(a) <> 0);
false

The behavior is correct in Maple 9.5 :

> a := Psi(2, 2+2*I);
a := Psi(2, 2 + 2 I)
> printf("%+.6e\n", Re(a));
+3.902435e-02
> evalb(Im(a) <> 0);
true

-- Regards,

Franky.

Yet another question (after scouring the Maple documentation and coming up empty-handed): The following commands: with(student) intercept(y=sin(x), y=0) Result in this output: {y=0, x=0} How can I return the other intercepts with the x-axis? For instance, the first positive and first negative non-zero intercepts? (i.e. the "first" x-intercepts traveling left and right away from the y-axis?) Thanks as always, Bryan

Quick question: in maple , sqrt(4) returns 2. Why not plus/minus 2, since -2 is also a valid square root of 4? Thanks! Bryan