Unanswered Questions

This page lists MaplePrimes questions that have not yet received an answer

Hello,

a follow up question.

I am solving some overdetermined system of ODEs in cylindrical coordinates r,phi,Z. I obtain some equations of the following type:

(diff(_F1(phi, Z), phi)*r + diff(diff(s_r(phi, Z), phi), phi))/r = -s_r(phi, Z)/r

As can be seen, the differentiated functions do not depend on r, which is an independent variable. Thus, the correct solution is to separate the equation and have 

_F1(phi,Z)=_F1(Z), s_r(phi,Z)=s_r(Z).

By using dsolve, I always obtain a solution containing r.

A similar problem that does no contain derivatives is solved by solve/identity.

Is there something similar for dsolve?

EDIT: I again put here more info and file. I solve some overdetermined system of differential equations.

[diff(s_r(r, phi, Z), r) = 0, diff(s_r(r, phi, Z), phi) = -diff(s_phi(r, phi, Z), r)*r^2,
diff(s_Z(r, phi, Z), r) = -diff(s_r(r, phi, Z), Z), diff(s_phi(r, phi, Z), phi) = -s_r(r, phi, Z)/r, 
diff(s_Z(r, phi, Z), phi) = -diff(s_phi(r, phi, Z), Z)*r^2, diff(s_Z(r, phi, Z), Z) = 0, 
diff(m(r, phi, Z), r) = s_Z(r, phi, Z)*B_phi(r, phi, Z) - s_phi(r, phi, Z)*B_Z(r, phi, Z),
 diff(m(r, phi, Z), phi) = s_r(r, phi, Z)*B_Z(r, phi, Z) - s_Z(r, phi, Z)*B_r(r, phi, Z),
 diff(m(r, phi, Z), Z) = s_phi(r, phi, Z)*B_r(r, phi, Z) - s_r(r, phi, Z)*B_phi(r, phi, Z), 
s_r(r, phi, Z)*diff(W(r, phi, Z), r) + s_phi(r, phi, Z)*diff(W(r, phi, Z), phi) + 
s_Z(r, phi, Z)*diff(W(r, phi, Z), Z) = 0]

After some time, I arrive at the equation in the original question. So the independance of the other functions on r is the consequence of the other equations.

Here is the file (shortened): mwquestion2.mw

Sorry if it's a dumb question, i'm new with Sage.

I'm doing the multiplication of two matrix and I obtain the expression 

 

[142√(2x−2√)+122√y] 

I want to represent this on 3d, but I realized that the expression has to be written with operators '*'

How can I add the operators to my expression?

 

   1 1 1 2 2 3 3  

(  1 3 3 2 2 1 2 )

how to input this permutation into permgroup ?

 

Who deleted www.mapleprimes.com/questions/228364-Graphics-Problem

and www.mapleprimes.com/posts/211539-Buggy-Lights

I thought I had alzheimers or something.  Who's deleting these posts and questions.  I realize they are similar but there is different information in there that was useful.  Why don't we delete all similar questions .. of course this would be horrible since some have slight differences - granted some don't have answers tagged on the end and are exact copies those are of course waranted and I would agree to be deleted becuase they were duplicates.

Please undelete those questions - there were referenses and extensive conversations in there that are not present in the saving question.

This is another problem I just found in Maple 2019.2 on windows 10. professional.

I wanted to close Maple, so did  File->Exit 

 

But Maple did nothing. It did not close.  Also Alt-F4 did not close Maple. I had to click on the little X on top right corner of the open window to close Maple.  

In earlier version this used to work to close Maple.

Do others see this as well?  To reproduce, simply start Maple, and do File->Exit.

Here is a movie also

 

This may be a bug.  In Maple 2019.0

plot3d(x^2+y^2)

Now grab the graphing window and drag to adjust it's size.  It doesn't adjust it's size until you let go of the mouse button.

If maple can auto discover derived equations, 

then most variable are unknown name which may or may not be physics variable.

so, how to guess which physics data suitable for these look like meaningless variable?

or these tools only visualize the relationship of existing known equations?

How to solve this DE with IC by using DTM.
D^m u(x,t)=u''(x,t)-u^2 (x,t), where n-1< m < n 

IC: u(x,0)=1+sin(x), and u'(x,0)=0

I have a problem for school that I need help with. 

Solve: by variation of parameters.

x3y’’’ – x2y’’ − 2xy’ + 6y = x2

Use any software, e.g. Maple, as an aid in computing roots of the auxiliary equation and the Wronskian-based determinants W, W1, W2, W3

I was able to hand calculate the roots as [m1 = 3, m2 = 2, m3 = -1]
 

NULL

with(VectorCalculus); Wronskian([x^3, x^2, 1/x], x, 'determinant')

Matrix(%id = 18446744078324919462), -12*x

(1)

with(VectorCalculus); Wronskian([0, x^2, 1/x], x, 'determinant')

Matrix(%id = 18446744078349556782), 0

(2)

with(VectorCalculus); Wronskian([x^3, 0, 1/x], x, 'determinant')

Matrix(%id = 18446744078266878782), 0

(3)

with(VectorCalculus); Wronskian([x^3, x^2, 0], x, 'determinant')

Matrix(%id = 18446744078331625094), 0

(4)

m3 := Matrix(3, 3, [0, x^2, 1/x, 0, 2*x, 1/x^2, x^2, 2, 2/x^3])

Matrix(%id = 18446744078385724822)

(5)

Determinant(m3, method = multivar)

Determinant(Matrix(%id = 18446744078385724822), method = multivar)``

(6)

roots(x^3-x^2-2*x+6, alpha, beta)

[]

(7)

``

solve(x^3-x^2-2*x+6 = 0, x)

-(1/3)*(71+9*58^(1/2))^(1/3)-(7/3)/(71+9*58^(1/2))^(1/3)+1/3, (1/6)*(71+9*58^(1/2))^(1/3)+(7/6)/(71+9*58^(1/2))^(1/3)+1/3+((1/2)*I)*3^(1/2)*(-(1/3)*(71+9*58^(1/2))^(1/3)+(7/3)/(71+9*58^(1/2))^(1/3)), (1/6)*(71+9*58^(1/2))^(1/3)+(7/6)/(71+9*58^(1/2))^(1/3)+1/3-((1/2)*I)*3^(1/2)*(-(1/3)*(71+9*58^(1/2))^(1/3)+(7/3)/(71+9*58^(1/2))^(1/3))

(8)

m3 := Matrix(3, 3, [0, x^2, 1/x, 0, 2*x, 1/x^2, 1/x, 2, 2/x^3])

Matrix(%id = 18446744078348970214)

(9)

Determinant(m3, method = multivar)

Determinant(Matrix(%id = 18446744078348970214), method = multivar)

(10)

``

m3 := Matrix(3, 3, [0, x^2, 1/x, 0, 2*x, 1/x^2, 1/x, 2, 2/x^3])

``

``


 

Download project_27.mw

 

 

 

 

I have just upgraded my laptop from Windows 7 to Windows 10.  On starting up Maple 2018,  I receive the attached message on screen.  This is after previously loading the worksheet successfully.   Today,  I am not able to do so.  I need to permanently register my firewall to allow Maple to run; can anyone help?

Melvin

 

After upgrading to 2019.2 (on windows 10, 64 bit) using the pull down menu option I noticed the following problem.

I create a new worksheet and run it. Here is an example

restart;

for i from 1 to 10 do
    print(i);
od;

(the restart is on seperate cell).

It runs OK first time. I close Maple. Restart Maple again, open the same worksheet. Now once I evaluate the restart; command, the kernel crashes. This happens each time.




Any suggestions what to do?  Only fix I found is to copy the content of the worksheet and create a new worksheet, then it runs. But then the new worksheet fails again when I close Maple and restart again. So this is not a solution.

I also noticed it crashes ONLY when running the restart; command/cell. Other commands work OK. So something to do with initialization. I have not changed anything at all on my end after upgrading to 2019.2 from 2019.1.

Evaluating the next cell first thing when opening the worksheet, then it runs OK and no crash.

This only started to happen after updating to Maple 2019.2. No problem with 2019.1. Clearly this version was not well tested.  How could one revert back to 2019.1 ??

Download bug.mw

To answer comment below. It also happens in document mode:

 

I have a .mapleinit file that amongst other things sets libname so my own packages are accessible and can be loaded using with(). Imagine my surprise when I found that a maple program that ran half an hour ago on Maple2019 bombed when it did not find its package. Closer investigation indicates that maybe .mapleinit is only searched for in currentdir()??

The updated Maple is `Maple 2019.2, APPLE UNIVERSAL OSX, Oct 30 2019, Build ID 1430966`

The old one is gone :-(. It was Maple2019.0. I should add that on this system, Maple is installed on an administrator account that is not my user account (luckily I have access to that). I am running macOS 10.14.6 Mojave.

I verified my older Maples still work as before.

Has anyone else seen this behaviour?

Mac Dude.

Edit: I need to be more specific here: When I double-click a worksheet (.mw file) then currentdir() will be set to the directory the .mw file is in. This has been long-standing Maple behaviour. But then currentdir() does not cover the home directory where .mapleinit sits. As a result, a blank sheet finds my .mapleinit whereas my files (that habitually I open with double-click or drag-drop) do no longer run the .mapleinit file.

 

can anyone help me to calculate the exact  value of the eigenvalues of this matrix:
 

Download mat.mw

 

is there any library or tools to design index of Grassmannian and its k and n for Schubert use?

is there any library to relate poset with index of Grassmannian and its k and n for Schubert use

with(RegularChains);
R := PolynomialRing([x, y, z]);
p1 := x^2+5-2*x*z;
p2 := z^3*y+x*y^2;
p3 := -8*z^3+3*y^2;
F := [p1, p2, p3];
MainVariable(p1, PolynomialRing([x, y, z]));
MainVariable(p1, PolynomialRing([z, y, x]));
prem(p1, p2, x);
prem(p2, p1, z);
prem(p3, p2, y);
premcustom := proc(Fparam,Gparam, xparam)
local R, G, F, lcg, lcr, dr, dg:
R := Fparam:
G := Gparam:
F := Fparam:
if degree(G,xparam) = 0 then
 print("return 0"):
 return 0:
elif degree(F, xparam) < degree(G, xparam) then
 print("return R"):
 return R:
else
 lcg := coeff(G, xparam, degree(G, xparam)):
 dg := degree(G, xparam):
 while degree(R, xparam) > degree(G, xparam) do
  lcr := coeff(R, xparam, degree(R, xparam)):
  dr := degree(R, xparam):
  R := lcg * R - lcr * G * (x^(dr - dg)):
 od:
end if:
return R:
end proc:
coeff(p1, x, degree(p1, x));
coeff(p2, x, degree(p2, x));
prem(p1,p2,x);
prem(p2,p1,z);

premcustom(p1,p2,x);
premcustom(p2,p1,z);
 
why premcustom looping?
 
is this coeff(p1, x, degree(p1, x)); wrong ?
 
should this to get real coefficient in number instead of variable and how?
 
is there equivalent function as dprem?
 
with(diffalg):with(diffalg):
FlessThanG := proc(Fparam, Gparam, PRing)
F := Fparam:
G := Gparam:
if leader(F,PRing) < leader(G,PRing) or (leader(F,PRing) = leader(G,PRing) and degree(F, leader(F,PRing)) < degree(G, leader(G,PRing))) then
 return True:
else
 return False:
end if:
end proc:
dprem := proc(Fparam, Gparam, x, PRing)
local R, theta, thetax, thetaG, F, G;
F := Fparam:
G := Gparam:
R := F:
thetax := leader(G, PRing);
while FlessThanG(R, G, PRing) or FlessThanG(G, R, PRing) do
theta := leader(R, PRing);
R := prem(R, thetaG, theta)
od:
return R:
end proc:
R := differential_ring(ranking = [[x,y,r]], derivations = [t], field_of_constants=[m,l], notation = diff):
p1 := m*x[2] + r*x;
p2 := m*y[2] + r*y - g;
p3 := x^2 + y^2 - l^2;
dprem(p1, p3, x, R);
 
 
is x[2] = diff(x,t$2) in diffalg ?
 
if not how to write in this way?
 
i had defined constant
but got error
 
Error, (in DifferentialAlgebra:-Tools:-LeadingDerivative) unknown symbol (approx. error location: [m*x[2 --> ] <-- +r*x])
 
 
5 6 7 8 9 10 11 Last Page 7 of 259