Hello to all

I am Jorge Gracia an exchange student working on my Thesis.

I started working with maple recently, and I am having some problems to achieve the results that I am looking for.

My work consist in the next Steps.

Using maple and taking some fatigue tests data, create a wav file with the some of the following formats:

Supported data formats are wave type 1 (PCM integer Data 8 Bit unsigned, 16 Bit signed or 32 Bit

signed) and wave type 3 (PCM float data (32 Bit IEEE float).

In the same way, if I already have a wav file, read the values of every point of the wave in order to understand and study the results.

If I have a normal plot, convert it into a wav file plot with the formats from above

With all the results, I will be able to work with this wav files in a fatigue test machine achieving more realistic results than with the standard waves.

All the information about this will be welcome.

Thank you for all.

Hello,

Thanks to the function Exports, it is possible to see the differents functions and modules which are in a module.

For example, I find in the module GetMultibody a module called GetLGM which permits to conduct some linear graph analysis. I'm interested by this module.

However, to my knowledge, I didn't see any help on this module.

How can I do so as to understand the inputs to give to the functions of this module in order to use the functions of this module ? In other words, how it is possible to use functions from a module when you don't have any help at your disposal ? 

Thanks a lot for your help.

Dear all;

I need you help for solving this problem, and thanks in advantage for your help.

I have a polynom like  P =x^6-4*x^3+x-2;  and i would like to find an approximate value of the roots in some interval [a,b] =[-2,2] using sturm sequence.

The method is based on:

1) first construct the sturm sequence:

For given polynom P =x^6-4*x^3+x-2;

Let S0=P;

S1=diff(p,x);

let   s:=quo(S0,S1,x);
       S2:=-rem(S0,S1,x);

.... S[k+1-rem(S[k-1],S[k]);

S[k] is the sturm sequence.

2) let f(a)= number of change of sign in the sturm sequence and f(b) the same . so f(b)-f(a) give the number of roots in the interval [a,b].

3) If f(b)-f(a) =0 so there are no roots

and if f(a)-f(b)=1 one can find the root

4) if f(a) -f(b) >2  :

given toterance tol=0.001; for example

if the abs(a-b)<2*epsilon we display a message that there are k roots at (b+a)/2

with our error tolerance

5) otherwise if c=(b+a)/2 is not a root of P_k(x)  for any k, ( where p_k is an element of the sturm sequence ) 

we divide the interval into equal halves [a,c] and [x,b] and we run step 2 on each interval

else if c is a root of one of these p_k(x) add any time account to c so that c lies close the middle of [a,b] and not a root

6) Give all the roots ( approximate the rrots with small error epsilon).

I kindly  appreciate your help

Hi,

Im trying to study some questions and I'm using maple to verify my answers.

Theres a few polynomial factoring questions and linear equation questions Im trying to get

maple to show its solutions steps using showsolution() no matter where I put it  the function wont work.

Ive switched between math/text functions. Im still pretty new to maple but I can't find any information on how to do it

on the web/youtube.

Thanks in advance!

For problem simplify(abs(1-b)+abs(1+b)), I want maple to take out the abs and get results for different ranges of b.

I have double indexed functiions f[j,k] of one variable and double indexed coefficients a[j,k].

I want my print do look like a[1,1]f₁₁+a[1,2]f₁₂ that is, the values of a[j,k] should appear beside the functions' names, like

7f₁₁+2f₁₂-3f₂₁ etc.

Thank yopu for any help

Hi, 

     I'm trying to solve this PDE, and Maple 2015 gives me a solution quickly. I can test the solution with pdetest() and this verifies that it works. However, when I try to verify this myself I don't get zero. Is there some trick pdetest() is using to that I am missing? Or is pdetest() wrong in this case?

Download PDESolving.mw

Hello,

I'm writing to ask how to equalize the coefficients of two multivariate polynomials. In particluar, I have two polynomials whose arguments are ln(E),ln(K),ln(L) (their levels, squared levels and interaction terms). The first one is:

(1/2*(p*a*b+(g-p)*b-g))*b*v*a*ln(E)^2-(-1+b)*v*(g-p+a*p)*b*a*ln(E)*ln(K)-b*p*(a-1)*v*a*ln(E)*ln(L)+v*a*b*ln(E)+(1/2*(p*(-1+b)*a+(g-p)*b+p))*(-1+b)*v*a*ln(K)^2+(-1+b)*v*p*(a-1)*a*ln(K)*ln(L)-v*a*(-1+b)*ln(K)+(1/2)*a*p*v*(a-1)*ln(L)^2-v*(a-1)*ln(L)

the second one is:

x_1*ln(E)+x_11*ln(E)^2+x_12*ln(E)*ln(K)+x_13*ln(E)*ln(L)+x_2*ln(K)+x_22*ln(K)^2+x_23*ln(K)*ln(L)+x_3*ln(L)`+x_33*ln(L)^2

I would like to know if it is possible to equalize the coefficients of the two polynomials and find the following system:

v*a*b = x_1, -v*(a-1) x_3, -v*a*(-1+b) = x_2, a*b*v*(b*rho*a-b*rho+g*(-1+b)) = x_11, v*rho*a*(a-1) = x_33, v*a*(rho*(-1+b)*a-rho*(-1+b)+b*g)*(-1+b) = x_22, -a*v*rho*(a-1)*b = x_13, -a*v*(a*rho-rho*u+g)*b*(-1+b) = x_12, a*v*u*rho*(a-1)*(-1+b) = x_23

I tried using "coeffs" and creating a sequence of values for x but then I don't know how to equalize them.

Thank you very much in advance for your time,

Elena

Hi everyone.

I'm trying to solve the following PDE

but I'm getting this error:

Error, (in pdsolve/numeric/process_IBCs) initial/boundary conditions can only contain derivatives which are normal to the boundary, got (D[1, 1](w))(x, 0)

The PDE represents the bending of a thin plate.

See File:

PDE-Problem.mw

into the "Ask a Question" window?

Nothing to add

Actually I want to ask something else.

(Maple 2015)

For the simple ODE with initial condition
dsolve({ diff(y(t),t) = y(t)^2 - y(t)^3, y(0) = 1/10 }, y(t));

dsolve produces two different answers, almost randomly (even after restart or after closing Maple and reloading the worksheet). Namely:

(1)

(2)

or

(2')

but this simplifies to (2), so it's not a "true" bug.

Notice that (2) is correct but (1) is incorrect even for t=0 (the initial condition!):

evalf(eval(RootOf(-ln(_Z)*_Z+ln(_Z-1)*_Z-ln(10)*_Z-ln(9/10)*_Z-I*Pi*_Z+_Z*t-10*_Z+1),t=0))=1/10;

Maple seems to prefer the wrong solution (1) but occasionally produces (2) e.g. in a new whorksheet!
In earlier versions it seems that only (1)  appears.

The same ODE with another IC

dsolve({ diff(y(t),t) = y(t)^2 - y(t)^3, y(0) = 1/100 }, y(t));

evalf(eval(%,t=0));

is always incorrect. It should be

but Lambert's function never shows up!
Let me mention that only the exact solutions are affected, numeric is ok.

Without an initial condition, dsolve always uses LambertW:

dsolve({ diff(y(t),t) = y(t)^2 - y(t)^3}, y(t));

Can you explain this behavior?

Dear Friends

My problem is related to exporting Maple code to Latex, to be more specific I want Latex output that look exactly same as it appears in Maple worksheet like like Red Maple prompt etc.

I have found another solution to this:

"In file menu go to “Export As”, then save file in Tex format. But in order to run this Tex file you need to have files like maple2e.sty, mapleenv.sty, maplestd2e.sty, mapletab.sty, mapleenv.def, mapleplots.sty, maplestyle.sty, mapleutil.sty(all part of package maplestd2e), these all file can be copied from C:drive where maple is installed, just copy all these file to Latex folder and you can run Tex file exported from Maple."

But this do not give exactly what I want.

Please see Maple file and its exported Tex version(Not possible to upload) along with pdf output for Maple worksheet.

Maple_Export_to_Latex.pdfMy_Paper_on_Zakharov-Kuznetsove_equation_(Simplified_determining_equations_for_classical_symmetries).mw

Let A and B be two lists of monomials. I want a new list C contained that monomials of A where not divide by any monomial of B. For example if A=[x,y,x^2y,xy^2,y^2] and B=[x^2,y^3] then C=[x,y,xy^2,y^2].

I am getting the following expression when I partially differentiate an expression:

PDE11 := diff(theta(z, p), z, z, p)+2*lambda(p)*theta(z, p)*(diff(lambda(p), p))+lambda(p)^2*(diff(theta(z, p), p))+lambda(p)^2*(sin(theta(z, p))-theta(z, p))+2*p*lambda(p)*(sin(theta(z, p))-theta(z, p))*(diff(lambda(p), p))+p*lambda(p)^2*(cos(theta(z, p))*(diff(theta(z, p), p))-(diff(theta(z, p), p)))

I differentiate the above equation to get each term in the form of :table([f=......])

(table([f = 1+sum(Lambda[n](0)/factorial(n), n = 1 .. infinity)]))(p)^2

It is difficult to understand the expression. Maple does not show any error. Can you please tell me what the error is?

Hi, I am trying to find adomian's polynomial of exp(y), but after execution it shows DF(e^y), D^2F(e^y) as it is. why it does'nt show its derivative?plz help meNewADP.2.mw