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.. 
 

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 ?

CH10A.MWS

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

Hi, 

I have created a Maple Workbook in which two worksheets are included. One worksheet does all the work and is password protected, while the other allows the user to insert input parameters and to view results.  In the password protected worksheet, I return the required results using the following statement:

> return matrix(G),matrix(A),matrix(DR),figplot

where the last returned item, figplot, is a figure (display(..., size = [600,600])).

When the results are requested in the other worksheet, Maple stacks the results on the same line, and when the width is not enough, the results overflow into the line next to it. I have two questions:

1- How can I print the returned entities (matrix(G), matrix(A), matrix(DR), figplot) on seperate lines instead of having them stacked in the same line?

2- Although I am specifying the size of the plot, when returned, the plot is very small. 

I tried to to return the output in a column matrix as follows:

> return matrix(3,1,[matrix(G), [matrix(A),matrix(DR)], figplot])

and it worked but the column matrix brackets encompass the results (looks ugly) and the figure shown is still very small and not according to the required size.

Thank you very much for your help

Can anyone help me to find a solution to  for the attached partial differential equation ?

I want to find a general solution to a partial differential equation by assuming that I know one solution, called , and trying to find another solution  by assuming that the general solution in the form of . I want to restrict the second solution to be in the form of  so that it satisfies the PDE, and is a function of  times . The latter makes error as the maple identifies that the function  depends on only one variable . Could you please help me to find a solution for  in the form ?

Also, I have trouble with defining the operator Do in the attached file.  When it operates on , maple gives D. It is not clear for me that whether this derivative is with respect to  or . I need is to define Do in a way so that the derivatives are correctly taken with respect to different separate variables.

Thank you for your help,

Ahmad

PDEs_KnownSolutoin.mw

Hi all, we know Maple provided discrim method to find discriminant of a polynomial 

I want to write a similar method with independent variable is ,... my code is below

Some examples

Maple already similar method? If not have, we can improve performance it?

Thank you very much.