Maple 2017 Questions and Posts

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

I would like to have my calculus results in its simplified form; below are some of the scenarios that can explain this.

 

 

As you can see the results are much simplified and in reduced form compared to what Maple gives. This has been tested for all integrals with roots and integrals with trigonometry functions. Is there any workaround in Maple that I would use to get it like the textbook result.

Thanks

 

When I resave a worksheet that is already saved on the computer desktop, it changes the position of the icon on the desktop: it is moved to the most left-upper column-row possible. Why that? Can anyone else reproduce that odd behaviour?

A worksheet can be closed using the Windows standard shortcut Ctrl-F4. But the Maple program itself cannot be closed using the Windows standard standard shortcut Alt-F4. Why that? Can anyone else reproduce that odd behaviour?

PS: I am using Maple 2017.1

This might be considered nit-picking, but nonetheless I think there is an issue: The metric tensor and its inverse are types (0,2) and (2,0) tensors, respectively. When once contracted with each other, the result is the Kronecker delta, which is then (necessarily) a tensor of type (1,1). I am therefore surprised to find that in Maple, this Kronecker delta is implemented as a type (0,2) tensor, via the command KroneckerDelta:

KroneckerDelta[mu,nu];

I don't think this makes any proper sense. I think that such an object of type (0,2) is, in fact, the metric itself. On a similar note, the (mixed tensor type) objects g_[mu,~nu] and g_[~mu,nu] are actually both Kronecker deltas, because they correspond to the metric having one index raised by contraction with the inverse metric itself. But, nonetheless, the following is the case:

g_[ mu,~nu],
g_[~mu, nu];

Relatedly, consider the following single contraction of the metric with its own inverse:

SumOverRepeatedIndices(g_[mu,nu]*g_[~nu,~sigma]);

Although this sum is formally quite correct, I think it should be given as just a Kronecker delta (of the correct mixed tensor type, that is).

Does Maple have build-in function, which when given an expression that depends on x and y, will separate it to a product of two functions, one that depedns on x only and the other that depends on y only?

The input mathematical expression is known to be seperable.

For example, If the input is

((3*y + y^2)*3*x)/(x + sin(x))

Then I'd the Maple function to take the above and return list or set of two parts  {(3*y + y^2)   ,     3*x/(x + sin(x) } (if it can't separate it, it can return null).

The API can be something as  

 f,g = find_product_functions(expression,[x,y])

Something like this is used on determining for example if RHS of first order ODE is separable in order to solve it more easily.  collect() does not really work for this. So Maple allready does this internally in its ODE solver when it checks if ODE is separable or not. But is the function available for users?

 

subs.mw

hi...i have a problem with subs rule

please help me.

thanks

eq81a:=m*diff(w(x,t),t$2)+c*diff(w(x,t),t)+E*Is*diff(w(x,t),x$4)+(P-f[p]*cos(Omega*t)-E*A/(2*l)*int(diff(w(x,t),x)^2,x=0..l))*diff(w(x,t),x$2)=0;

m*(diff(diff(w(x, t), t), t))+c*(diff(w(x, t), t))+E*Is*(diff(diff(diff(diff(w(x, t), x), x), x), x))+(P-f[p]*cos(Omega*t)-(1/2)*E*A*(int((diff(w(x, t), x))^2, x = 0 .. l))/l)*(diff(diff(w(x, t), x), x)) = 0

(1)

 

bc81a:=B[1](w(0,t))=0,B[2](w(0,t))=0,B[3](w(l,t))=0,B[4](w(l,t))=0;

B[1](w(0, t)) = 0, B[2](w(0, t)) = 0, B[3](w(l, t)) = 0, B[4](w(l, t)) = 0

(2)

 

 

nondimRule1:=w(x,t)=l*w[n](x/l,t/T);

w(x, t) = l*w[n](x/l, t/T)

(3)

nondimRule2:=x=l*x[n],t=T*t[n],Omega=Omega[n]/T,P=P[n]*E*Is/l^2,f[p]=f[n]*E*Is/l^2,c=c[n]*sqrt(E*Is*m)/l^2,A=2*alpha*Is/l^2;

x = l*x[n], t = T*t[n], Omega = Omega[n]/T, P = P[n]*E*Is/l^2, f[p] = f[n]*E*Is/l^2, c = c[n]*(E*Is*m)^(1/2)/l^2, A = 2*alpha*Is/l^2

(4)

intRule1:=Int(D[1](w[n])(x[n],t[n])^2,l*x[n]=0..l)=l*int(D[1](w[n])(x[n],t[n])^2,x[n]=0..1);

Int((D[1](w[n]))(x[n], t[n])^2, l*x[n] = 0 .. l) = l*(int((D[1](w[n]))(x[n], t[n])^2, x[n] = 0 .. 1))

(5)

dropnRule:=w[n]=w,x[n]=x,t[n]=t,c[n]=c,P[n]=P,f[n]=f,Omega[n]=Omega;

w[n] = w, x[n] = x, t[n] = t, c[n] = c, P[n] = P, f[n] = f, Omega[n] = Omega

(6)

 

eq81b:=convert(expand(l^3/(E*Is)*subs(int=Int,nondimRule2,intRule1,dropnRule,value(subs(nondimRule1,eq81a)))),diff);

Error, invalid input: diff received T*t, which is not valid for its 2nd argument

 

 

``

TRule;=solve(coeff(lhs(eq81b),diff(w(x,t),t$2))=1,{T})[1];

TRule

 

Error, `=` unexpected

 

``


 

Download subs.mw

 

 

One of the things I like in Maple is that I can return a local symbol from a proc() in some expression and it will not "conflict" with same symbol in the global space and will show the same.

I just do not know how Maple manages to do this.

For example:

foo:=proc(n)
   local x;
   x^n;
end proc;

And now if I do

x:=99;
foo(3);

Will return  x^3. This is even thought I had defined x:=99; before the call.

So there is one global `x` with value 99 and the `x` in the expression returned `x^3` did not get confused with the global `x`. Yet they look the same.

How does Maple manages to do this? In Mathematica, it always return local symbols with $nnn assigned to them to distinguish them from global symbols. (attaches the Module ID). For example, in Mathematica the same example above gives

Notice that the `x` returned from a proc() look different from inside the Module. It is not the same as the x in the global space.

Maple seems to be able to do the same thing, but using the same looking symbol. So it must be keeping track of things internally? It must know that the x in x^3 is not the same x in the x:=99 ofcourse.

Any idea how Maple does this?

 

With Maple 2017,when the following codes are being executed in the document mode, the error information

"Error, attempting to assign to `RealDomain:-`^`` which is protected"

will be displayed and echoed for several times, which won't appear in the previous versions of Maple.

 

 

///////////////////////////////////

 

restart

with(RealDomain):

with(Grid):

a := 1:

Seq(1/a, i = 1 .. 2)

 

///////////////////////////////////

 

 

Is this a new bug in Maple 2017?

restart;
sol:=dsolve(diff(y(x),x)= x/(sqrt(x^2-16))*1/(2*y(x)),y(x));

Gives

But the solution can also be written as

I just do not know how to transform the first solution to the second simpler one. I tried:

sol:=map(x->rhs(x),[sol]);
simplify(sol,radical);

Also tried simplify(sol,sqrt); simplify(sol,radical,symbolic); simplify(sol,size);

The simpler solution can be found as follows

restart;
sol:=dsolve(diff(y(x),x)= x/(sqrt(x^2-16))*1/(2*y(x)),y(x),'implicit');

But the term in the middle above is

Therefore the solution is really

eq:=y(x)^2-sqrt(x^2-16)-_C1 = 0;
solve(eq,y(x));

Which gives

What command to simplify the long solution to the shorter one obtained from the implicit?

 

The solution by Maple below is correct, but non-the-less, a little strange.

restart;
dsolve(diff(y(x),x)=3*x^2*(y(x)^2+1),y(x));

Gives

Ofcourse 3*constant is still constant. But it is a little strange and have no reason for it to be there.  When I solve it by hand

What could made Maple put the 3 in there? Again, solution is 100% correct, but it could be simpler.

Maple 2017.1

 

 

  1. I have a set of equations described on different lines and I want them to be executed together and their answers to be displayed in seperate lines after that. How do I do that?
  2. This is particulalry important because, at times I will have to document multiple equations sequentially and later on have their answers reffere to. This way its more structured and I can group equations falling into certain category.

Please find the attached Minimum working example. Any advise is appreciated. File link: [MWE.mw]

 

Thanks

 

 

Hello all, 

 

This is my first time with Maple, I have been a student of Mathematica for 7 years. I purchased Maple to learn a new software and I have heard great things about it. I somehow dont feel the flexibility of Mathematica in Maple documentations. It seems to be a bit constrained and not very straight forward in some aspects. Please correct me if I am wrong and also point out to tutorials or documents that I should be looking at before nose diving into Maple.

worksheet example here: 

I have faced 2 simple problems which I think is a bug in some form, or I may be wrong. Please advise.

  1. How do I insert Equation 7 before Equation 6? The worksheet wont let me do it.
  2. Why are 'and' and 'in' bolded automatically in SECTION format?

the following code results in an error message: Error, (in forget) lexically scoped parameter out of context

If I click on this error message, it brings me to a page which I visited too often.

if I uncomment the irrelevant minimize command, the error message disappears.

How can I prevent this error without giving irrelevant commands?

kind regards,

Harry Garst

hello.i have a problem for solving this equation.i dont why my past post about this is deleted.!!!

please help me

thanks,,,

9.mw
 

restart:

A1:= 27159:  n:= 0.59:  A2:= 70941:  h0:= 3e-4:   
L:= 0.8:  dpx := -98100:  uc:= 0.007:  k:=2.7:

ODE:= (A3,y)->
   (h0^(n+1)*L/sqrt(n)*(A1*exp(sqrt(n)*y/L)-A2*exp(-sqrt(n)*y/L))/k+dpx*y*h0^(n+1)/k+A3*(h0)^n/k)^(1/n)
;

proc (A3, y) options operator, arrow; (h0^(n+1)*L*(A1*exp(sqrt(n)*y/L)-A2*exp(-sqrt(n)*y/L))/(sqrt(n)*k)+dpx*y*h0^(n+1)/k+A3*h0^n/k)^(1/n) end proc

(1)

ODEINT:= proc(A3)
option remember;
local y;
   evalf(Int(ODE(A3,y), y= 0..1, epsilon= 1e-7)) - uc
end proc:

ReINT:= proc(A3x, A3y)
   Digits:= 15:
   Re(ODEINT(A3x + I*A3y))
end proc:

ImINT:= subs(Re= Im, eval(ReINT)):

Digits:= 7:
a3:= fsolve([ReINT, ImINT]);

fsolve([ReINT, ImINT])

(2)

A3:= Complex(a3[]);

Complex(fsolve([ReINT, ImINT])[])

(3)

Solve as IVP:

Digits:= 15:
sol:= dsolve({diff(u(y),y) = ODE(A3,y), u(0)=0}, numeric, range=0..1,  output=listprocedure):

Warning,  computation interrupted

 

NULL

``

NULL

NULL

plots:-odeplot(
   sol, [[y, Re(u(y))], [y, Im(u(y))]], y= 0..1,
   legend= [real, imag], labels= [y, u(y)]
);

Verify that boundary condition at u(1) is satisfied:

 

 

 

abs(eval(u(y), sol(1)) - uc);

sol(.5);

"\"

fy3 := eval(u(y), sol); with(CurveFitting); fy33 := PolynomialInterpolation([[0, fy3(0)], [.1, fy3(.1)], [.2, fy3(.2)], [.3, fy3(.3)], [.4, fy3(.4)], [.5, fy3(.5)], [.6, fy3(.6)], [.7, fy3(.7)], [.8, fy3(.8)], [.9, fy3(.9)], [1, fy3(1)]], y)

DEBI := int(fy33, y = 0 .. 1)

NULL

``

plot(DEBI, y = 0 .. 1)

``

``

``

``

``


 

Download 9.mw

 

Dear sir, 

 

I request your kind help with regard to the above cited subject.

 

 

With warm regards.

 

Mr.M.Anand

Associate Professor in Mathematics.

Aurora's Technological & Research Institute, 

Hyderabad-INDIA

 To Whom It Might Concern,

I purchased Maple on Jun 5, 2017 for my Window laptop. However, my harddrive in my laptop was broken and i lost all my stuff including Maple. I am using my husband's laptop right now. How can re install Maple? I still have the transaction number (removed). Would you let me know what I can do to re-install the software? Also, my husband's laptop is a Mac. 

I am looking forward to hearing from you.

Nancy Phu.

First 32 33 34 35 36 37 38 Page 34 of 40