Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I am learning Typesetting so I can get better Latex.

I'd like to have all derivatives generated as y'(x) and not dy/dx regardless of what the derivative variable is. (It is more clear this way).

I found I can't even give Maple a list of the variables. But must do it each time one at a time.

I have thousands of ODE's I need to typeset, (they use different variables (depending on the textbook), and not always x or t, but the set of variables is not large).

so having to check and keep changing this setting each time is awkward (but I think is doable, as I know what the independent variable is for each ode so I can call Typesetting at start of each ode processing, once I know what the variable is).   

It wil be much easier if there is a way to tell Maple to do this for all variables. May be there is a good reason why Maple does not seem to allow such an option?

Here is an example

Typesetting:-Settings(usedot=false,prime=t,typesetprime=true):
ode:=diff(y(t),t$2)+diff(y(t),t)+y(t)= 0:
Physics:-Latex(ode);

Which gives the Latex I want

But if the variable was x instead of t, I have to do this

Typesetting:-Settings(usedot=false,prime=x,typesetprime=true):
ode:=diff(y(x),x$2)+diff(y(x),x)+y(x)= 0:
Physics:-Latex(ode);

I tried doing prime=[x,t] but Maple did not like this.

Is there a way or trick to tell Maple to use typesetprime=true for any letter?  I looked at https://fr.maplesoft.com/support/help/Maple/view.aspx?path=Typesetting%2fSettings but do not see a way so far.

Maple 2020.1 with Physics 776

I am trying the new Physics:-Latex command in version 774.

Couple of isssues I found:

1) One test I have failed with this error

Error, (in anonymous procedure called from anonymous procedure called from anonymous procedure called from anonymous procedure called from Physics:-*) invalid boolean expression: NULL
 

Using the normal latex() command works OK. I hope it is not too hard to fix since this new version seems to have fixed the fraction problem in latex(). Please see

restart;
sol:=dsolve(t*(t-2)^2*diff(diff(y(t),t),t)+t*diff(y(t),t)+y(t) = 0,y(t));
latex(sol);  #OK

Physics:-Latex(sol); #error

 

2) Maple now generates Latex using maple macro, which I do not know what it is.  Comparing

restart;
ode:=t*(t-2)^2*diff(diff(y(t),t),t)+t*diff(y(t),t)+y(t) = 0;
latex(ode)

which gives 

t \left( t-2 \right) ^{2}{\frac {{\rm d}^{2}}{{\rm d}{t}^{2}}}y
 \left( t \right) +t{\frac {\rm d}{{\rm d}t}}y \left( t \right) +y
 \left( t \right) =0

compare to

Physics:-Latex(ode)

which gives

t \left(t -2\right)^{2} \Mapleoverset{\ldots}{y}\left(t \right)
      +t \Mapleoverset{.}{y}\left(t \right)+y \left(t \right) = 0

What is \Mapleoverset ? googling around, these look like part of internal Maple style sheet for Latex?

Does this mean when using the new Physics Latex command, now one needs to use the package maplestd2e.sty?  With the normal latex() command, this was not needed. 

Just trying to understand the new setup.

What are the things you most like to see improved/add to next version of Maple? 

This is my list for a starter:

1.  Improve the debugger. Debugger is very useful but needs more work. At least be able to see code listing in larger view as one steps in for example. See Matlab debugger for inspiration.

2.  Improve Latex. It still does not do fractions well. Posted about this before.

3. Eliminate hangs when using timelimit(). On long runs, random hangs happen when timelimit() do not expire as requested. Posted about this before.

 

 

I'm having hard time to do the following substitution. Please help

diff(x(t), t, t) := -R*diff(psi(t), t)*diff(phi(t), t)*sin(phi(t)) + R*diff(psi(t), t, t)*cos(phi(t)) <---- this is my expression xdoubledot

I want to substitute it in following equation

`Eq__&lambda;1` := diff(x(t), t, t) - R*diff(psi(t), t, t)*cos(phi(t)) + R*diff(psi(t), t)*diff(phi(t), t)*sin(phi(t))  <----this is my equation

 

I have turned on physics and algebra mobdules. it still not substitiutes. even though you see diff here I have dots in my maple output I don't know why it changes when I copy paste here.

I tried subs, solve,simplify commands. None of them works

Thank you for your help in advance

 

Hello,

I'm new to Maple and I'm struggling with dsolve output.

I solve this equation
 

ode:=diff(y(t),t)=-(4/3)*(a/eta)*(y(t)^3+1):
lambdaE:=dsolve({ode,y(0)=1},y(t),numeric,method=classical[rk4],parameters=[a,eta],output=listprocedure):
lambdaE(parameters=[theta[5,1],eta0]):

tCompute:=20*(lambda-~1)/(50/60): #vector with times I need 
# lambda is a 371 lines column vector

Then I need the value of the solution at each time to compute components of a matrix

J:=Matrix(371,5):
for i from 1 to 371 do
	J[i,1]:=2*(lambda[i]-1/lambda[i]^2);
	J[i,2]:=4*(lambda[i]^2-2/lambda[i]-3)*(lambda[i]-1/lambda[i]^2);
	J[i,3]:=6*(lambda[i]^2-2/lambda[i]-3)^2*(lambda[i]-1/lambda[i]^2);
	J[i,4]:=2*(1-1/lambda[i]^3);
	J[i,5]:=2*(lambdaE(tCompute[i])^2/lambda[i]-(1/lambdaE(tCompute[i])*lambda[i])); #here is where I need the value
end do;

However the output in the 5th column of J is something with y(t) and t inside

1.997498088*[t = 0.0300605520000000, y(t) = 0.999829465489350]^2 - 2.002505046/([t = 0.0300605520000000, y(t) = 0.999829465489350])

Would you have a solution to get only the numerical value ?

Thanks a lot in advance

Manon

This is the same example in previous post,

But here my question is how to produce smother plot of date. I tried the interpolation package but it not working.. 
 

restart:

with(LinearAlgebra): with(CurveFitting): with(Interpolation): with(plots): with(Statistics):
Omega:=0.1:no:=0:Delta:=0:tau0:=1:c:=0:gamma1:=0:
j:=1:k1:=0:t:=1:xr:=1:Gamma:=0.1:

Error, invalid input: with expects its 1st argument, pname, to be of type {`module`, package}, but received Interpolation

 


for i from -10 to 10 do
d[i]:=i:
b:=sum(-gamma1*tau0+I*tau0*Delta-2*a*(k-1)*xr,k=0..1):
a:=1+I*c;
c1:=sqrt(conjugate(a))/tau0:
c2:=0.5*((conjugate(b)/sqrt(conjugate(a)))):
lambda1[i]:=2*(Gamma-I*d[i])*(t-k1)/j+(Gamma-I*d[i])^2:
lambda2[i]:=t*(1-j)-k1/j+1/sqrt(2)*(Gamma-I*d[i]):
lambda3[i]:=c1*(t-k1)/j+c1*(Gamma-I*d[i])+c2:
J1[i]:=sqrt(Pi)/sqrt(2)*(1-erf(lambda2[i]));
J1mod[i]:=(Re(J1[i]))^2+(Im(J1[i]))^2:

g1[i]:=0.5*sqrt(Pi)*tau0*exp(c2^2)
    *exp(-conjugate(a)*((k-1)*xr)^2)/(sqrt(conjugate(a)));

F2[i]:=map(fnormal,evalf(sqrt(2)*Int(unapply(exp(-x^2)*erf(-sqrt(2)*c1*x+lambda3[i]),x),-lambda2[i]..5,epsilon=1e-5)));

J2[i]:=sum(g1[i]*J1[i]+g1[i]*F2[i],k=1..1);

J2mod[i]:=abs(J2[i])^2:
h1[i]:=map(fnormal,evalf(exp(-d[i]^2))):
f[d[i]]:=map(fnormal, evalf(h1[i]*(no*J1mod[i]+Omega^2*J2mod[i]))):
#f[d[i]]:=-d[i]^2*log(exp(1))+2*log(Omega)+log(J2mod[i]):
#map(simplify~, %[2], zero):
end do:

with(CurveFitting):
ptsN1 := CodeTools:-Usage(
           [seq([d[i], evalf(f[d[i]])], i = -10 .. 10)]
                         ):

memory used=3.44KiB, alloc change=0 bytes, cpu time=0ns, real time=0ns, gc time=0ns

 

 

P1:=plot(ptsN1, color = black, linestyle = solid, thickness = 2,
     color = black, axes = boxed, labelfont = ["HELVETICA", 14],
     labels = ["d", "S_M(d)"]);

 

         
 

points := [seq(x, x = -10 .. 10, 1)]:
data   := [seq(evalf(f[d[i]]), i = -10 .. 10, 1)]:
f := Interpolate(points, data):
plot(f,-10..10);

Error, (in plot) procedure expected, as range contains no plotting variable

 

``

``


I noticed this package is not introduced in Maple 17.. what is the alternative way to produce smother curve

Download interploation.mw

 

y(t)/diff(y(t),t)

 

how to dsolve when highest derivative at denominator?


For a from 1 to 2 do

convert(cat("k",a),symbol) := a;

od;

Error : can not do assignment 

how to do ?

Hello

I need Nu[a] in the label of y-axis.

i am writing      labels = [eta, 'Nu[a]'(eta)] but I am not getting. Please help me for writing correct code.

 

Hello

I am trying to reproduce a simple result from Linear Control Theory using Laplace Transform.  Here are the steps:

assume(omega>0);assume(zeta>0 and zeta<1);
tf:=omega^2/(s^2+2*zeta*omega*s+omega^2);
y:=tf*1/s;
yt:=inttrans:-invlaplace(y,s,t);
dyt:=diff(yt,t);

So far, so good.  Now I need to find the values of t (time) such that the dyt=0.  

 

solve(dyt=0,t);

but Maple returns only the trivial solution, that is, zero.   How can I find the next one (Tp - peak time = Pi/(omega*sqrt(1-zetaˆ2)))?  

 

Many thanks

 

Dear All, 

When I checked the dynamic responses of second-order systems, I saw this:

While MATLAB step response shows a smooth curve (the blue one), Maple came up with wiggles (the blue one).

What might have worked on Maple's side to make the curve wiggling?

The Maple worksheet as well as MATLAB M script are attached. The MATLAB M script extension (.m) is not allowed, therefore it was modified to '.txt'. 

Sallen_Key.mw

SallenKey_check.txt

We have just released an update to Maple, Maple 2020.1.1. This update includes the following:

  • Correction to a problem that occurred when printing or exporting documents to PDF. If the document included a 3-D plot, nearby text was sometimes missing from the printed/exported document.
  • Correction to an issue that prevented users from installing between-release updates to the Physics package

This update is available through Tools>Check for Updates in Maple, and is also available from our website on the Maple 2020.1.1 download page. If you are a MapleSim user, you can obtain this update from the corresponding MapleSim menu or MapleSim 2020.1.1  download page.

In particular, please note that this update fixes the problems reported on MaplePrimes in Maple 2020.1 issue with print to PDF  and installing the Maplesoft Physics updates. As always, we appreciate the feedback.

Hi,

Thank you all for participating to my questions before.

I was wondering if we can change an equation into quartic form or quadratic form in maple.

I have succesfully got an equation like the following:

`S__2 ` = sin(alpha - phi)*sin(-beta + alpha)*(gamma*H^2*sin(beta - varepsilon)*sin(alpha) - h^2*sin(beta)*sin(alpha - varepsilon)*gamma + h^2*sin(beta)*sin(alpha - varepsilon)*psi)/(2*sin(-beta - delta - phi + alpha)*sin(beta)^2*sin(alpha - varepsilon)*sin(alpha)) - S__1

From the paper I studied, they change it to form an equation like this

 

`S__2 ` = 1/2*gamma*H^2*sin(beta - alpha)*(M__3(X)^3 + M__2(X)^2 + M__1*X)/(sin(beta)^2*(D__3(X)^3 + D__2(X)^2 + D__1*X + D__0)) - S__1

 

Where value of M, D and X is sets of long equation. Can someone teach me how to assign maple to change this kind of equation to another form of equation. It's good enough if I can learn how to learn the basic.

 

Thank you.

 

Regards

Faiz Farhan

First 10 11 12 13 14 15 16 Last Page 12 of 1728