MaplePrimes Questions

So I'm new to maple, more used to working in matlab. Anyway, I'm trying to solve the pendulum equation with a friction term, but the "dsolve" command spits out something I have no clue how to handle. What do I do with the "theta = _a & where [loads of bullshit]"? It's supposed to give me an integral of some sort that I can work with, not a strange substitution. Halp! 

 

I got bit by Maple when I tried to combine a number of 2-d plots into one 3-d plot, using plottools:-transform to space the plots out in the third dimentsion. I have done such things before. This time, the first time I executed the final plot statement the GUI got clobbered and somehow also the whole worksheet (I may have saved it). Now, when I load it nothing shows up and worse, I cannot enter anything in the worksheet, it is completely screwy. I looked at the XML code in Emacs and it validates and I cannot find anything obviously wrong with it, i.e. no dangling tags or similar.

I used Maple 2017. In Maple 2015 the same worksheet overprints itself & makes a big mess on the screen. It is similarly unuseable.

Can one of the experts have a look at this? Although rare I have experienced similar issues on occasion before & would like to know how to get this into a sane state.

TIA,

Mac Dude

triangle_analysis.mw

 

I need to find all integer solutions of these multivariable equations (and many more similar equations)

I have the boundary conditions that b-1 ≥ a1 > a2 ≥ 0

How can I do this in Maple?

For context, b is a base and a1 and a2 are digits in a b-adic number where a1 is bigger than a2. Solutions to these equations correspond to numbers that have a particular property.

 

I have posted this question on math.stackexchange, but saught an algebraic approach rather than a computational one. Here it is shown that a1 and a2 are linear functions of b. This is where the families of solutions come in. Mathloves answer is detailed and potentially useful in regards to formatting the answer.

https://math.stackexchange.com/questions/3001095/integer-solutions-of-a-variable-coefficient-polynomial

 

Thank you for any help / advice. Ben

I simply used solve to solve an equation but it doesn't work properly , Here is my code :

I want to solve system of nonlinear ODEs using RKF with shooting technique. If some body have a idea please help me coding I am new In Maple. If somebody have a relevant code please share with me. I try to understand and then apply on my problem. 

I want to solve attached link paper and image.Thanks

Hi,

I am trying to compute the following integral:

int(a*b*c/sqrt((x-a)^2 + (y-b)^2 + (z-c)^2), [a=-1..1, b=-1..1, c=-1..1])

where x, y, z are Real parameters.

The integral looks so simple but Maple can't compute it as is.

Thanks

Hi

I would like to solve the integrodifferential equation and then look to the  stability of the origin.

Is it  stable, uniformly stable, asymptotically stable and uniformly asumptotically stable.

Please see the following code.

Code.mw

Thanks

 

 

I have to use the Chinese Remainder Theorem but once I have found m^e I cant seem to figure out a way that works to get just the message without the exponent. I assumed that the exponent is 3 since the message was encrypted 3 different times. How can I do this with Maple?

CRA_project.mw

Hi all. how can i detect and remove the lists with same members? for example in the code below, i want to remove duplicate of the lists with the same members, for example [1,1,5],[1,5,1],[5,1,1]. i want just one of these sets to remain and the others to be removed. it doesn't matter which. i want just one set with members 1,5,1 which their summation is 7. and the same thing for other lists. thanks in advance.
 

restart:with(combinat):with(Iterator):

MM:=composition(7, 3)

{[1, 1, 5], [1, 2, 4], [1, 3, 3], [1, 4, 2], [1, 5, 1], [2, 1, 4], [2, 2, 3], [2, 3, 2], [2, 4, 1], [3, 1, 3], [3, 2, 2], [3, 3, 1], [4, 1, 2], [4, 2, 1], [5, 1, 1]}

(1)

 


 

Download code.mw

I know this is a trivial worksheet (in document mode) but why does the second subs fail? Thanks  Gravity.mw
 

``

Justification for using mg as the force exerted by mass m in gravity field g

Using Newton's Law of Gravity, find ``

work := `assuming`([int(G*M*m/r^2, r = R .. R+h)], [0 < R, 0 < h])

G*M*m*h/(R*(R+h))

(1)

``

R is >> h, so let R + h = R

work := subs(R+h = R, work)

G*M*m*h/R^2

(2)

``

G, M and R are all constants, so let GM/R² = g

work := subs(G*M/R^2 = g, work)

G*M*m*h/R^2

(3)

``

``


 

Download Gravity.mw

 

I am having problems when trying to sue de colorscheme opiton. I always get an error message. The example below is from the Maple homepage:

plot3d(sin(x)*cos(y), x = 0 .. 4*Pi, y = 0 .. 4*Pi, view = [default, default, -3 .. 3], colorscheme = ["xgradient", ["Green", "Purple", "Blue"]]);
Error, (in plot3d) unexpected option: colorscheme = ["xgradient", ["Green", "Purple", "Blue"]]


Any hint?

Hi I'm done this and I really have no Idea what's wrong with maple :

 

I have no idea how to solve the equation ...
Thanks in advance !

Hello, i am trying to solve the PDE

 (E*I) d4/dx4 [w(x,t)] + (rho*A) d2/dt[w(x,t)]=0

where w(x,t) = W(x) * cos (ω *t)

W(x) is the amplitude of the beam, the PDE is solved for W(x) by seperation of variables.

my first question is how to simplify the equations by replacing the e^ terms with alpha as done in the lecture, so that things look cleaner. And the next question is how do i proceed to input the Boundary conditions to solve the problem? I have looked in quite a lot of places but there is no example of a second or third order boundary condition(if that makes sense)

 beam_vibration_displacement.mwMech_vibrations.pdf

Hi,

I would like to ask you a question on the following program. Where does the error come from?

Thank you for your help in advance.

 

som:=0:

for b1 from 10 to 20 by 1 do
for b2 from 1 to 10 by 0.1 do
for alpha from 0.5 to 0.9 by 0.1 do
for beta from 0.1 to 0.4 by 0.1 do
for c from 1 to 1 by 1 do
for f from 1 to 10 by 1 do
for g from 8 to 470 by 1 do
for lambdai from 0.2 to 0.2 by 0.1 do
for lambdaj from 0.2 to 0.2 by 0.1 do
for gammai from 0.4 to 0.4 by 0.1 do
for gammaj from 0.4 to 0.4 by 0.1 do

aiSQ:=(alpha*b1)/(alpha*b2+beta*b2+c);
ajSQ:=(beta*b1)/(alpha*b2+beta*b2+c);
UiSQ:=(1/2)*alpha*b1^2*(alpha^2*b2+2*alpha*beta*b2+c*alpha+beta^2*b2+2*beta*c)/(alpha*b2+beta*b2+c)^2;
UjSQ:=(1/2)*beta*b1^2*(alpha^2*b2+2*alpha*beta*b2+2*c*alpha+beta^2*b2+beta*c)/(alpha*b2+beta*b2+c)^2;
USQ:=(1/2)*b1^2*(alpha+beta)*(alpha*b2+beta*b2+2*c)/(alpha*b2+beta*b2+c)^2;

ai:=(c*b1*alpha+c*p*f+p*f*beta*b2-alpha*b2*p*f)/(c*(alpha*b2+beta*b2+c));
aj:=(-p*f*beta*b2+alpha*b2*p*f+c*beta*b1+c*p*f)/(c*(alpha*b2+beta*b2+c));
aineg:=(-p*f*lambdai*b2*beta+c*b1-c*b1*lambdai+b2*p*f*lambdaj+b2*alpha*p*f*lambdaj-p*f*lambdai*b2+c*b1*alpha-c*p*f*lambdai)/(c*(b2*alpha-b2*lambdai+b2*beta-b2*lambdaj+c+2*b2));
ajneg:=-(b2*alpha*p*f*lambdaj-p*f*lambdai*b2-p*f*lambdai*b2*beta+b2*p*f*lambdaj+c*p*f*lambdaj-c*b1-c*beta*b1+c*lambdaj*b1)/(c*(b2*alpha-b2*lambdai+b2*beta-b2*lambdaj+c+2*b2));
ui:=alpha*(b1*(aineg+ajneg)-(b2/2)*(aineg+ajneg)^2)-(c/2)*aineg^2-p*f*(aineg-ai);
uj:=beta*(b1*(aineg+ajneg)-(b2/2)*(aineg+ajneg)^2)-(c/2)*ajneg^2-p*f*(ajneg-aj);
u:=(b1*(aineg+ajneg)-(b2/2)*(aineg+ajneg)^2)+p*f*(aineg-ai)+p*f*(ajneg-aj);
eqti:=gammai*(u-tj-USQ)-((1-gammai)/(1-lambdai))*(ui-UiSQ);
eqtj:=gammaj*(u-ti-USQ)-((1-gammaj)/(1-lambdaj))*(uj-UjSQ);
solt:=solve({eqti, eqtj}, {ti, tj});
ti:=subs(solt,ti);
tj:=subs(solt,tj);
dai:=diff(ai,p);
daj:=diff(aj,p);
daineg:=diff(aineg,p);
dajneg:=diff(ajneg,p);
dti:=diff(ti,p);
dtj:=diff(tj,p);
eqp:=(b1-b2*(ai+aj))*(dai+daj)-dti-dtj+f*(aineg-ai)+p*f(daineg-dai)+f*(ajneg-aj)+p*f(dajneg-daj)-g*p=0;
p:=solve(eqp,p);

Uip:=alpha*(b1*(ai+aj)-(b2/2)*(ai+aj)^2)-(c/2)*ai^2-p*f*(aineg-ai)+(1-lambdai)*ti;
Ujp:=beta*(b1*(ai+aj)-(b2/2)*(ai+aj)^2)-(c/2)*aj^2-p*f*(ajneg-aj)+(1-lambdaj)*tj;
Up:=(b1*(ai+aj)-(b2/2)*(ai+aj)^2)+p*f*(aineg-ai)+p*f*(ajneg-aj)-((g^2)/2)*p-ti-tj;

CSQ:=b1-b2*(aiSQ+ajSQ);
Cabat:=b1-b2*(ai+aj);
Cneg:=b1-b2*(aineg+ajneg);

if (CSQ>0 and Cabat>0  and Cneg>0 and beta<alpha and p>0 and p<1)
then
print(b1,b2,alpha,beta,c,f,g,lambdai,lambdaj,`aiSQ=`,aiSQ,`ajSQ=`,ajSQ,`UiSQ=`,UiSQ,`UjSQ=`,UjSQ,`USQ=`,USQ,`ai=`,ai,`aj=`,aj,`aineg=`,aineg,`ajneg=`,ajneg,`ti=`,ti,`tj=`,tj,`p=`,p,`Uip=`,Uip,`Ujp=`,Ujp,`Up=`,Up);
som:=som+1;
fi;
od;od;od;od;od;od;od;od;od;od;od;
som;

10, 1, 0.5, 0.1, 1, 1, 8, 0.2, 0.2, aiSQ=, 3.12500000000000, ajSQ=, 

  0.625000000000000, UiSQ=, 10.3515625000000, UjSQ=, 2.85156250000000, USQ=, 

  30.4687500000000, ai=, 3.14314584993738, aj=, 0.667340316520551, aineg=, 

  4.06068541500626, ajneg=, 2.80826596834794, ti=, 19.2559726892491, tj=, 

  7.44890809994178, p=, 0.0483889331663441, Uip=, 25.8431760998162, Ujp=, 

  8.71735374786572, Up=, 2.73962847065372
Error, (in solve) a constant is invalid as a variable, 7.44890809994178, 19.2559726892491
                                      1

is possible to solve this pde via maple?

m1.mwm1.mw
 

restart

sys := [-(-r^2+1)*(diff(theta(r, z), z))+(diff(theta(r, z), r)+r*(diff(theta(r, z), r, r)))/r+diff(theta(r, z), z, z)+(diff(theta(r, z), r))*(diff(sigma(r, z), r))+(diff(sigma(r, z), z))*(diff(theta(r, z), z))+(diff(theta(r, z), r))^2+(diff(theta(r, z), z))^2 = 0, -(-r^2+1)*(diff(sigma(r, z), z))+(diff(sigma(r, z), r)+r*(diff(sigma(r, z), r, r)))/r+diff(sigma(r, z), z, z)+(diff(theta(r, z), r)+r*(diff(theta(r, z), r, r)))/r+diff(theta(r, z), z, z) = 0]; IBCs := {sigma(1, z) = 1, sigma(r, 0) = 1, theta(1, z) = 1, theta(r, 0) = 1, (D[1](sigma))(0, z) = 0, (D[1](theta))(0, z) = 0, (D[2](sigma))(r, 1) = 0, (D[2](theta))(r, 1) = 0}

[-(-r^2+1)*(diff(theta(r, z), z))+(diff(theta(r, z), r)+r*(diff(diff(theta(r, z), r), r)))/r+diff(diff(theta(r, z), z), z)+(diff(theta(r, z), r))*(diff(sigma(r, z), r))+(diff(sigma(r, z), z))*(diff(theta(r, z), z))+(diff(theta(r, z), r))^2+(diff(theta(r, z), z))^2 = 0, -(-r^2+1)*(diff(sigma(r, z), z))+(diff(sigma(r, z), r)+r*(diff(diff(sigma(r, z), r), r)))/r+diff(diff(sigma(r, z), z), z)+(diff(theta(r, z), r)+r*(diff(diff(theta(r, z), r), r)))/r+diff(diff(theta(r, z), z), z) = 0]

 

{sigma(1, z) = 1, sigma(r, 0) = 1, theta(1, z) = 1, theta(r, 0) = 1, (D[1](sigma))(0, z) = 0, (D[1](theta))(0, z) = 0, (D[2](sigma))(r, 1) = 0, (D[2](theta))(r, 1) = 0}

(1)

NULL


 

Download m1.mw

 

4 5 6 7 8 9 10 Last Page 6 of 1688