MaplePrimes Questions

Hello everyone,
My name's Rafael, I am a master's student in Dentistry in Brazil.
I need your help.
I have several points with X, Y and Z coordinates. I would like to create a curve that represents these points in a 3-dimensional plane.
I'm new using Maple and managed to use the "plots [pointplot3d]" function to represent the points in the chart. However I have a list in excel with several point (over 100) and would like to know how to import these points to use this function.
Thank you very much


How can you get maple to evaluate i^i?

when i type in
I^I

i just get

I^I

and similarly when i raise numbers to complex powers i get results like 2^(2I+6)

 

Recently, my questions are deleted. last month I ask 3 questions but that all questions are deleted without my knowledge. yesterday also i ask one question "how to plot3d graph" that also deleted. "https://www.mapleprimes.com/questions/225618-How-To-Plot-The-3d-Graph?sq=225618". In the mapleprime, I am following the ethics and rules. 

Hi,

This more a warning to focus your attention on a specific point than a true question.
 

I submit you this test case which works in Maple 2015 and Maple 2016 but not in Maple 2018.

In a few words:

  • let X and Y two independant random variables with respective distributions Normal(mu__x, sigma__x) and  
    Normal(mu__y, sigma__y)
     
  • let Z := q -> cos(q)
     
  • You can easily verify that Maple can compute the formal expression of Mean(Z(X)) and Variance(Z(X))
    (which means that it could compute Z(X+Y) for X+Y is just another gaussian RV)
     
  • What I found is that:
    1. Mean(Z(X+Y)) returns same expressions in Maple 2015 and Maple 2016, but a different one in Maple 2018.
      Luckily the later is more readable than the former ones, and closer to the one of Mean(Z(U))  where U=X+Y is the RV of distribution Normal(mu__x+mu__y, sqrt(sigma__x^2+sigma__y^2))
      This suggest that the integration algorithm has evolved somewhere in between Maple 2016 and Maple 2018
       
    2. While Maple 2015 and Maple 2016 return an evaluated result for Variance(Z(X+Y)) Maple 2018 fails.
       

Can this "failure" be fixed by some adhoc option of Variance?
Or could it come from a "regression" in the implementation of this procedure (or of the underlying int procedure) in Maple 2018?

PS: I did not try to compute Variance(Z(X+Y)) from an explicit double integration


Stat_2015.mw

Why does the following not get simplified?

simplify(conjugate(b) * I - conjugate(a) + conjugate(b*I+a));

It seems Maple does not know conjugation is linear, since the following

simplify(conjugate(b) * I                + conjugate(b*I  ));
simplify(                 - conjugate(a) + conjugate(    a));

do get simplified to 0.

The function does simplify to 0, but it assumes my and are real, which they are not. 

Why does the following code 

evalb(a * conjugate(a) = abs(a) ^ 2);

not produce true?

In ScientificConstants, the density of carbon is given as 2.2 g/cm^3. That happens to be the value for graphite. Carbon comes in different states with different density, diamond being one of them, amorphous being another. Is this somewhere in ScientificConstants where just I cannot find it, or is it not there?

TIA,

M.D.

PS: I do know the values and where to find them (Particle Data Group). I'd like to use ScientificConstants so I can write a more general piece of code.

My worksheet is shown below. I am copying an example from a textbook and know that h = 26.4 is correct.

However, if I use !!! (execute entire worksheet) then I get h = -13

What an i doing wrong

Thanks

Hi 

I am having problems with the syntax/symbols i maple. I down know what "setting" has changed but + and - and so on, has change to K and C. anyone know how to fix this?

Regards

Morten

 

I have just found that the standard routine of Maple for performing mod 2 multivariate factorization, that is
Factor(p) mod 2
takes a very long computing time (I stopped after a long while) for the polynomial

p := K[4]*(K[4]^2*K[5]*K[6]^3*K[7]^2+K[4]*K[5]^2*K[6]^4*K[7]+K[4]^3*K[6]^2*K[7]^2+K[4]^2*K[5]*K[6]^3*K[7]+K[4]^2*K[5]*K[6]^2*K[7]^2+K[4]*K[5]^3*K[6]^3+K[5]^3*K[6]^4+K[4]^2*K[5]^3*K[6]+K[4]^2*K[5]^2*K[6]^2+K[4]^3*K[5]^2+K[4]^2*K[5]^3);

Note that the factorization of this polynomial is very fast on other computer algebra systems. Since mod 2 multivariate factorization (in fact, normalization of mod 2 multivariate rational functions) is essential for some algorithms that I'm currently implementing in Maple, I would ask if you know how to improve (bug, additional parameters?) the function Factor or if there are other functions, packages in Maple which implement a more feasible mod 2 multivariate factorization.

Many thanks in advance.

Define a piecewise function of  so that the the value of the function is  if  or if , it is equal to the vertical distance between the curves if . Numerically integrate the area under this curve on the non-zero interval and explain how the numerical answer is what you would anticipate the answer to be given that this curve is related to semi-circles.

The curve is defined by

Hi everyone,

Was hoping someone could point me in the right direction. A homework question asks us to create a function using arrow (->) notation which take a Maple list of complex numbers as its input and returns the largest modulus from that list.

I've tried writing this without arrow notation --

Say L1 is a list of complex numbers:

max(abs~(L1))

and this seems to work fine. But I'm particularly struggling with how to introduce the arrow notation into this.

Any help would be much appreciated!

Hi. Can anyone help me with this error, it says 

Error, (in unknown) cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up

 

restart

Shootlib := "F:\MAPLE\3";

"F:MAPLE3"

(1)

libname := Shootlib, libname;

"F:MAPLE3", "C:\Program Files\Maple 2015\lib", "."

(2)

with(Shoot);

[shoot]

(3)

with(plots):

``

FNS := {`φp`(eta), f(eta), fp(eta), fpp(eta), fppp(eta), phi(eta), theta(eta), `θp`(eta)}

{`φp`(eta), f(eta), fp(eta), fpp(eta), fppp(eta), phi(eta), theta(eta), `θp`(eta)}

(4)

ODE := {diff(`φp`(eta), eta) = (-f(eta)*`φp`(eta)+K*phi(eta)*(1-phi(eta))^2)*Sc, diff(f(eta), eta) = fp(eta), diff(fp(eta), eta) = fpp(eta), diff(fpp(eta), eta) = fppp(eta), diff(fppp(eta), eta) = (-fppp(eta)-f(eta)*fpp(eta)+k1*(2*fp(eta)*fppp(eta)-fpp(eta)^2)+lambda*fp(eta)+(1+Fr)*fp(eta)^2)/(k1*f(eta)), diff(phi(eta), eta) = `φp`(eta), diff(theta(eta), eta) = `θp`(eta), diff(`θp`(eta), eta) = (-f(eta)*`θp`(eta)+y*f(eta)*fp(eta)*`θp`(eta))/(1/Pr-y*f(eta)^2)}

{diff(`φp`(eta), eta) = (-f(eta)*`φp`(eta)+K*phi(eta)*(1-phi(eta))^2)*Sc, diff(f(eta), eta) = fp(eta), diff(fp(eta), eta) = fpp(eta), diff(fpp(eta), eta) = fppp(eta), diff(fppp(eta), eta) = (-fppp(eta)-f(eta)*fpp(eta)+k1*(2*fp(eta)*fppp(eta)-fpp(eta)^2)+lambda*fp(eta)+(1+Fr)*fp(eta)^2)/(k1*f(eta)), diff(phi(eta), eta) = `φp`(eta), diff(theta(eta), eta) = `θp`(eta), diff(`θp`(eta), eta) = (-f(eta)*`θp`(eta)+y*f(eta)*fp(eta)*`θp`(eta))/(1/Pr-y*f(eta)^2)}

(5)

blt := 4.; 1; IC := {`φp`(0) = Psi, f(0) = 0, fp(0) = 1, fpp(0) = 0, fppp(0) = alpha, phi(0) = 0, theta(0) = 1, `θp`(0) = beta}

4.

 

{`φp`(0) = Psi, f(0) = 0, fp(0) = 1, fpp(0) = 0, fppp(0) = alpha, phi(0) = 0, theta(0) = 1, `θp`(0) = beta}

(6)

BC := {fp(blt) = 0, fpp(blt) = 0, phi(blt) = 1, theta(blt) = 0};

{fp(4.) = 0, fpp(4.) = 0, phi(4.) = 1, theta(4.) = 0}

(7)

``

infolevel[Shoot:-shoot] := 1:

.7

 

1.0

 

.2

 

.2

 

1.0

 

.3

 

.5

 

.3

(8)

S := shoot(ODE, IC, BC, FNS, [alpha = 3.0125, beta = -2.31, Psi = -0.])

shoot: Step #  1

shoot: Parameter values :  alpha = 3.0125 beta = -2.31 Psi = -0.

Error, (in unknown) cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up

 

p := odeplot(S, [eta, fp(eta)], 0 .. 15)

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

 

display(p)

Error, (in plots:-display) expecting plot structure but received: p

 

p2 := odeplot(S, [eta, theta(eta)], 0 .. 10)

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

 

display(p2)

Error, (in plots:-display) expecting plot structure but received: p2

 

p3 := odeplot(S, [eta, phi(eta)], 0 .. 10)

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution

 

display(p3)

Error, (in plots:-display) expecting plot structure but received: p3

 


 

Download Darcy.mw

Below are five subsindets commands.

I believe I understand the actions of B and C, but I fail to understand the actions, individually and taken together, of  E, F and G.

 

B := subsindets(u(i, j)^2*v(i, j)+u(i-1, j), 'specfunc(u)', proc (f) options operator, arrow; subsop(1 = op(1, f)+1, f) end proc);

C := subsindets(u(i, j)^2*v(i, j)+u(i-1, j), 'specfunc(anything, u)', proc (f) options operator, arrow; subsop(1 = op(1, f)+1, f) end proc);

E := subsindets(u(i, j)^2*v(i, j)+u(i-1, j), 'specfunc(symbol, u)', proc (f) options operator, arrow; subsop(1 = op(1, f)+1, f) end proc);

F := subsindets(u(i, j)^2*v(i, j)+u(i-1, j), 'specfunc(`+`, u)', proc (f) options operator, arrow; subsop(1 = op(1, f)+1, f) end proc);

G := subsindets(u(i, j)^2*v(i, j)+u(i-1, j), 'specfunc({`+`, symbol}, u)', proc (f) options operator, arrow; subsop(1 = op(1, f)+1, f) end proc);

Where can I find a thorough explanation of specfunc with examples?

First 768 769 770 771 772 773 774 Last Page 770 of 2427