So I have an integral that computes perfectly in wolfram alpha but not in maple...

I will post it here

int(1/((4.532055545*10^9/f^4.14-2.311250000*10^5/f^2+(111*(1-0.2163331531e-4*f^2+2.340001656*10^(-10)*f^4))/(1+0.1081665766e-4*f^2)))*(6*10^(-21)*abs(1/f^(4/3)))^2, [f = 50 .. 1500])

the answer should be 3.05364*10^-46

If you try that exact line of code in maple, it will not compute (is stuck on evaluating)


Best Regards to all,
Zeus

I'm displaying a series of point plots as an animation, and would like to update a displayed parameter as well.  I have a nested list L[t] where there's a set of points for each t, and for each t there's also a numerical value M that I'd like to display.  (In my real problem, L[t] is the number of particles in each of several states, and M is the rms deviation from am algebraic probability distribution.)  

The closest I've gotten (for a simple L and M) is the following, but it displays all of the M values in the legend at once:

with(plots); with(Statistics);

L := [[1, 2, 3, 4, 5, 6], [2, 4, 6, 8, 10, 12]];

M := [1, 2];

display([seq(PointPlot(L[t], legend = M[t]), t = 1 .. 2)], insequence = true)

I don't need this to be in the legend.  Is there a way to display only the current value of M for each t?  Thanks very much.

How to prove or disprove the flatness of the surface x = (u-v)^2, y =  u^2-3*v^2, z = (1/2)*v*(u-2*v), where u and v are real-valued parameters? Here is my try:

I think Gaussian curvature should be used to this end. Dr. Robert J. Lopez is my hope.

Download flat.mw

I've been trying to make a smooth plot of some ODEs. It should show C rapidly increasing at the innitiation, until they get into a quassi steady state,  and then all three variables increase much slower. This should look like a roughly straight line that elbows sharply into a smooth curve.

Any attempt to DEplot3d it i've made either just shows the time before the quassi steady state is reached, so shows the straight line; or smooths that time together with the next period, making the straight line look like a part of the smooth curve.

Model := [diff(B[1](t), t) = k[a1]*C(t)*(R-B[1](t)-B[2](t))-k[d1]*B[1](t), diff(B[2](t), t) = k[a2]*C(t)*(R-B[1](t)-B[2](t))-k[d2]*B[2](t), diff(C(t), t) = (-(k[a1]+k[a2])*C(t)*(R-B[1](t)-B[2](t))+k[d1]*B[1]+k[d2]*B[2](t)+k[m]*((I)(t)-C(t)))/h];
DissMod := subs((I)(t) = 0, Model);
AssMod := subs((I)(t) = C[T], Model);


Pars := [k[a1] = 6*10^(-4), k[d1] = 7*10^(-3), k[a2] = 5*10^(-4), k[d2] = 10^(-2), R = .5, k[m] = 10^(-4), C[T] = 100, h = 10^(-6)]

StateSol := DEplot3d(subs(Pars, AssMod), [B[1](t), B[2](t), C(t)], t = 0 .. 1000, number = 3, B[1] = 0 .. .5, B[2] = 0 .. .5, [[B[1](0) = 0, B[2](0) = 0, C(0) = 0]], scene = [B[1](t), B[2](t), C(t)], maxstep = .1, maxfun = 0, method = l)

Is there a way to display the list of coordinates for implicitplot and plot, rather than the displaying image itself?

In Maple Classic this happens automatically when the assignment operator is used  ( := ) together with a semi-colon, eg.   A:=implicitplot(F, x=-1..1, y=-1..1);

But is there a way to do this in the non-classic version?  Thanks very much.

- Graham

Hello fellow maple struglers!

This is my first question, so hopefully I did it in the correct manner. My problem (I think) is of a simple nature, but I can't seem to get it to work. The equation I have is the following:

All is nice and well and as you all can see it depends on two variables. Now what I want to do is to express the function in a shape where it depends on the ratio t/d only. Simple calculus shows that for this equation it is rather easy, where the result is:

(62,8*t/d+0,0945)/pi

Now I would like to know whether it is possible to automate this in maple, such that I can plot Wr_over_p as a function of t/d. I know I can simply use the form above and make the plot, but further down in my maple file I have an equation which is much more complex, and I want to plot thatone as a function of t/d as well.

I hope I have made my problem clear to some of you guys, but if you require some further explaination than please let me know. Thanks in advance for your help, I hope there is someone who can help me with this!

I have two deformed planes, that i would like to draw with 3dplot, as well as drawing a curve marking their intersection.

the curves are given by the expressions:

C = -(k[d2]*B[2]+I*k[m]+k[d1]*B[1])/((B[1]+B[2]-R)*k[a1]+(B[1]+B[2]-R)*k[a2]-k[m]),

C = k[d1]*B[1]/(k[a1]*(R-B[1]-B[2]))

evaluated at

Pars := [k[a1] = 6*10^(-4), k[d1] = 7*10^(-3), k[a2] = 5*10^(-4), k[d2] = 10^(-2), R = .5, k[m] = 10^(-4), C[T] = 100, h = 10^(-6)]

with the variables B[1],B[2] and C within the bounds [0..0.5],[0..0.5],[0..100].

My method was to try and use solve to find a formula for the intersection curve- but i couldn't get 3dplot to plot it!

Consider the following code snippet:

with(DifferentialGeometry):
DGsetup([x],M);
RemoveFrame(M);

By itself it is, of course, pretty meaningless, but that is not the point. My problem is that after Maple has executed it, it will annoyingly continue using M> for each new execution group that is inserted. Having removed the frame with the command RemoveFrame, I would have expected that not to be the case. How can the frame be truly removed?

hi. i have a system of ODEs. is there any answer ? rho is constant !

Download dsolve.mw

Hi, I have this procedure (from Maple 5) but I am using it in Maple 15. My problem is this program can not run. I think there is some commond incorrect, but not sure which ones. Please help me in this problem. Thanks a lot.

cocycle.mw

Download cocycle.m

in Mathematica, there is the option called BoxRatios

"is an option for Graphics3D that gives the ratios of side lengths for the bounding box of the three‐dimensional picture."

It is sort-of like aspect ratio, but for 3D. It is set by default so make 3D plot looks "nice". I can't seem to find equivalent Maple option. The closest is the option "s=" for plot3d, but this just turns of/on "constrained scaling" and does not allow one to modify the "BoxRatios"

Let me give an example. Here is 3D plot in Mathematica and the same in Maple. I'd like to get the Maple 3D to look similar to Mathematica 3D in terms of the "aspect ratio". Maple on the z-axis is using the same size as in the x and y axis, and even though this is realistic, it does not make the plot as nice. I want to change this ratio.

T0[x_, y_, m_] :=20/Pi Sum[ (-1)^(n + 1)/n Exp[- (n Pi/10) y] Sin[ (n Pi/10) x], {n,1, m}]
Plot3D[T0[x, y, 70], {x, 0, 10}, {y, 0, 10}, PlotRange -> All,  AxesLabel -> {x, y, z}]

In Maple:

T0:= (x,y,m)-> 20/Pi*sum( (-1)^(n+1)/n*exp(- n*Pi/10*y)*sin(n*Pi/10*x),n=1..m);
plot3d(T0(x,y,50),x=0..10,y=0..10,scaling=unconstrained);

So Maple is using 1:1:1 box ratio. Mathematica default is 1:1:0.4, and I wanted to see if I can change Maple to be the same.

I get same plot in Maple using scaling=unconstrained or scaling=constrained. So this option is not very useful for what I want.

Is there a way to change the "BoxRatios" as defined above in Maple? There must be, right? Do I need to use different package?

Dear All

In following I tried to find symmetries of certain partial differential equation taken from paper "Group classification and exact solutions of generalized modified Boussinesq equation". But the determining equations are not matching with equations obtained in paper.

Download [1116]_Symmetries_Determination.mw[1116]_Group_classification_and_exact_solutions_of_generalized_modified_Boussinesq_equation.pdf

Regards

Still a little unclear what this error means tho

Download dispatchTOshillCORE.mw

yep the errors recieved using some packages are very very specific for maple, for example, the one i got today using the ODE package was profoundly helpful:

Error, (in ODEtools/info) unable to handle derivatives as {diff(1/(ln(X)-Psi(1-f(X))-Psi(f(X))), [`$`(X, n-k[1]-k[2])]), diff(1/f(X), [`$`(X, n-k[1]-k[2])]), diff(Psi(1-f(X)), [`$`(X, k[2])]), diff(Psi(f(X)), [`$`(X, k[2])]), diff(f(X), [`$`(X, k[1])])} while solving w.r.t f(X)