Maple 2017 Questions and Posts

These are Posts and Questions associated with the product, Maple 2017

HOW I can find the function that satisfies these boundary conditions?

please help me.

Thanks

FUNCTION.mw

Download FUNCTION.mw

 

equation of ellipse

eqn := (1/32)*(x-16)^2+(1/2025)*(y+0)^2 = 1

only interested in the positive section

plots[implicitplot]([(1/32)*(x-16)^2+(1/2025)*(y+0)^2 = 1], x = 0 .. 18, y = 0 .. 60, scaling = constrained)

equation for arch length of function

int(sqrt(1+(diff(eqn, x))^2), x = 0 .. 16)

Doesn't seem like the best way to solve this.

Wanting to generate a general fomula I can then use in excel to calculate the arch length of the positive only section of the ellipse

any ideas of what the assumptions could be to have maple solve this?

 

eqn2 := (x-h)^2/a^2+(y-k)^2/b^2 = 1;
                               2          2    
                        (x - h)    (y - k)     
                eqn2 := -------- + -------- = 1
                            2          2       
                           a          b        
                             "(->)"

     [[                                            (1/2)]  
     [[          / 2  2    2  2      2        2  2\     ]  
     [[    a k + \a  b  - b  h  + 2 b  h x - b  x /     ]  
     [[y = ---------------------------------------------], 
     [[                          a                      ]  

       [                                            (1/2)]]
       [          / 2  2    2  2      2        2  2\     ]]
       [    a k - \a  b  - b  h  + 2 b  h x - b  x /     ]]
       [y = ---------------------------------------------]]
       [                          a                      ]]

eqn3 := (a*k+sqrt(a^2*b^2-b^2*h^2+2*b^2*h*x-b^2*x^2))/a;
                                                     (1/2)
                   / 2  2    2  2      2        2  2\     
             a k + \a  b  - b  h  + 2 b  h x - b  x /     
     eqn3 := ---------------------------------------------
                                   a                      
diff(eqn3, x);
                           2        2                
                        2 b  h - 2 b  x              
          -------------------------------------------
                                              (1/2)  
            / 2  2    2  2      2        2  2\       
          2 \a  b  - b  h  + 2 b  h x - b  x /      a
`assuming`([int(sqrt(1+(diff(eqn3, x))^2), x = c .. d)], [a > 0, b > 0, c >= 0, d > c]);
 

Hello everyone,

I want to solve the matrix equation UTCU=M, where I know the matrices C and M and I need to find que matrix U.

I already did so many tests but without success.

Could someone help me?

Thanks!

I am unable to do the following multiplication of two vectors(B and C given below), which is possible in Mathematics but not in Maple(Tried all the possible ways)

 

B := Vector(2, {1 = 2, 2 = 3});
C := Vector(1, {1 = 4});


B . C;
Error, (in LinearAlgebra:-DotProduct) vectors must have the same dimension
B*C;
Error, (in rtable/Product) use *~ for elementwise multiplication of Vectors or Matrices; use . (dot) for Vector/Matrix multiplication
`~`[`*`](B, C);
Error, dimension bounds must be the same for all container objects in an elementwise operation
 

I have the following two questions/issues:

1.) It seems to me that in the Tetrads package, some tetrad is automatically calculated in terms of the metric. Is there any way to reverse the roles, i.e., to have the vierbein determine the metric?, as is anyway a more natural setup, I think. The point is that I would like to work with an arbitrary, non-fixed vierbein, and from it having the metric, Christoffel symbols, etc., calculated 'on the fly', and then ultimatively construct the socalled minimal spin connection from these quantities.

2.) Consider the following expressions:

Simplify(e_[a,mu ]*Christoffel[~mu,nu,rho]);
Simplify(e_[a,nu ]*Christoffel[mu,~nu,rho]);
Simplify(e_[a,rho]*Christoffel[mu,nu,~rho]);

As the Christoffel symbols do not transform as tensors, I find this notation unfortunate, as one cannot simply move the vierbein through a (partial) derivative. Conversion between world indices and Lorentz indices should, I think, only be performed for world indices of tensorial type.

Is there something wrong with dsolve?

ode := diff(y(x), x) = sqrt(2*32.2*y(x)):
ics := y(0) = 0:
dsolve({ics, ode});

maple output:   y(x) = 0

The answer should be

y(x) = 16.1 x^2

Wolfram got it

 

I tried using restart, with(DEtools), still no luck. Though I don't think its necessary to call with(DEtools) for this simple equation.

To make the problem simpler, use dsolve on  dy/dx = √y , y(0)=0.

ode := diff(y(x), x) = sqrt(y(x)):
ics := y(0) = 0:
dsolve({ics, ode});

maple returns  y(x) = 0, which is incorrect.  Should be y(x) = x^2/4

I uploaded the worksheet just in case 'its just me'.

diffeqseperable.mw

 

Download test.mw
Udp: document attached.

Q__gr := Typesetting[delayDotProduct](Vector[row](8, {(1) = 0, (2) = 5, (3) = 10, (4) = 15, (5) = 20, (6) = 25, (7) = 30, (8) = 35}), Unit('m'^3/'day'), true)

Vector[row](%id = 18446746861178193550)

(1)

`η__gr` := Vector[row](8, {(1) = 0, (2) = 9.324, (3) = 17, (4) = 23.232, (5) = 27.192, (6) = 29.6, (7) = 29.516, (8) = 24.592})

Vector[row](%id = 18446746861214593262)

(2)

`pointsη` := [seq([Q__gr[i], `η__gr`[i]], i = 1 .. 8)]

[[0, 0], [5*Units:-Unit(m^3/d), 9.324], [10*Units:-Unit(m^3/d), 17], [15*Units:-Unit(m^3/d), 23.232], [20*Units:-Unit(m^3/d), 27.192], [25*Units:-Unit(m^3/d), 29.6], [30*Units:-Unit(m^3/d), 29.516], [35*Units:-Unit(m^3/d), 24.592]]

(3)

CurveFitting[PolynomialInterpolation](`pointsη`, q)

0.2348698413e-7*q^7/Units:-Unit(m^3/d)^7-0.2980622223e-5*q^6/Units:-Unit(m^3/d)^6+0.1482222222e-3*q^5/Units:-Unit(m^3/d)^5-0.3663955556e-2*q^4/Units:-Unit(m^3/d)^4+0.4666528890e-1*q^3/Units:-Unit(m^3/d)^3-.3165382223*q^2/Units:-Unit(m^3/d)^2+2.655161905*q/Units:-Unit(m^3/d)

(4)

eta := proc (q) options operator, arrow; 0.2348698413e-7*q^7/Units:-Unit(m^3/d)^7-0.2980622223e-5*q^6/Units:-Unit(m^3/d)^6+0.1482222222e-3*q^5/Units:-Unit(m^3/d)^5-0.3663955556e-2*q^4/Units:-Unit(m^3/d)^4+0.4666528890e-1*q^3/Units:-Unit(m^3/d)^3-.3165382223*q^2/Units:-Unit(m^3/d)^2+2.655161905*q/Units:-Unit(m^3/d) end proc

proc (q) options operator, arrow; 0.2348698413e-7*q^7/Units:-Unit(m^3/d)^7-0.2980622223e-5*q^6/Units:-Unit(m^3/d)^6+0.1482222222e-3*q^5/Units:-Unit(m^3/d)^5-0.3663955556e-2*q^4/Units:-Unit(m^3/d)^4+0.4666528890e-1*q^3/Units:-Unit(m^3/d)^3-.3165382223*q^2/Units:-Unit(m^3/d)^2+2.655161905*q/Units:-Unit(m^3/d) end proc

(5)

eta(12*Unit('m'^3/'day'))

.8415811058*Units:-Unit(m^3/d)^7/Units:-Unit(m^3/d)^7-8.900090268*Units:-Unit(m^3/d)^6/Units:-Unit(m^3/d)^6+36.88243199*Units:-Unit(m^3/d)^5/Units:-Unit(m^3/d)^5-75.97578241*Units:-Unit(m^3/d)^4/Units:-Unit(m^3/d)^4+80.63761922*Units:-Unit(m^3/d)^3/Units:-Unit(m^3/d)^3-45.58150401*Units:-Unit(m^3/d)^2/Units:-Unit(m^3/d)^2+31.86194286*Units:-Unit(m^3/d)/Units:-Unit(m^3/d)

(6)

simplify(19.76619849)

Error, (in Units:-Standard:-+) the units `m^21/s^7` and `m^18/s^6` have incompatible dimensions

 

``


 

Download test.mw

 

Hi!

I got a function (from CurveFitting), that produce a polynomial with the some units inside:

 

As i can see, it can be easy simplified, but i get:


 

Is where any trick to do it, without stripping units with convert(unit_free) or something?

Thank you!

 

  1. (k-2)*(k^2+5)*(k^3-k^2+7*k+8)/(6*k*(k^2-3*k+8))
  2. (k+2)*(k^2+5)*(k^3-5*k^2+13*k-8)/(6*k*(k^2-3*k+8))
  3. k^3+3*k^2+11*k-3
  4. k^2+2*k+9

Hello everybody

I am currently having trouble 'manually simplifying' an equation that I differentiated using maple. Basically, the equation that I get involves hessian-matrices. Since the resulting equation is rather long, I would like to replace every such matrix by a sign, say H.
This is easy to do for non-matrix equations, i.e.

subs(sin(x)sin(y) = z, sin(x)sin(y) + xy)

will give me z + xy, which looks a bit easier.

However, I couldn't figure out how to do a similar thing with matrices. For example, the following code

subs(Matrix(3, 3, [[x, y, z], [y, z, x], [z, x, y]]) = A, Matrix(3, 3, [[x, y, z], [y, z, x], [z, x, y]])+B)

 will not work.
I already thought about converting the matrices first into lists (because for some reason it works for lists). However, I would also like to do the same as above for for functions of matrices, i.e. set F(A) = G, where A is a matrix (this also fails, probably for the same reason).

Also, the same thing seems to fail for vectors and general arrays, so I guess the actual problem might be the array type.

I also tried alternative ways, such as eval or algsubs, but the later even gives an error since it cannot handle matrices at all.

This feels a bit like a noob question, but I spent almost 2 days now searching for an answer or a workaround, so my apologies if I missed something trivial.

All the best
Adrian

 

 

Hi, I try to solve the below integral. when I press enter key maple dosen't show answer and show the integral again.

int(r*r[bc]*r[tc], r = r[bc] .. r[tc]);

but when I write intgral this way and use " i " as subscript ,maple solve it.

int(r[i]*r[bc]*r[tc], r[i] = r[bc] .. r[tc]);

I just want to know why?

what is difference between first and second integral?

and also is there any way (or any packages) to solve these integrals?

(I read https://www.maplesoft.com/applications/view.aspx?SID=6846&view=html article befor)

thanks.

  1. (2*k^3-6*k^2+7*k+15-k*sqrt(k^6-12*k^5+64*k^4-198*k^3+448*k^2-636*k+369))/(-k^4+2*k^3-2*k^2+10*k+15)
  2. (k^3+5*k+(-k^2+k)*sqrt(k^4-10*k^3+37*k^2-60*k+180)+30)/(-k^4+k^3+k^2+5*k+30)

hi

please help me for simplify (factor) this equations.

thanks

vel.mw
 

simplify(-(1/226609908940800)*(106722*Br*NT*ln(h)*NB-106722*Br*NT*ln(R0)*NB-106722*Br*NT^2*ln(h)+106722*Br*NT^2*ln(R0)-106722*NB^2*Gr*ln(h)+106722*NB^2*Gr*ln(R0)+106722*Gr*NT*ln(h)*NB-106722*Gr*NT*ln(R0)*NB)*r^12/(L^6*NB*(ln(h)-ln(R0)))-(1/226609908940800)*(-524288*Br*NT*ln(h)*R0*NB-524288*Br*NT*ln(h)*h*NB+524288*Br*NT*ln(R0)*R0*NB+524288*Br*NT*ln(R0)*h*NB+524288*Br*NT^2*ln(h)*R0+524288*Br*NT^2*ln(h)*h-524288*Br*NT^2*ln(R0)*R0-524288*Br*NT^2*ln(R0)*h+524288*NB^2*Gr*ln(h)*R0+524288*NB^2*Gr*ln(h)*h-524288*NB^2*Gr*ln(R0)*R0-524288*NB^2*Gr*ln(R0)*h-524288*Gr*NT*ln(h)*R0*NB-524288*Gr*NT*ln(h)*h*NB+524288*Gr*NT*ln(R0)*R0*NB+524288*Gr*NT*ln(R0)*h*NB)*r^11/(L^6*NB*(ln(h)-ln(R0)))-(1/226609908940800)*(1920996*Br*NT*ln(h)*R0*h*NB-15367968*L^2*Br*NT^2*ln(h)+15367968*L^2*Br*NT^2*ln(R0)-15367968*L^2*Gr*NB^2*ln(h)+15367968*L^2*Gr*NB^2*ln(R0)-960498*NB^2*Gr*ln(h)*R0^2-960498*NB^2*Gr*ln(h)*h^2+960498*NB^2*Gr*ln(R0)*R0^2-1920996*Br*NT*ln(R0)*R0*h*NB+1920996*Gr*NT*ln(h)*R0*h*NB-1920996*Gr*NT*ln(R0)*R0*h*NB+960498*NB^2*Gr*ln(R0)*h^2-960498*Br*NT^2*ln(h)*h^2+960498*Br*NT^2*ln(R0)*R0^2-960498*Br*NT^2*ln(h)*R0^2+960498*Br*NT^2*ln(R0)*h^2+1920996*Br*NT^2*ln(R0)*R0*h-960498*Gr*NT*ln(R0)*h^2*NB+960498*Gr*NT*ln(h)*h^2*NB-960498*Gr*NT*ln(R0)*R0^2*NB+960498*Gr*NT*ln(h)*R0^2*NB-1920996*NB^2*Gr*ln(h)*R0*h+1920996*NB^2*Gr*ln(R0)*R0*h-15367968*L^2*Gr*NB*NT*ln(R0)-960498*Br*NT*ln(R0)*h^2*NB+15367968*L^2*Gr*NB*NT*ln(h)+960498*Br*NT*ln(h)*h^2*NB-960498*Br*NT*ln(R0)*R0^2*NB-15367968*L^2*Br*NB*NT*ln(R0)+960498*Br*NT*ln(h)*R0^2*NB-1920996*Br*NT^2*ln(h)*R0*h+15367968*L^2*Br*NB*NT*ln(h))*r^10/(L^6*NB*(ln(h)-ln(R0)))-(1/226609908940800)*(-63438848*L^2*Br*NB*NT*ln(h)*R0-63438848*L^2*Br*NB*NT*ln(h)*h+63438848*L^2*Br*NB*NT*ln(R0)*R0+63438848*L^2*Br*NB*NT*ln(R0)*h+63438848*L^2*Br*NT^2*ln(h)*R0+63438848*L^2*Br*NT^2*ln(h)*h-63438848*L^2*Br*NT^2*ln(R0)*R0-63438848*L^2*Br*NT^2*ln(R0)*h+63438848*L^2*Gr*NB^2*ln(h)*R0+63438848*L^2*Gr*NB^2*ln(h)*h-63438848*L^2*Gr*NB^2*ln(R0)*R0-63438848*L^2*Gr*NB^2*ln(R0)*h-63438848*L^2*Gr*NB*NT*ln(h)*R0-63438848*L^2*Gr*NB*NT*ln(h)*h+63438848*L^2*Gr*NB*NT*ln(R0)*R0+63438848*L^2*Gr*NB*NT*ln(R0)*h)*r^9/(L^6*NB*(ln(h)-ln(R0)))-(1/226609908940800)*(11116875*Br*NT^2*R0^4-11116875*Br*NT^2*h^4-11116875*NB^2*Gr*h^4+1536796800*p*ln(h)*NB+11116875*NB^2*Gr*R0^4-1536796800*p*ln(R0)*NB+11116875*Gr*NT*h^4*NB-11116875*Br*NT*R0^4*NB-11116875*Gr*NT*R0^4*NB-192099600*L^2*Gr*NB*NT*ln(R0)*R0*h+192099600*L^2*Gr*NB*NT*ln(h)*R0*h-192099600*L^2*Br*NB*NT*ln(R0)*R0*h+192099600*L^2*Br*NB*NT*ln(h)*R0*h-5336100*Br*NT^2*ln(R0)*R0*h^3-11116875*Gr*NT*R0^3*h*NB-11116875*Br*NT*R0^3*h*NB+11116875*Br*NT*R0*h^3*NB+11116875*Gr*NT*R0*h^3*NB+5336100*Gr*NT*ln(R0)*R0*h^3*NB+24012450*Gr*NT*ln(R0)*R0^2*h^2*NB-24012450*Gr*NT*ln(h)*R0^2*h^2*NB-5336100*Gr*NT*ln(h)*R0^3*h*NB-5336100*Gr*NT*R0*h^3*ln(r)*NB+5336100*Gr*NT*R0^3*h*ln(r)*NB-96049800*L^2*Gr*NB*NT*ln(R0)*h^2-96049800*L^2*Gr*NB*NT*ln(R0)*R0^2+96049800*L^2*Gr*NB*NT*ln(h)*h^2+96049800*L^2*Gr*NB*NT*ln(h)*R0^2+192099600*L^2*Gr*NB^2*ln(R0)*R0*h-192099600*L^2*Gr*NB^2*ln(h)*R0*h+5336100*Br*NT*ln(R0)*R0*h^3*NB+24012450*Br*NT*ln(R0)*R0^2*h^2*NB-24012450*Br*NT*ln(h)*R0^2*h^2*NB-5336100*Br*NT*ln(h)*R0^3*h*NB-5336100*Br*NT*R0*h^3*ln(r)*NB+5336100*Gr*NT*ln(R0)*h^4*NB-5336100*Gr*NT*ln(h)*R0^4*NB-5336100*Gr*NT*h^4*ln(r)*NB+5336100*Br*NT*R0^3*h*ln(r)*NB+192099600*L^2*Br*NT^2*ln(R0)*R0*h-192099600*L^2*Br*NT^2*ln(h)*R0*h-96049800*L^2*Br*NB*NT*ln(R0)*h^2-96049800*L^2*Br*NB*NT*ln(R0)*R0^2+96049800*L^2*Br*NB*NT*ln(h)*h^2+96049800*L^2*Br*NB*NT*ln(h)*R0^2+5336100*Gr*NT*R0^4*ln(r)*NB-5336100*NB^2*Gr*ln(R0)*R0*h^3-24012450*NB^2*Gr*ln(R0)*R0^2*h^2+24012450*NB^2*Gr*ln(h)*R0^2*h^2+5336100*NB^2*Gr*ln(h)*R0^3*h+5336100*NB^2*Gr*R0*h^3*ln(r)-5336100*NB^2*Gr*R0^3*h*ln(r)+96049800*L^2*Gr*NB^2*ln(R0)*h^2+96049800*L^2*Gr*NB^2*ln(R0)*R0^2-96049800*L^2*Gr*NB^2*ln(h)*h^2-96049800*L^2*Gr*NB^2*ln(h)*R0^2-1536796800*L^4*Gr*NB*NT*ln(R0)+5336100*Br*NT^2*R0*h^3*ln(r)+1536796800*L^4*Gr*NB*NT*ln(h)-5336100*Br*NT^2*R0^3*h*ln(r)-5336100*Br*NT*ln(h)*R0^4*NB+5336100*Br*NT*ln(R0)*h^4*NB-5336100*Br*NT*h^4*ln(r)*NB+5336100*NB^2*Gr*ln(h)*R0^4-5336100*NB^2*Gr*ln(R0)*h^4-5336100*NB^2*Gr*R0^4*ln(r)+5336100*NB^2*Gr*h^4*ln(r)-1536796800*L^4*Gr*NB^2*ln(h)+1536796800*L^4*Gr*NB^2*ln(R0)+1536796800*L^4*Br*NT^2*ln(R0)-5336100*Br*NT^2*R0^4*ln(r)+5336100*Br*NT^2*h^4*ln(r)-1536796800*L^4*Br*NT^2*ln(h)+11116875*Br*NB*NT*h^4+11116875*NB^2*Gr*R0^3*h-11116875*NB^2*Gr*R0*h^3+11116875*Br*NT^2*R0^3*h-11116875*Br*NT^2*R0*h^3+5336100*Br*NT^2*ln(h)*R0^4-5336100*Br*NT^2*ln(R0)*h^4+96049800*L^2*Br*NT^2*ln(R0)*h^2+5336100*Br*NT*R0^4*ln(r)*NB-96049800*L^2*Br*NT^2*ln(h)*h^2+96049800*L^2*Br*NT^2*ln(R0)*R0^2-1536796800*L^4*Br*NB*NT*ln(R0)-96049800*L^2*Br*NT^2*ln(h)*R0^2+24012450*Br*NT^2*ln(h)*R0^2*h^2+1536796800*L^4*Br*NB*NT*ln(h)+5336100*Br*NT^2*ln(h)*R0^3*h-24012450*Br*NT^2*ln(R0)*R0^2*h^2)*r^8/(L^6*NB*(ln(h)-ln(R0)))-(1/226609908940800)*(-5138546688*L^4*Br*NB*NT*ln(h)*R0-5138546688*L^4*Br*NB*NT*ln(h)*h+5138546688*L^4*Br*NB*NT*ln(R0)*R0+5138546688*L^4*Br*NB*NT*ln(R0)*h+5138546688*L^4*Br*NT^2*ln(h)*R0+5138546688*L^4*Br*NT^2*ln(h)*h-5138546688*L^4*Br*NT^2*ln(R0)*R0-5138546688*L^4*Br*NT^2*ln(R0)*h+5138546688*L^4*Gr*NB^2*ln(h)*R0+5138546688*L^4*Gr*NB^2*ln(h)*h-5138546688*L^4*Gr*NB^2*ln(R0)*R0-5138546688*L^4*Gr*NB^2*ln(R0)*h-5138546688*L^4*Gr*NB*NT*ln(h)*R0-5138546688*L^4*Gr*NB*NT*ln(h)*h+5138546688*L^4*Gr*NB*NT*ln(R0)*R0+5138546688*L^4*Gr*NB*NT*ln(R0)*h)*r^7/(L^6*NB*(ln(h)-ln(R0)))-(1/226609908940800)*(341510400*L^2*Gr*NB*NT*ln(R0)*R0*h^3+1536796800*L^2*Gr*NB*NT*ln(R0)*R0^2*h^2-1536796800*L^2*Gr*NB*NT*ln(h)*R0^2*h^2-341510400*L^2*Gr*NB*NT*ln(h)*R0^3*h-341510400*L^2*Gr*NB*NT*ln(r)*R0*h^3+341510400*L^2*Gr*NB*NT*ln(r)*R0^3*h-12294374400*L^4*Gr*NB*NT*ln(R0)*R0*h+12294374400*L^4*Gr*NB*NT*ln(h)*R0*h+341510400*L^2*Br*NB*NT*ln(R0)*R0*h^3+1536796800*L^2*Br*NB*NT*ln(R0)*R0^2*h^2-1536796800*L^2*Br*NB*NT*ln(h)*R0^2*h^2-341510400*L^2*Br*NB*NT*ln(h)*R0^3*h-341510400*L^2*Br*NB*NT*ln(r)*R0*h^3+341510400*L^2*Br*NB*NT*ln(r)*R0^3*h+341510400*L^2*Gr*NB*NT*ln(R0)*h^4-341510400*L^2*Gr*NB*NT*ln(h)*R0^4-341510400*L^2*Gr*NB*NT*ln(r)*h^4+341510400*L^2*Gr*NB*NT*ln(r)*R0^4-341510400*L^2*Gr*NB^2*ln(R0)*R0*h^3-12294374400*L^4*Br*NB*NT*ln(R0)*R0*h+12294374400*L^4*Br*NB*NT*ln(h)*R0*h-1536796800*L^2*Gr*NB^2*ln(R0)*R0^2*h^2+1536796800*L^2*Gr*NB^2*ln(h)*R0^2*h^2+341510400*L^2*Gr*NB^2*ln(h)*R0^3*h+341510400*L^2*Gr*NB^2*ln(r)*R0*h^3-341510400*L^2*Gr*NB^2*ln(r)*R0^3*h-6147187200*L^4*Gr*NB*NT*ln(R0)*h^2-6147187200*L^4*Gr*NB*NT*ln(R0)*R0^2+6147187200*L^4*Gr*NB*NT*ln(h)*h^2+6147187200*L^4*Gr*NB*NT*ln(h)*R0^2+12294374400*L^4*Gr*NB^2*ln(R0)*R0*h-12294374400*L^4*Gr*NB^2*ln(h)*R0*h-341510400*L^2*Br*NT^2*ln(R0)*R0*h^3-1536796800*L^2*Br*NT^2*ln(R0)*R0^2*h^2+1536796800*L^2*Br*NT^2*ln(h)*R0^2*h^2+341510400*L^2*Br*NT^2*ln(h)*R0^3*h+341510400*L^2*Br*NT^2*ln(r)*R0*h^3-341510400*L^2*Br*NT^2*ln(r)*R0^3*h+341510400*L^2*Br*NB*NT*ln(R0)*h^4-341510400*L^2*Br*NB*NT*ln(h)*R0^4-341510400*L^2*Br*NB*NT*ln(r)*h^4+341510400*L^2*Br*NB*NT*ln(r)*R0^4+12294374400*L^4*Br*NT^2*ln(R0)*R0*h-12294374400*L^4*Br*NT^2*ln(h)*R0*h-6147187200*L^4*Br*NB*NT*ln(R0)*h^2-6147187200*L^4*Br*NB*NT*ln(R0)*R0^2+6147187200*L^4*Br*NB*NT*ln(h)*h^2+6147187200*L^4*Br*NB*NT*ln(h)*R0^2+626102400*L^2*Br*NT^2*R0^3*h-626102400*L^2*Br*NT^2*R0*h^3-626102400*L^2*Br*NB*NT*R0^4+626102400*L^2*Br*NB*NT*h^4+626102400*L^2*Gr*NB*NT*h^4-626102400*L^2*Gr*NB^2*R0*h^3+626102400*L^2*Gr*NB^2*R0^3*h-626102400*L^2*Gr*NB*NT*R0^4-6147187200*L^4*Gr*NB^2*ln(h)*h^2+6147187200*L^4*Gr*NB^2*ln(R0)*R0^2-6147187200*L^4*Gr*NB^2*ln(h)*R0^2-341510400*L^2*Br*NT^2*ln(R0)*h^4+341510400*L^2*Br*NT^2*ln(r)*h^4+341510400*L^2*Br*NT^2*ln(h)*R0^4-626102400*L^2*Gr*NB*NT*R0^3*h+626102400*L^2*Gr*NB*NT*R0*h^3-626102400*L^2*Br*NB*NT*R0^3*h+626102400*L^2*Br*NB*NT*R0*h^3+6147187200*L^4*Br*NT^2*ln(R0)*h^2-341510400*L^2*Br*NT^2*ln(r)*R0^4+6147187200*L^4*Br*NT^2*ln(R0)*R0^2-6147187200*L^4*Br*NT^2*ln(h)*R0^2-6147187200*L^4*Br*NT^2*ln(h)*h^2-341510400*L^2*Gr*NB^2*ln(R0)*h^4-393419980800*c2*ln(h)*ln(r)*L^2*NB+393419980800*c2*ln(R0)*ln(r)*L^2*NB-341510400*L^2*Gr*NB^2*ln(r)*R0^4+341510400*L^2*Gr*NB^2*ln(r)*h^4+341510400*L^2*Gr*NB^2*ln(h)*R0^4+6147187200*L^4*Gr*NB^2*ln(R0)*h^2-626102400*L^2*Gr*NB^2*h^4+626102400*L^2*Br*NT^2*R0^4-626102400*L^2*Br*NT^2*h^4-393419980800*c1*ln(h)*L^2*NB+327849984000*c2*ln(h)*L^2*NB+98354995200*L^2*p*NB*ln(h)+393419980800*c1*ln(R0)*L^2*NB-327849984000*c2*ln(R0)*L^2*NB-98354995200*L^2*p*NB*ln(R0)+626102400*L^2*Gr*NB^2*R0^4)*r^6/(L^6*NB*(ln(h)-ln(R0)))-(1/226609908940800)*(-3540779827200*c3*L^2*NB*ln(h)+3540779827200*c3*L^2*NB*ln(R0)+12294374400*L^4*Br*NT^2*ln(h)*R0^3*h+55324684800*L^4*Br*NT^2*ln(h)*R0^2*h^2-55324684800*L^4*Br*NT^2*ln(R0)*R0^2*h^2-12294374400*L^4*Br*NT^2*ln(R0)*R0*h^3-12294374400*L^4*Gr*NB^2*ln(r)*R0^3*h+12294374400*L^4*Gr*NB^2*ln(r)*R0*h^3+12294374400*L^4*Gr*NB^2*ln(h)*R0^3*h+55324684800*L^4*Gr*NB^2*ln(h)*R0^2*h^2-55324684800*L^4*Gr*NB^2*ln(R0)*R0^2*h^2-12294374400*L^4*Gr*NB^2*ln(R0)*R0*h^3+12294374400*L^4*Gr*NB*NT*ln(r)*R0^4-12294374400*L^4*Gr*NB*NT*ln(r)*h^4-12294374400*L^4*Gr*NB*NT*ln(h)*R0^4+12294374400*L^4*Gr*NB*NT*ln(R0)*h^4+12294374400*L^4*Br*NB*NT*ln(r)*R0^4-12294374400*L^4*Br*NB*NT*ln(r)*h^4-12294374400*L^4*Br*NB*NT*ln(h)*R0^4+12294374400*L^4*Br*NB*NT*ln(R0)*h^4-12294374400*L^4*Br*NT^2*ln(r)*R0^3*h+12294374400*L^4*Br*NT^2*ln(r)*R0*h^3-18441561600*L^4*Gr*NB*NT*R0^3*h+18441561600*L^4*Gr*NB*NT*R0*h^3-18441561600*L^4*Br*NB*NT*R0^3*h+18441561600*L^4*Br*NB*NT*R0*h^3+12294374400*L^4*Br*NB*NT*ln(R0)*R0*h^3+12294374400*L^4*Gr*NB*NT*ln(r)*R0^3*h-12294374400*L^4*Gr*NB*NT*ln(r)*R0*h^3-12294374400*L^4*Gr*NB*NT*ln(h)*R0^3*h-55324684800*L^4*Gr*NB*NT*ln(h)*R0^2*h^2+55324684800*L^4*Gr*NB*NT*ln(R0)*R0^2*h^2+12294374400*L^4*Gr*NB*NT*ln(R0)*R0*h^3-18441561600*L^4*Gr*NB^2*R0*h^3+18441561600*L^4*Gr*NB*NT*h^4+18441561600*L^4*Br*NB*NT*h^4-18441561600*L^4*Br*NB*NT*R0^4-18441561600*L^4*Br*NT^2*R0*h^3+18441561600*L^4*Br*NT^2*R0^3*h-12294374400*L^4*Gr*NB^2*ln(R0)*h^4-14163119308800*L^4*c2*ln(r)*NB*ln(h)+14163119308800*L^4*c2*ln(r)*NB*ln(R0)-3540779827200*c4*ln(r)*L^2*NB*ln(h)+3540779827200*c4*ln(r)*L^2*NB*ln(R0)-18441561600*L^4*Gr*NB*NT*R0^4+18441561600*L^4*Gr*NB^2*R0^3*h-12294374400*L^4*Br*NT^2*ln(r)*R0^4+12294374400*L^4*Br*NT^2*ln(r)*h^4+12294374400*L^4*Br*NT^2*ln(h)*R0^4-12294374400*L^4*Br*NT^2*ln(R0)*h^4-12294374400*L^4*Gr*NB^2*ln(r)*R0^4+12294374400*L^4*Gr*NB^2*ln(r)*h^4+12294374400*L^4*Gr*NB^2*ln(h)*R0^4+55324684800*L^4*Br*NB*NT*ln(R0)*R0^2*h^2-55324684800*L^4*Br*NB*NT*ln(h)*R0^2*h^2-12294374400*L^4*Br*NB*NT*ln(h)*R0^3*h-12294374400*L^4*Br*NB*NT*ln(r)*R0*h^3+12294374400*L^4*Br*NB*NT*ln(r)*R0^3*h+18441561600*L^4*Gr*NB^2*R0^4-18441561600*L^4*Gr*NB^2*h^4+18441561600*L^4*Br*NT^2*R0^4-3540779827200*L^4*p*NB*ln(R0)+5311169740800*c4*ln(h)*L^2*NB-5311169740800*c4*ln(R0)*L^2*NB-18441561600*L^4*Br*NT^2*h^4-14163119308800*L^4*c1*NB*ln(h)+7081559654400*L^4*c2*ln(h)*NB+3540779827200*L^4*p*NB*ln(h)+14163119308800*L^4*c1*NB*ln(R0)-7081559654400*L^4*c2*ln(R0)*NB)*r^4/(L^6*NB*(ln(h)-ln(R0)))-(1/226609908940800)*(-226609908940800*L^6*c2*ln(r)*NB*ln(h)+226609908940800*L^6*c2*ln(r)*NB*ln(R0)-226609908940800*L^6*c1*NB*ln(h)+226609908940800*L^6*c1*NB*ln(R0)-56652477235200*L^4*c4*ln(r)*NB*ln(h)+56652477235200*L^4*c4*ln(r)*NB*ln(R0)-56652477235200*L^4*c3*NB*ln(h)+56652477235200*L^4*c4*ln(h)*NB+56652477235200*L^4*c3*NB*ln(R0)-56652477235200*L^4*c4*ln(R0)*NB)*r^2/(L^6*NB*(ln(h)-ln(R0)))-(1/226609908940800)*(-226609908940800*L^6*c4*ln(r)*NB*ln(h)+226609908940800*L^6*c4*ln(r)*NB*ln(R0)-226609908940800*L^6*c3*NB*ln(h)+226609908940800*L^6*c3*NB*ln(R0))/(L^6*NB*(ln(h)-ln(R0)))):

NULL


 

Download vel.mw

 

I am trying to input data via a data table. Have several problems here.

I used an array because I want the row and column numbers to start at 0.

1st When the table appers after that the document runs hediously slow as in a second or two to enter a digit or letters appear after typing.  Like something is absorbing the computer resources. But I have a fast machine.

2nd Any data I  enter to the table vanishes but does get stored.

3rd I tried to turn it all into a procedure but cant get that to work.
 

restart

with(DocumentTools:-Components)

[Button, CheckBox, CodeEditRegion, ComboBox, DataTable, Dial, Label, ListBox, MathContainer, Meter, Microphone, Plot, RadioButton, RotaryGauge, Shortcut, Slider, Speaker, State, TextArea, ToggleButton, VideoPlayer, VolumeGauge]

(1)

with(DocumentTools:-Layout)

[Cell, Column, DocumentBlock, Equation, Font, Group, Image, InlinePlot, Input, Output, Row, Section, Table, Textfield, Title, Worksheet]

(2)

with(DocumentTools)

[AddIcon, AddPalette, AddPaletteEntry, Components, ContentToString, CreateTask, Do, GetDocumentProperty, GetProperty, InsertContent, InsertTask, Layout, RemovePalette, RemovePaletteEntry, RemoveTask, Retrieve, RunWorksheet, SetDocumentProperty, SetProperty, Tabulate]

(3)

ary := Array(0 .. 3, 0 .. 3)

Array(%id = 18446746457454449478)

(4)

``

 

 

``

DT := DataTable(identity = "DataTable0", variable = 'ary', rowheader, columnheader, columnnames = [beta^0, beta, beta^2, beta^3], rownames = [alpha^0, alpha, alpha^2, alpha^3])

xml := Worksheet(Group(Input(Textfield(DT))))

DocumentTools:-InsertContent(xml)

PN1 := copy(ary, 0 .. (), 0 .. ())

Array(%id = 18446746457454471998)

(5)

Matrix(PN1)

Matrix(%id = 18446746457454477550)

(6)

PN1[0, 0]

6

(7)

BiPolyNum := proc (a := 4, b := 4) local ary, DT; description "Creates Bi Polynumbers"; ary := Array(0 .. a, 0 .. b); DT := DataTable(identity = "DataTable0", variable = 'ary', rowheader, columnheader, columnnames = [1, beta, beta^2, beta^3], rownames = [1, alpha, alpha^2, alpha^3]); DocumentTools:-InsertContent(xml); return copy(ary) end proc

proc (a := 4, b := 4) local ary, DT; description "Creates Bi Polynumbers"; ary := Array(0 .. a, 0 .. b); DT := DocumentTools:-Components:-DataTable(identity = "DataTable0", variable = 'ary', rowheader, columnheader, columnnames = [1, beta, beta^2, beta^3], rownames = [1, alpha, alpha^2, alpha^3]); DocumentTools:-InsertContent(xml); return copy(ary) end proc

(8)

``

``

f := BiPolyNum()

Array(%id = 18446746457454464046)

(9)

f

Array(%id = 18446746457454464046)

(10)

Matrix(f)

Matrix(%id = 18446746457454466582)

(11)

``


 

Download DataTable_Experiment.mw

Why is maple showing 1D math when evaluating?

I have a trigonometric equation that outputs with a solution in terms of _B1 which I want to remove.

restart: solve({7*cos(2*t)=7*cos(t)^2-5, t>=0, t<=2*Pi}, t, allsolutions, explicit);

output:

{t = arccos((1/7)*sqrt(14))},

{t = 2*Pi-arccos((1/7)*sqrt(14))},

{t = 2*arccos((1/7)*sqrt(14))*_B1-2*_B1*Pi+2*Pi*_Z1-arccos((1/7)*sqrt(14))+Pi}

Is there anyway to get rid of the _B1, or somehow evaluate it by a substitution?

 

Even numerically the answer still retains the _B1.

{t = 1.006853685}, {t = 5.276331623}, {t = -4.269477938*_B1+6.283185308*_Z1+2.134738969}

 

Also it would be nice to remove the _Z1 subscript too, as the domain of the equation is [0, 2pi].

I tried removing the 'AllSolutions' command , but then I am missing two solutions:

solve({7*cos(2*t)=7*cos(t)^2-5., t>=0 and t<=2*Pi}, t, Explicit);

 {t = 1.006853685}, {t = 2.134738969}

There should be 4 solutions in the domain [0, 2pi].

First 18 19 20 21 22 23 24 Last Page 20 of 40