Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

how to plot exp(f)*erf(g) where f and g are complex function

(-6.328281880*10^(-49)-1.071713312*10^(-47)*I)*exp(-5996.664400+2.000000000*10^14*del+(1.547600000*10^16*I)*del)*(erf(11.+77.38000000*I+1.000000000*10^14*del)-1.*erf(-4.+77.38000000*I+1.000000000*10^14*del))

del from -100e-15 to 100e-15

In this paper we will demonstrate the many differences of implementation in the modeling of mechanical systems using embedded components through Maplesoft. The mechanical systems are used for different tasks and therefore have different structure in its design; as to the nature of the used functional elements placed on them, they vary greatly. This diversity is reflected in approaches and practices in modeling.

The following cases focus on mechanical components of the units manufacturing and processing machines. We can generate graphs for analysis using different dynamic pair ametros; all in real-time considerations in its manufacturing costs from the equations of conservation of energy.
Therefore modeling with Maplesoft ensures the smooth optimum performance in mechanical systems, highlighting the sustainability criteria for other areas of engineering.

 

XXXIII_Coloquio_SMP_2015.pdf

XXXIII_Coloquio_UNASAM_2015.mw

(in spanish)

L.AraujoC.

 

 

I am creating a plot but the numbers on the horizontal axis overlap.  Is there a way to change the directions of the numbers displayed under the horizontal axis (tickmark labels) from horizontal to vertical?

The picture below shows how the numbers are bunched together under the horizontal axis.  I'd like them to be displayed vertical

 

 

Thank you

1. Take a group for example. First we can set up a group G by

G:=<<a,b>|<a2=1,b3=1,(a.b)2=1>>. Actually,G=S3. So how to simplify a long production,

such as "a.b.b.a.b.b.a.b"?

2. How to define a finitely presented algebra over some field, such as the enveloping algebra

or quantum group of a Lie algebra? And moreover how to do the similar computation about 

simplifying a long production?

Hi

Anyone could help me in solving the following system of equations to get constants C1, C2, C3 and C4. MALPE give me this "soution may have been lost".  The MAPLE sheet is also attached.

 

restart:

Eq1:=simplify(C3*exp(-(1/4)*(C2*(x^2-2*0)+sqrt(C2*(x^2-2*0)^2+4*M*(x^2-2*0)*w1*(x^2-2*0)))/w1)+C4*exp((1/4)*(-C2*(x^2-2*0)+sqrt(C2*(x^2-2*0)^2+4*M*(x^2-2*0)*w1*(x^2-2*0)))/w1)-U) = 0;

C3*exp(-(1/4)*(C2*x^2+(x^4*(4*M*w1+C2))^(1/2))/w1)+C4*exp(-(1/4)*(C2*x^2-(x^4*(4*M*w1+C2))^(1/2))/w1)-U = 0

(1)

Eq2:=simplify(exp(-(1/4)*(C2+sqrt(C2^2+4*M*w1))*(x^2-2*0)/w1)*C3*x+exp((1/4)*(-C2+sqrt(C2^2+4*M*w1))*(x^2-2*0)/w1)*C4*x+C2-V-z) = 0;

exp(-(1/4)*(C2+(C2^2+4*M*w1)^(1/2))*x^2/w1)*C3*x+exp(-(1/4)*(C2-(C2^2+4*M*w1)^(1/2))*x^2/w1)*C4*x+C2-V-z = 0

(2)

Eq3:=simplify((-2*w2*w5*ln(C3*exp(-(1/2)*sqrt(w2*w4*(w2*w4+w3*w6))*C2*(x^2-2*0)/(w2*w4*w5))-C4)*sqrt(w2*w4*(w2*w4+w3*w6))+w2*w5*(-w2*w4+sqrt(w2*w4*(w2*w4+w3*w6)))*ln(exp(-(1/2)*sqrt(w2*w4*(w2*w4+w3*w6))*C2*(x^2-2*0)/(w2*w4*w5)))+C1*w3*w6*sqrt(w2*w4*(w2*w4+w3*w6)))/(sqrt(w2*w4*(w2*w4+w3*w6))*w3*w6)-1)= 0;

(-ln(exp(-(1/2)*(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5)))*w2^2*w4*w5+C1*w3*w6*(w2*w4*(w2*w4+w3*w6))^(1/2)-2*w2*w5*ln(C3*exp(-(1/2)*(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5))-C4)*(w2*w4*(w2*w4+w3*w6))^(1/2)+ln(exp(-(1/2)*(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5)))*(w2*w4*(w2*w4+w3*w6))^(1/2)*w2*w5-(w2*w4*(w2*w4+w3*w6))^(1/2)*w3*w6)/((w2*w4*(w2*w4+w3*w6))^(1/2)*w3*w6) = 0

(3)

Eq4:= simplify((-C2*x^2*w2*w4-.50*C2*x^2*w3*w6+sqrt(w2*w4*(w2*w4+w3*w6))*C2*x^2+2.*w2*w4*w5*ln(w3^4*w6^2*(C3^2*exp(-1.0*sqrt(w2*w4*(w2*w4+w3*w6))*C2*x^2/(w2*w4*w5))-2*C3*exp(-.5*sqrt(w2^2*w4^2+w2*w3*w4*w6)*C2*x^2/(w2*w4*w5))*C4+C4^2)/(w2*w4*(w2*w4+w3*w6)*C2^2))-5.544000000*w2*w4*w5-w3^2*w6)/(w3^2*w6)) = 0;

(-C2*x^2*w2*w4-.5000000000*C2*x^2*w3*w6+(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2+2.*w2*w4*w5*ln(w3^4*w6^2*(C3^2*exp(-(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5))-2.*C3*exp(-.5*(w2*w4*(w2*w4+w3*w6))^(1/2)*C2*x^2/(w2*w4*w5))*C4+C4^2)/(w2*w4*(w2*w4+w3*w6)*C2^2))-5.544000000*w2*w4*w5-w3^2*w6)/(w3^2*w6) = 0

(4)

solve({Eq1, Eq2, Eq3,Eq4}, {C1, C2, C3,C4});

Warning, solutions may have been lost

 

``

``

Download solution_lost.mw

 

 

Hi all,

I have this system

> system1D := H = alpha*gamma[2, 2]*d[2, 1]-beta*d[1, 2]*gamma[1, 2]^2-gamma*d[1, 2]*gamma[2, 1]^2+alpha*gamma[2, 2]^2*d[2, 2]-beta*d[2, 2]*gamma[2, 2]^2-gamma*d[2, 2]*gamma[2, 2]^2, E = alpha*gamma[2, 1]*d[1, 1]-beta*d[1, 2]*gamma[1, 1]-gamma*d[1, 1]*gamma[2, 1]+alpha*gamma[2, 1]^2*d[1, 2]-beta*d[2, 2]*gamma[2, 1]-gamma*d[2, 1]*gamma[2, 2], B = alpha*gamma[1, 1]*d[2, 1]-beta*d[1, 1]*gamma[1, 1]^2-gamma*d[1, 1]*gamma[1, 1]^2+alpha*gamma[1, 1]^2*d[2, 2]-beta*d[2, 1]*gamma[2, 1]^2-gamma*d[2, 1]*gamma[1, 2]^2, D = alpha*gamma[1, 2]*d[2, 1]-beta*d[1, 1]*gamma[1, 2]^2-gamma*d[1, 2]*gamma[1, 1]^2+alpha*gamma[1, 2]^2*d[2, 2]-beta*d[2, 1]*gamma[2, 2]^2-gamma*d[2, 2]*gamma[1, 2]^2, A = alpha*gamma[1, 1]*d[1, 1]-beta*d[1, 1]*gamma[1, 1]-gamma*d[1, 1]*gamma[1, 1]+alpha*gamma[1, 1]^2*d[1, 2]-beta*d[2, 1]*gamma[2, 1]-gamma*d[2, 1]*gamma[1, 2], C = alpha*gamma[1, 2]*d[1, 1]-beta*d[1, 1]*gamma[1, 2]-gamma*d[1, 2]*gamma[1, 1]+alpha*gamma[1, 2]^2*d[1, 2]-beta*d[2, 1]*gamma[2, 2]-gamma*d[2, 2]*gamma[1, 2], F = alpha*gamma[2, 1]*d[2, 1]-beta*d[1, 2]*gamma[1, 1]^2-gamma*d[1, 1]*gamma[2, 1]^2+alpha*gamma[2, 1]^2*d[2, 2]-beta*d[2, 2]*gamma[2, 1]^2-gamma*d[2, 1]*gamma[2, 2]^2, G = alpha*gamma[2, 2]*d[1, 1]-beta*d[1, 2]*gamma[1, 2]-gamma*d[1, 2]*gamma[2, 1]+alpha*gamma[2, 2]^2*d[1, 2]-beta*d[2, 2]*gamma[2, 2]-gamma*d[2, 2]*gamma[2, 2], H = alpha*delta[2, 2]*d[2, 1]-beta*d[1, 2]*delta[1, 2]^2-gamma*d[1, 2]*delta[2, 1]^2+alpha*delta[2, 2]^2*d[2, 2]-beta*d[2, 2]*delta[2, 2]^2-gamma*d[2, 2]*delta[2, 2]^2, E = alpha*delta[2, 1]*d[1, 1]-beta*d[1, 2]*delta[1, 1]-gamma*d[1, 1]*delta[2, 1]+alpha*delta[2, 1]^2*d[1, 2]-beta*d[2, 2]*delta[2, 1]-gamma*d[2, 1]*delta[2, 2], B = alpha*delta[1, 1]*d[2, 1]-beta*d[1, 1]*delta[1, 1]^2-gamma*d[1, 1]*delta[1, 1]^2+alpha*delta[1, 1]^2*d[2, 2]-beta*d[2, 1]*delta[2, 1]^2-gamma*d[2, 1]*delta[1, 2]^2, D = alpha*delta[1, 2]*d[2, 1]-beta*d[1, 1]*delta[1, 2]^2-gamma*d[1, 2]*delta[1, 1]^2+alpha*delta[1, 2]^2*d[2, 2]-beta*d[2, 1]*delta[2, 2]^2-gamma*d[2, 2]*delta[1, 2]^2, A = alpha*delta[1, 1]*d[1, 1]-beta*d[1, 1]*delta[1, 1]-gamma*d[1, 1]*delta[1, 1]+alpha*delta[1, 1]^2*d[1, 2]-beta*d[2, 1]*delta[2, 1]-gamma*d[2, 1]*delta[1, 2], C = alpha*delta[1, 2]*d[1, 1]-beta*d[1, 1]*delta[1, 2]-gamma*d[1, 2]*delta[1, 1]+alpha*delta[1, 2]^2*d[1, 2]-beta*d[2, 1]*delta[2, 2]-gamma*d[2, 2]*delta[1, 2], F = alpha*delta[2, 1]*d[2, 1]-beta*d[1, 2]*delta[1, 1]^2-gamma*d[1, 1]*delta[2, 1]^2+alpha*delta[2, 1]^2*d[2, 2]-beta*d[2, 2]*delta[2, 1]^2-gamma*d[2, 1]*delta[2, 2]^2, G = alpha*delta[2, 2]*d[1, 1]-beta*d[1, 2]*delta[1, 2]-gamma*d[1, 2]*delta[2, 1]+alpha*delta[2, 2]^2*d[1, 2]-beta*d[2, 2]*delta[2, 2]-gamma*d[2, 2]*delta[2, 2];


> subs({A = 0, B = 0, C = 0, D = 0, E = 0, F = 0, G = 0, H = 0, delta[1, 1] = 1, delta[1, 2] = 0, delta[2, 1] = 0, delta[2, 2] = 0, gamma[1, 1] = 1, gamma[1, 2] = 0, gamma[2, 1] = 0, gamma[2, 2] = 0, delta[1, 1]^2 = 0, delta[1, 2]^2 = 0, delta[2, 1]^2 = 1, delta[2, 2]^2 = 0, gamma[1, 1]^2 = 0, gamma[1, 2]^2 = 1, gamma[2, 1]^2 = 0, gamma[2, 2]^2 = 0}, {system1D});

The problem is: there is any simple way to use command "subs" when some expression such that delta[1,1]=1, gamma[1,1]=1, gamma[1,2]^2=1 have value and others are zero.

Can someone please advice and help me on this?

thanks

witribm

Hello to you all,

I have this DE

but when I try to change the variable

I get this

>algsubs(r(t) = 1/u(t), diff(diff(r(t), t), t)-k/r(t) = -K/r(t)^2);

Why is it not done for the rest of the terms.  Is there a more easy way to do it.  Take into account that I cannot use dsolve because it is in an course in integral.

 

Thank in advance for your help.

 

--------------------------------------
Mario Lemelin
Maple 2015 Ubuntu 14.04 - 64 bits
Maple 2015 Win 10 Pro - 64 bits messagerie : mario.lemelin@cgocable.ca téléphone :  (819) 376-0987

Does anyone know an easy way to convert %d_ to D? I tried the convert command, but no effect.  Thanks.

 

Michael

Hi!

For my lecture advanced dynamics, I got the question to make a sketch of the coordinate system and the particle’s motion, with the following position vector:

r(t) = (a cos(ωt))i + (b sin(ωt))j + (ct)k

In this function, ω is in [rad/s], a and b are in [m], t in [s] and c is [m/s].

I have no idea what this should look like, and I was thinking about plotting this with Maple, but thus far I'm not succeeding in plotting it. Which function or method can I use to visualize this in Maple?

Thanks a lot!

 

Hi all,

I have this system

> system2hD := 0 = -d[1, 2], 0 = -d[2, 1], 0 = -delta*d[1, 2], 0 = -d[1, 1]-delta*d[1, 1], 0 = -delta*d[1, 1]-d[2, 2];
> solve({system2hD, delta <> 0, delta <> 1}, {d[1, 1], d[1, 2], d[2, 1], d[2, 2]});
      {d[1, 1] = 0, d[1, 2] = 0, d[2, 1] = 0, d[2, 2] = 0}

But, I hope to get the answers as {d[1, 1] = -delta*d[1, 1], d[1, 2] = 0, d[2, 1] = 0, d[2, 2] = -delta*d[1, 1]}. I think, I should not use the command "solve".

Can someone please advice and hepl me on this?

thanks

witribm

Apparently inconsistent behaviour of the BesselJ() function.

Examples: BesselJ(-3, 0)  ... gives 0 (correct)

but BesselJ(-3.0, 0), BesselJ(-3, 0.0)  and BesselJ(-3, 0.0) all give Float(infinity) (wrong! - should be 0.0)

The problem seems to occur for all negative integer values of the first argument (the order) when the second argument is 0 or 0.0.

Hello all,

Yesterday, I upgraded to Mac OS X El Capitan.

Now, when working with Maple 2015, I feel the gui is very slow and sometimes irresponsive, when trying to scroll through my worksheets as well as through help pages (e.g., "help plot3d"). When I do the same within Maple 18, it works without any problems.

Does anybody else have the same issue?

Cheers,

Franky

Below is a calculation with 200 Volt and 100 Ω.
The output is 1 Ampere at -45 degrees.
This is useful for electrical circuit calculations, and I can get the right result on a TI-nspire.

However I can't make maple 2015 do it. It's resembling polar coordinates but I can't seem to find any info I can understand.

"((200&angle;0&deg;))/((100&angle;45&deg;)+(100&angle;45&deg;))=\

1&angle;-45&deg;"

below is the same equation but in a picture so it can be read by humans:

http://imgur.com/3i2Su0B

 

Hopefully theres some simple comands that can help make this easy.

Thanks

For example, I have a commutation table of finite Lie algebra. I'd like to extract some information from it (mainly, I'm interested in maximal dimension of abelian Lie subalgebra).

if eval(simplify( subs(a=2,subs(b= 5,f))))) - n = 0 then

      hello := union(hello, {f})

end if:

after subs(a=2,subs(b= 5,f))) to verify , some of hello return -n

negative

why it is possible to pass the if condition statement?

if it is normal case, how to write correct if statement to ensure no -n 

First 1232 1233 1234 1235 1236 1237 1238 Last Page 1234 of 2255