MaplePrimes Questions

I’m using Maple 2015.2 to plot irregular spaced data using the surfplot command. The shape of the plot is as expected, however I’ve been unsuccessful in changing the color of the surface. I would like the surface to be a consistant color instead of the “rainbow” color scheme. I’ve tried both the color and colorscheme options to no avail.  I’m sure I’m making a very basic mistake, however I’m having difficulty finding it.

 

See attached worksheet

 feildstrength.mw

Best regards,

Ron

[EDIT]: Worksheet didn't attach the first time

why is simplify is making errors - me or maple?

Ec := (1/2)*Kc*(2*Pi*h0/Pi(a0^2+h0^2)-2*Pi*h/Pi(a^2+h^2))^2;
2
1 / 2 Pi h0 2 Pi h \
- Kc |------------- - -----------|
2 | / 2 2\ / 2 2\|
\Pi\a0 + h0 / Pi\a + h //
simplify(%);
2
2 Kc (-h0 + h)

 

Simplify is erasing the varibles a0 and a. Is there a secret to using simplify?

In Maple 2015.1 we have

restart;

solve([sin(2*x)/cos(x+3*Pi/2)=1,  x>-4*Pi, x<-5*Pi/2], x, allsolutions, explicit);

solve([sin(2*x)/cos(x+3*Pi/2)=1, x>0, x<2*Pi], x, allsolutions, explicit);

 

 

In the first example, the error message is not clear (actually there exists a unique root  x=-11*Pi/3), in the second example, one root  (x=5*Pi/3) is lost.

 

Hello,

I am trying to write 

D= d/dx + y1*d/dy + y*y1/x*d/dy1 

So that I can apply D to some function L(x,y)

But I don't know how to write the derivative like this without them operation on a function. Can somebody help me? 

Thank you

 

I've got this huge chunk of code which leads to an optimiazation at the very last line (Bestangles:=minimize(maximize()-minimize))). This minization is taking a very long time (havent solved it yet) and I would very much like to reduce that time. As I've understood maple does optimization by differentiating and then finding all extremes and comparing. Would this mean that since I minimize and optimize within a minimization command, it differentiates a ton of times? And if this is the case, can I somehow do the differentiation beforehand, since it is the same function being differentiate all the time? Or is there some other way I can improve the code? 
Thanks alot!

Heres the full code:

 

 

I am very beginner about Maple. How to get the general solution from the follwoing equations by Maple. Please help me. Its Urgent. Please help me out.

why the the software can't plot the function like x^(4/3)*sin(1/x) or x^(1/3)

it could only plot where x>0,but the value is does exist where x<0.

Thanks in advance for your help.

If I have an N-dimensional vector V and a polytope defined as the set of solutions to the equation A*x = b, where A is a d X N matrix, and b is a d X 1 vector, how can I project V onto the surface defined as above?

Thanks!

 

Consider a taper steel plate of uniform thickness t := 25mm as shown in the Fig. In addition to its self weight, the plate is subjected to a point load P := 100N at its mid point. Find the global force vector [F] , global stiffness matrix [K] , displacement in each element (1 and 2) , stresses in each element  (1 and 2) and reaction force at the support.Take E := 2*10^5N/mm2; rho := 8.2*10^(-5)kg/m3;

restart

t__1 := 150:

t__3 := 75:

w := 25:

l := 600:

t__2 := (t__1-t__3)/l*((1/2)*l)+t__3 = 225/2

A__1 := t__1*w = 3750``

A__2 := t__2*w = 5625/2``

A__3 := t__3*w = 1875``

Revised areas:

A__1e := (A__1+A__2)*(1/2) = 13125/4``

A__2e := (A__2+A__3)*(1/2) = 9375/4``

  E := 2*10^11:m2; F__1 := R__1:is support reaction N; F__2 := 100:N;``

rho__1 := 82*10^(-6) = 41/500000  N/mm2

rho__2 := 82*10^(-6) = 41/500000 N/mm2

l := 600:``

Number of elements,

n__e := 2:

l__e := 300 = 300````

q__0 := 100:N/m ; l := 1: m; n__e := 4:  elementsl  l__e := l/n__e: m;

We shall consider a two element system as shown in the Fig.
For element 1 Stiffness matrix K is

                                           Vector[row](2, {(1) = 1, (2) = 2})
K__1 := A__1e*E/l__e.(Matrix(2, 2, {(1, 1) = 1, (1, 2) = -1, (2, 1) = -1, (2, 2) = 1})) = Matrix([[2625000000000000, -2625000000000000], [-2625000000000000, 2625000000000000]])  Vector(2, {(1) = 1, (2) = 2})

For element 2 Stiffness matrix K is

                                         Vector[row](2, {(1) = 2, (2) = 3})
K__2 := A__2e*E/l__e.(Matrix(2, 2, {(1, 1) = 1, (1, 2) = -1, (2, 1) = -1, (2, 2) = 1})) = Matrix([[1875000000000000, -1875000000000000], [-1875000000000000, 1875000000000000]])  Vector(2, {(1) = 2, (2) = 3})

Global stiffness matrix obtained by adding all the elemental stiffness matrices and given b

           Vector[row](3, {(1) = 0, (2) = 0, (3) = 0})

K__g := Matrix(3, 3, {(1, 1) = K__1[1, 1], (1, 2) = K__1[1, 2], (1, 3) = 0, (2, 1) = K__1[2, 1], (2, 2) = K__1[1, 2]+K__2[1, 1], (2, 3) = K__2[1, 2], (3, 1) = 0, (3, 2) = K__2[2, 1], (3, 3) = K__2[2, 2]}) = Matrix([[K__1[1, 1], K__1[1, 2], 0], [K__1[2, 1], K__1[1, 2]+K__2[1, 1], K__2[1, 2]], [0, K__2[2, 1], K__2[2, 2]]])  Vector(3, {(1) = 0, (2) = 0, (3) = 0})

For element 1 Load matrix F is

  F__1e := (1/2)*`&rho;__1`*A__1e*l__e*(Vector(2, {(1) = 1, (2) = 1})) = Vector[column]([[861/25600], [861/25600]]) Vector(2, {(1) = 1, (2) = 2})

``

For element 2 Load matrix F isNULL

F__2e := (1/2)*A__2e*l__e*`&rho;__2`*(Vector(2, {(1) = 1, (2) = 1})) = Vector[column]([[123/5120], [123/5120]]) 

``

 

Download wrong_answers.mwwrong_answers.mwwrong_answers.mw

Ramakrishnan V

rukmini_ramki@hotmail.com

Greetings all,

 

I was wondering if any Maple users out there are using the lastest Microsoft tablet, the Surface Pro 4 with Maple?

I would be interested in any opinions and specifically if the pen is effective when used as an input device with Maple

 

Thanks

hello everyone

can any one tell me what is this anti reduction method. In the paper of serdal palmuk,the link is given bellow

http://www.hindawi.com/journals/mpe/2009/202307/

in this paper question #4 is first solved by anti reduction method for  exact solution.

but i dont understand this method,

if anybody know this then please also tell me how to solve this,

and in the next  (6 & 7 ) examples "in the pourus media equation" they first find its particular exact solution.i also dont understand this,so please tell me

actually i know how to solve ODE to find its exact solution but  i dont know how we find exact solutions of partial differtial equations,

so please help me to solve this problem

thanks

 

I wish to define a function with a finite number of inputs, but I do not know that number ahead of time (in other words the user will specify n and my function operates on vectors of size n). How can this be done?

How to find the integral
,

assuming k and n  integer?
It is known (McCrea W. H., Whipple F. J. W.Random paths in two and three dimensions, Proc. Roy. Soc. Edinburgh. 1940. V. 60. P. 281–298) that

G(n,n)=2/Pi*sum(1/(2*k-1),k=1..n).

The general case is reduced to the case k=n.
This is not a creature of pure reason: the one appears in electric circuits
(see M. Skopenkov, A. Paharev, A. Ustinov, Through resistor net, Mat. pros. Issue 18 (2014), 33-65, in Russian, http://www.mccme.ru/free-books/matpros/pdf/mp-18.pdf).
I found G(8,8) = 182144/(45045*Pi) in 657.797 s and G(9,9) = 3186538/(765765*Pi) in 4157.687 s on my comp by

restart; s := time():(1/2)*VectorCalculus:-int((1-cos(9*Pi*x)*cos(9*Pi*y))/(sin((1/2)*Pi*x)^2+sin((1/2)*Pi*y)^2), [x, y] = Rectangle(0 .. 1, 0 .. 1)); time()-s;
Mathematica 10.3.0 does G(9,9) in 250.391 s on my comp.

 

I have run into a "funny" feature of 2-D input: It seems to convert something like k/2 into this k*`^`(2,-1). While this would often not be an issue (it is correct after all), it becomes a problem when used in an argument list to a procedure. It becomes even more of a problem when, by chance, I have overloaded `^` to act on specific types that I have defined.

Let me try to explain briefly. I have a package called "Lattice" that does whatever it does (not of relevance here). I am writing a little manual for this package, for which I use 2-D input so I can write it in Maple and have the examples right in it and "live".

Here is what happens:

with(Lattice) # load the package

QFh:=Quad(0,kf/2) # Define an element for Lattice

Error, invalid input: Lattice:-`^` expects its 1st argument, element, to be of type Element, but received 2

Copy-pasting kf/2 into a 1-d worksheet, I get

QFh := Quad(0, kf*Lattice[`^`](2, -1));

So it uses Lattice[`^`] which actually appears to bypass the overload I have in the Lattice package. `^` is defined like this in Lattice:

`^`:=proc(element::Element,n::algebraic) option overload; # Element is a defined type in Lattice
...
end proc;

How can I possibly rewrite `^` to fall-back to Maple's ^ operator when called as Lattice[`^`] ?? I know there is a function overload() but have no experience with it. Would it even help?

Or am I missing something completely here? I do not use 2-D input for my usual work, but in this case I want and need to use it. The reason for its bizarre rewrite of "/2" is beyond me. Note that I can replace /2 by *0.5; but that causes problems later on for algebraic work as 1/2 is not 0.5 in Maple. I tried *1/2 but that has the same problem.

Has anyone a clean solution for this? I assume this effect is not limited to my own package but would affect others as well.

M.D.

PS: I ran into this using Maple 15 but I doubt it is specific to this particular version.

Hello,

My code records the values I need, however, I need to implement a modulo of 2*Pi on my result for theta. But this leads to a graph with no plots and I'm not sure how to fix it. Any help is greatly aprreciated! Thank you in advance!

Kind regards,

Gam

with(plots):

a := 1.501*10^9:

Th := sqrt(4*Pi^2*a^3/(G*(Mh+Msat))):

HyperionOrbit := proc (`&theta;IC`, `&omega;IC`, n) local a, Mh, Msat, G, e, beta, M, Eqns, ICs, soln; option remember; global `&omega;H`, Th; a := 1.501*10^9; Mh := 5.5855*10^18; Msat := 5.6832*10^26; G := 6.67259/10^11; e := .232; beta := .89; M := Mh+Msat; Eqns := diff(theta(t), t) = omega(t), diff(omega(t), t) = -G*Msat*beta^2*(xH(t)*sin(theta(t))-yH(t)*cos(theta(t)))*(xH(t)*cos(theta(t))+yH(t)*sin(theta(t)))/(xH(t)^2+yH(t)^2)^2.5, diff(xH(t), t) = vxH(t), diff(vxH(t), t) = -G*M*xH(t)/(xH(t)^2+yH(t)^2)^(3/2), diff(yH(t), t) = vyH(t), diff(vyH(t), t) = -G*M*yH(t)/(xH(t)^2+yH(t)^2)^(3/2); ICs := xH(0) = a*(1+e), yH(0) = 0, vxH(0) = 0, vyH(0) = sqrt(G*M*(1-e)/(a*(1+e))), theta(0) = `&theta;IC`, omega(0) = `&omega;IC`; soln := dsolve({Eqns, ICs}, numeric, maxfun = 0, output = array([seq(i, i = 0 .. n*Th, Th)])); plots:-odeplot(soln, [modp(theta(t), 2*Pi), omega(t)/`&omega;H`], 0 .. n*Th, labels = ["&theta;(t)","&omega;(t)/&omega;H"], axes = boxed, style = plottools:-point, size = [.25, .75]) end proc:

plots:-display(HyperionOrbit(.5, 1.8*`&omega;H`, 10));

Download Poincare_section_Boyd_plot_fixing_theta.mw

bia Man

First 1190 1191 1192 1193 1194 1195 1196 Last Page 1192 of 2429