Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

H := 1/(2*m)*diff(f(q,P,t), q)^2 = -diff(f(q,P,t),t);
assume f = f1(q) + f2(t);
H := 1/(2*m)*diff(f1(q), q)^2 = -diff(f2(t),t);
dsolve(H);

set each side of equation be a constant E which is anticipated integration constant
how to solve to S = q*sqrt(2*m*E) - E*t

and

p = diff(S,q) = sqrt(2*m*E);
Q = diff(S,E) = q*m/sqrt(2*m*E) - t;

There is a desire to explore the process of filtration combustion. To do this, you must solve a system of differential equations in partial derivatives.
I write down all the equations.
Boundary conditions in Maple 2015.0 interpreted incorrectly.
I need to write like that:

 

 

given that:


It turns out so:

 

 

ie somewhere lost derivatives

,....

 

 

As in Maple record boundary conditions correct?

Thanking you in advance.

can maple 18 install in mobile phone Lumia 950?

Hello everyone. 

I can not calculate it here


restart

with(VectorCalculus):

with(plots):

v1 := x^2+y+`-`(sin(z))

x^2+y-sin(z)

(1)

v2 := VectorCalculus:-`+`(VectorCalculus:-`+`(x^2, 1/y), VectorCalculus:-`-`(VectorCalculus:-`*`(2, z)))

x^2+1/y-2*z

(2)

v3 := VectorCalculus:-`+`(VectorCalculus:-`+`(VectorCalculus:-`*`(3, x), y^2), z)

y^2+3*x+z

(3)

vv := VectorField(`<,>`(v1, v2, v3), 'cartesian'[x, y, z])

Vector[column](%id = 18446744074235504150)

(4)

fieldplot3d(vv, x = -1 .. 1, y = -1 .. 1, z = -1 .. 1, arrows = `3-D`, grid = [5, 5, 5], axes = normal)

 

Flux(vv, Sphere(`<,>`(0, 0, 0), r))

Error, (in sprintf) Maple was unable to allocate enough memory to complete this computation.  Please see ?alloc

 

 

``


Probably, the problem with the calculation of surface integrals. I would be grateful for any help.

Download 1.mw

 

Hi all,

I want to know how we can have the result which is made in Maple in Latex in a good style... I must prepare a report of it today, but when I copy the result from Maple to Latex, there're a lot of things to edit, also the results that are obtained from Maple are a lot, so it's very difficult for me to type it in Latex...

Please help me

I have the solution to a heat PDE, v and the error esitmate u + cos(x+t) = v

 

I want to plot log v(1,t) as function of log u(1,t) in maple, but I seem to get an error:

Error, (in plot) unexpected option: ln(u(1, t))

 

I am attaching my code below.

How to fix this problem?

Thanks in advance.

PDE+cos.mw

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.

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]f11+a[1,2]f12 that is, the values of a[j,k] should appear beside the functions' names, like

7f11+2f12-3f21 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?

 

restart;

eq := I*exp(-(2*I)*k*t)*k*sin(theta)*r^2*cos(theta)^3+4*exp(-(2*I)*k*t)*r*cos(theta)^3+2*(diff(Vr(t, r, theta), theta, theta))*cos(theta)*exp(-I*k*(sin(theta)*r+t))-6*(diff(Vr(t, r, theta), theta))*sin(theta)*exp(-I*k*(sin(theta)*r+t))-4*Vr(t, r, theta)*cos(theta)*exp(-I*k*(sin(theta)*r+t))-4*exp(-(2*I)*k*t)*r*cos(theta);

I*exp(-(2*I)*k*t)*k*sin(theta)*r^2*cos(theta)^3+4*exp(-(2*I)*k*t)*r*cos(theta)^3+2*(diff(diff(Vr(t, r, theta), theta), theta))*cos(theta)*exp(-I*k*(sin(theta)*r+t))-6*(diff(Vr(t, r, theta), theta))*sin(theta)*exp(-I*k*(sin(theta)*r+t))-4*Vr(t, r, theta)*cos(theta)*exp(-I*k*(sin(theta)*r+t))-4*exp(-(2*I)*k*t)*r*cos(theta)

(1)

sol := pdsolve(eq);

Vt(t, r, theta) = _F2(t, r)/cos(theta)^2+sin(theta)*_F1(t, r)/cos(theta)^2-((1/2)*I)*(cos(theta)^2*k^2*r^2-2)*exp(I*(sin(theta)*r-t)*k)/(k^3*r^2*cos(theta)^2)

(2)

pdetest(sol, eq);

0

(3)

eq2 := eval(eq, Vr(t,r,theta) = rhs(sol)):
eq2 := simplify(%);

-((1/2)*I)*exp(-(2*I)*k*t)*k*r^2*cos(theta)^3+2*exp(-(2*I)*k*t)*r*sin(theta)*cos(theta)-3*(diff(Vt(t, r, theta), theta))*sin(theta)*exp(-I*k*(sin(theta)*r+t))-2*Vt(t, r, theta)*cos(theta)*exp(-I*k*(sin(theta)*r+t))+(diff(diff(Vt(t, r, theta), theta), theta))*cos(theta)*exp(-I*k*(sin(theta)*r+t))

(4)

evalb(eq2 = 0);

false

(5)

 

 

 

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?

 

 

First 1227 1228 1229 1230 1231 1232 1233 Last Page 1229 of 2254