MaplePrimes Questions

When I type the command

sum(2^n, n=0..N)

Maple calculates it correctly. However, if I then go back to that expression and click on it to change the value of N, Maple tries to execute something before I hit enter. It won't stop unless I press the Interrupt button, and after that I have to insert some other calculation, e.g. 1+1, for it to work once again.

This happens even on completely new files and for any explicit values of N.

Why is this happening?

Hello! I am trying to plot a parabola given by the equation x=y^2. I used the following code

plot(y^2, y = -1 .. 1, color = blue, thickness = 3, title = "T1.1 x=y^2 [ECE]")

but my plot looks like this

I want it to have the x-axis horizontally displayed and the y-axis to be vertically displayed. Something like this

File:X=y^2.svg - Wikimedia Commons

Any advice? 

Thank you!

I think I am doing this question right. Is there anything to improve the code? 

How can I get exactly graph plotting as the sample show? All blue area under the curve. 

 


 

" ref: https://www*maplesoft*com/applications/view*aspx"?SID=153480&view=html"  Data :=  [ [0,0.0,0.0],      [5,0.2,0.1],      [10,0.25,0.2],      [18,0.3,0.3],      [25,0.4,0.4],      [32,0.6,0.45],      [38,0.72,0.50],      [43,0.6,0.40],      [47,0.3,0.30],      [50,0.0,0.0] ]: xdata := [0,5,10,18,25,32,38,43,47,50];  ddata := [0.0,0.2,0.25,0.3,0.4,0.6,0.72,0.6,0.3,0.0];  vdata := [0.0,0.1,0.2,0.3,0.4,0.45,0.5,0.4,0.3,0.0];  "

[0, 5, 10, 18, 25, 32, 38, 43, 47, 50]

 

[0., .2, .25, .3, .4, .6, .72, .6, .3, 0.]

 

[0., .1, .2, .3, .4, .45, .5, .4, .3, 0.]

(1)

with(plots); with(plottools); n := nops(xdata)-2; riverbed := plot([seq([xdata[i], -ddata[i]], i = 1 .. n+2)], style = line, color = blue); readings := seq(plot([[xdata[i], 0], [xdata[i], -ddata[i]]], style = line), i = 1 .. n+2); for i to n do x := xdata[i+1]; d := ddata[i+1]; v := vdata[i+1]; s := sprintf("%4.2f m/s", v); xplot[i] := textplot([x, 0, cat("x", i)], align = above); vplot[i] := textplot([x, -d-0.5e-1, s], font = [HELVETICA, 14]) end do; xtickmarks := seq(xplot[i], i = 1 .. n); velocities := seq(vplot[i], i = 1 .. n); display([readings, riverbed, velocities, xtickmarks], view = [0 .. 51, -1 .. 0], labels = ["river cross section x (m)", "depth (m)"]); n := nops(xdata); for i from 2 to n-1 do x := xdata[i]; d := ddata[i]; xL := xdata[i-1]; xR := xdata[i+1]; poly[i] := polygon([[(xL+x)*(1/2), 0], [(xL+x)*(1/2), -d], [(xR+x)*(1/2), -d], [(xR+x)*(1/2), 0]], color = cyan) end do; polys := seq(poly[i], i = 2 .. n-1); plots[display]([readings, riverbed, polys, velocities, xtickmarks], view = [0 .. 51, -1 .. 0], labels = ["river cross section x (m)", "depth (m)"])

``


 

Download Untitled.mw

When Maple converts sin(x)^n to Latex, the result remain  sin(x)^n.  But in Mathematical typesetting, this is normally written as sin^n(x).   Ofcourse this is only for Latex. In Maple code this not valid.

Is it possible to change Maple's Latex to make it do this automatically? Mathematica does this automatically. Here is an example

restart;
expr:=sin(x)^3+cos(3*x)^5;              
Physics:-Latex(expr)
 
            \sin \left(x \right)^{3}+\cos \left(3 x \right)^{5}

Which when compiled gives

Compare to Latex generated by Mathematica

Which compiles to 

Which is more standard in books and papers, than Maple's version.

Both Maple's Physics:-Latex and latex() command do the same thing.

Is there a way to make it generate the improved version for latex?

Maple 2020.1

 

 

sometimes I get intermediate expressions generated from other operations that contain terms such as exp(x)^n in them. As an example, exp(x)^3.  In Mathematica, it automatically replaces these by exp(3*x). But in Maple I need to force this change.

For purposes of Latex only, I like to change these terms to exp(3*x) before converting the whole expression to Latex, as it is looks much better that way.

expr:=exp(x)^3;
Physics:-Latex(expr)

                \left({\rm e}^{x}\right)^{3}

expr:=exp(3*x);
Physics:-Latex(expr)

                {\rm e}^{3 x}


I found that doing simplify(expr,exp)  does the trick. It changes exp(x)^n to exp(n*x). But I am worried about applying this whole simplification command to the whole expression, which can be very large, and do not want to change it all yet.

I just want to change any occurance of exp() there, and nothing more.

I tried using subsindent to do that, but it does not work on terms in denominator

restart;
expr:=exp(x)^3*sin(x)+3/(exp(x)^n);
subsindets(expr,'exp(anything)^anything',f->simplify(f,exp))

I tried

subsindets(expr,'1/exp(anything)^anything',f->simplify(f,exp))

and it did not work.

I am still not good at subsindent. How to make it change all exp(x)^n to exp(n*x) everywhere?

Hello

I need to detect if different expressions contain a radical. For example

aaa := X3*(alpha[2, 8]*(sqrt(X2/alpha[1, 7])*alpha[3, 8] + X1*alpha[3, 6])*X2 + X3*(alpha[2, 8] + alpha[3, 9]/2))/(X2*alpha[2, 8])

As can be noticed there is a square root in the expression.   

I have tried type(expr,sqrt) and has but to no avail. (I am not sure if I use them as they should though).   

Many thanks

 

Ed

 

Hi:

Every once-in-a-while Maple crashes due to a bug in numerical integration that I reported ages ago, but was never fixed (or if it was, I was never told). Anyway, this time, in addition to crashing, it wiped out my worksheet, and froze the machine. When I recovered, Maple had lost all memory of my preferences, recent files, etc - it was as if I had a clean installation.

I restored the link to the usual style file I use, but it still refuses to give me back the style I want - that is - the classical worksheet. Every command I enter now requires F5 or I will get 2D input which I hate. Does anyone know where the preferences, probably an .INI file, are kept, so I can restore it from an older version of Maple.

Thank you

Hi, i am using solve and solve command to find the root but when i used fsolve command to separate only real root, could separate all roots, can anyone correct me, please


 

restart

f := 9.765625000*10^(-6)*(-6671.221362*(x^2+2)^5*sqrt(2)*arctan((1/2)*x*sqrt(2))*x-555.9351135*(x^2+2)^6/((1/2)*x^2+1)-10479.13001*(x^2+2)^5*sqrt(2)*x-(374220*(0.297116730e-1*x^9+.269385824*x^7+.99643086*x^5+5.18951288*x^3+4.42867382*x))*x-1111.870227*x^10-12601.19538*x^8-62147.39274*x^6-485504.8775*x^4-828649.1585*x^2-788850.2769)/(x^2+2)^6-(0.1171875000e-3*(-555.9351135*(x^2+2)^6*sqrt(2)*arctan((1/2)*x*sqrt(2))-873.2608343*(x^2+2)^6*sqrt(2)-(374220*(0.29711673e-2*x^10+0.33673228e-1*x^8+.16607181*x^6+1.29737822*x^4+2.21433691*x^2+2.107985348))*x))*x/(x^2+2)^7+(3.484800000*sqrt(2)*(x^2+2)*arctan((1/2)*x*sqrt(2))*x+.8712000000*(x^2+2)^2/((1/2)*x^2+1)+(5.473911040*(x^2+2))*sqrt(2)*x+5.227200000*x^2-22.99200001)/(16*(x^2+2)^2)-(.8712000000*sqrt(2)*(x^2+2)^2*arctan((1/2)*x*sqrt(2))+1.368477760*sqrt(2)*(x^2+2)^2-36*x*(-0.484000000e-1*x^2+.638666667))*x/(4*(x^2+2)^3)

0.9765625000e-5*(-6671.221362*(x^2+2)^5*2^(1/2)*arctan((1/2)*x*2^(1/2))*x-555.9351135*(x^2+2)^6/((1/2)*x^2+1)-10479.13001*(x^2+2)^5*2^(1/2)*x-374220*(0.297116730e-1*x^9+.269385824*x^7+.99643086*x^5+5.18951288*x^3+4.42867382*x)*x-1111.870227*x^10-12601.19538*x^8-62147.39274*x^6-485504.8775*x^4-828649.1585*x^2-788850.2769)/(x^2+2)^6-0.1171875000e-3*(-555.9351135*(x^2+2)^6*2^(1/2)*arctan((1/2)*x*2^(1/2))-873.2608343*(x^2+2)^6*2^(1/2)-374220*(0.29711673e-2*x^10+0.33673228e-1*x^8+.16607181*x^6+1.29737822*x^4+2.21433691*x^2+2.107985348)*x)*x/(x^2+2)^7+(1/16)*(3.484800000*2^(1/2)*(x^2+2)*arctan((1/2)*x*2^(1/2))*x+.8712000000*(x^2+2)^2/((1/2)*x^2+1)+5.473911040*(x^2+2)*2^(1/2)*x+5.227200000*x^2-22.99200001)/(x^2+2)^2-(1/4)*(.8712000000*2^(1/2)*(x^2+2)^2*arctan((1/2)*x*2^(1/2))+1.368477760*2^(1/2)*(x^2+2)^2-36*x*(-0.484000000e-1*x^2+.638666667))*x/(x^2+2)^3

(1)

ip := solve(f = 0, x)

.6540411301, 3126.002498+5414.398621*I, .4137989369+1.038962897*I, .6364817315+1.870977651*I, -.6364817315+1.870977651*I, -.4137989369+1.038962897*I, -.6540411301, -6252.010299, -.4137989369-1.038962897*I, -.6364817315-1.870977651*I, .6364817315-1.870977651*I, .4137989369-1.038962897*I, 3126.002498-5414.398621*I

(2)

cp := fsolve(numer(f) = 0, x)

.6540411302

(3)

``


 

Download help_fsolve_real_root.mw

Hi there.

As we all know if we multiply two polynomials f(x) and g(x) of degrees m and n respectively we get polynomial h(x)= f(x)*g(x) of degree m+n and with m+n+1 coefficients in general. Function modp1(('Multiply')(...)) doing this very well. But sometimes we don't need full resulting h(x) - just subset of monomials and subset of coefficients of h(x) - so we don't need to calculate all m+n+1 coefficients of h(x) and waste time and resources for that.

I would request some additional rework of modp1 package: by adding to modp1(('Multiply')(...)) two optional parameters - degrees of first and last calculating coefficients of h(x).

For example:

h:=modp1(Multiply(f, g,n-1,n+1), p) could calculate only monomials with n-1, n and n+1 degrees and set other monomials to zero.

Or maybe it should be new function:

h:=modp1(Multiply_Truncate(f, g,n-1,n+1), p)

 

Is it possible?

It would be great and very efficient in many tasks.

Thank you.


 

restart;

M__h := 0.352e-1;

0.352e-1

 

0.34e-1

 

0.8354e-1

 

0.96e-2

 

.123

 

0.7258e-1

 

0.214e-1

 

0.219e-1

 

.123

 

.7902

 

.11

 

0.136e-3

 

0.5e-1

 

0.8910e-1

 

0.45e-1

 

.7

 

.7214

 

1.354

 

0.235e-1

(1)

pdes := [diff(B(t, x), t) = M__h-beta__1*B(t, x)*G(t, x)/N__h+beta__2*B(t, x)*G(t, x)/N__h-mu__h*B(t, x)+sigma__h*E(t, x)*(diff(B(t, x), x, x)), diff(C(t, x), t) = beta__1*B(t, x)*G(t, x)/N__h-u[1]*C(t, x)/(1+C(t, x))-mu__h*C(t, x)*(diff(C(t, x), x, x)), diff(DD(t, x), t) = beta__2*DD(t, x)*G(t, x)/N__h-u[1]*DD(t, x)/(1+DD(t, x))-mu__h*DD(t, x)-delta__1*DD(t, x)*(diff(DD(t, x), x, x)), diff(E(t, x), t) = u[1]*C(t, x)/(1+C(t, x))+u[1]*DD(t, x)/(1+DD(t, x))-(mu__h+sigma__h)*E(t, x)*(diff(E(t, x), x, x)), diff(F(t, x), t) = M__b-beta__3*F(t, x)*C(t, x)/N__b+beta__4*F(t, x)*DD(t, x)/N__b-mu__b*F(t, x)*(diff(F(t, x), x, x)), diff(G(t, x), t) = beta__3*F(t, x)*C(t, x)/N__b+beta__4*F(t, x)*DD(t, x)/N__b-mu__b*G(t, x)*(diff(G(t, x), x, x))];

[diff(B(t, x), t) = 0.352e-1-0.891056911e-1*B(t, x)*G(t, x)-0.96e-2*B(t, x)+0.8910e-1*E(t, x)*(diff(diff(B(t, x), x), x)), diff(C(t, x), t) = .6791869919*B(t, x)*G(t, x)-0.45e-1*C(t, x)/(1+C(t, x))-0.96e-2*C(t, x)*(diff(diff(C(t, x), x), x)), diff(DD(t, x), t) = .5900813008*DD(t, x)*G(t, x)-0.45e-1*DD(t, x)/(1+DD(t, x))-0.96e-2*DD(t, x)-0.235e-1*DD(t, x)*(diff(diff(DD(t, x), x), x)), diff(E(t, x), t) = 0.45e-1*C(t, x)/(1+C(t, x))+0.45e-1*DD(t, x)/(1+DD(t, x))-0.9870e-1*E(t, x)*(diff(diff(E(t, x), x), x)), diff(F(t, x), t) = .7214-.1739837398*F(t, x)*C(t, x)+.1780487805*F(t, x)*DD(t, x)-1.354*F(t, x)*(diff(diff(F(t, x), x), x)), diff(G(t, x), t) = .1739837398*F(t, x)*C(t, x)+.1780487805*F(t, x)*DD(t, x)-1.354*G(t, x)*(diff(diff(G(t, x), x), x))]

(2)

bcs := [(D[2](B))(t, 0) = 0, (D[2](B))(t, 1) = 0, (D[2](C))(t, 0) = 0, (D[2](C))(t, 1) = 0, (D[2](DD))(t, 0) = 0, (D[2](DD))(t, 1) = 0, (D[2](E))(t, 0) = 0, (D[2](E))(t, 1) = 0, (D[2](F))(t, 0) = 0, (D[2](F))(t, 1) = 0, (D[2](G))(t, 0) = 0, (D[2](G))(t, 1) = 0, B(0, x) = 100, C(0, x) = 70, DD(0, x) = 50, E(0, x) = 70, F(0, x) = 100, G(0, x) = 70]

[(D[2](B))(t, 0) = 0, (D[2](B))(t, 1) = 0, (D[2](C))(t, 0) = 0, (D[2](C))(t, 1) = 0, (D[2](DD))(t, 0) = 0, (D[2](DD))(t, 1) = 0, (D[2](E))(t, 0) = 0, (D[2](E))(t, 1) = 0, (D[2](F))(t, 0) = 0, (D[2](F))(t, 1) = 0, (D[2](G))(t, 0) = 0, (D[2](G))(t, 1) = 0, B(0, x) = .100, C(0, x) = .70, DD(0, x) = .50, E(0, x) = .70, F(0, x) = .100, G(0, x) = .70]

(3)

sol := pdsolve(pdes, bcs, numeric);

module () local INFO; export plot, plot3d, animate, value, settings; option `Copyright (c) 2001 by Waterloo Maple Inc. All rights reserved.`; end module

(4)

sol:-plot3d([B(t, x), C(t, x)], t = 0 .. 20, x = 0 .. 20)

Error, (in pdsolve/numeric/plot3d) unable to compute solution for t>HFloat(0.25):
Newton iteration is not converging

 

``


 

Download spatial_1.mw

Not sure part A is done correctly. flow (m^3 / s) = Area (M^2) * velocity (m/s)

And have no idea how to do part B in maple

Thank you!

How to solve the system equation with Matrix? Basically is part (b)

And not quite familiar with Matrix multiplication in maple? How to do part(c)?

Thank you!!!!

Hi all,

We want to find a curve fit for an integer sequence.

We have n such that n^2+n+17 is a prime number.

See oeis.org/A007635 and comments.

Use the Maple CurveFitting package.

I tried with(CurveFitting).

We do not know if this is best represented by a polynomial or exponential curve fit.

 

n2_and_n_and_17_in_OEIS_007635.mw

n2_and_n_and_17_in_OEIS_007635.pdf

Regards,

Matt

The low-level command is "RTABLE". I have used it a million times. I use it with Maple 13 to replace expressions by tables.
For example,

print(`Hello! `*RTABLE(1,x^2+y^2));

will print: Hello! x^2+y^2. However, the main reason I am using it is to avoid Maple messing up the order of things like sums,
products, etc. However, this command has disappeared. It does not exist in Maple 2020. So, I cannot run my programs at the university which now has Maple 2020. Anybody knows if this command can still be used in another form the same way?

Thank you!
mapleatha

First 486 487 488 489 490 491 492 Last Page 488 of 2427