Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Here was my search strategy to find the answer:

keyboard shortcut execute worksheet on http://www.mapleprimes.com/search/?term=keyboard%20shortcut%20execute%20worksheet

only found one thing http://www.mapleprimes.com/users/erik10/answers

Googled it found 

http://www.maplesoft.com/support/help/maple/view.aspx?path=worksheet%2Freference%2Fhotwin

http://www.maplesoft.com/support/help/Maple/view.aspx?path=worksheet/expressions/executeworksheet

http://scipp.ucsc.edu/~haber/webpage/maple_cheatsheet.pdf

Don't see a keyboard shortcut. 

What is the keyboard shortcut to execute entire worksheet? Is there a name for this action of executing the whole worksheet? 

 

What did you search and click to find the answer? 

Hello everybody,

 

I find a answer on that subject at this link: http://www.mapleprimes.com/maplesoftblog/95508-Monte-Carlo-Pi

But for starting easy, I try this:

with(RandomTools[MersenneTwister]);
directpi := proc (n)
local i, x, y, N;
N := 0;
Digits := 2;
for i to n do
x := GenerateFloat();
y := GenerateFloat();
if x^2+y^2 < 1 then
N := N+1:
end if:
end do:
N;
end proc;


evalf(directpi(4000)/(1000.0), 4);

                             3.176

I would like to know if this code is appropriate and if it is possible to have a random number between -1 to 1 so it will be easier to plot the points in the square with the circle inside.

 

Thank you.

--------------------------------------
Mario Lemelin
Maple 18 Ubuntu 13.10 - 64 bits
Maple 18 Win 7 - 64 bits messagerie : mario.lemelin@cgocable.ca téléphone :  (819) 376-0987

I am Maple 15 student edition user and the packages is not loaded when I use the execute all command from the GUI. I have to manually press the enter key to have by package loaded before I use a function defined in the package. How can I work around this issue?

Thanks

P.S: You can see the Maple document attached to this post.

NumericalSimulation.mw

A := Matrix(3, 4, {(1, 1) = 1, (1, 2) = 0, (1, 3) = 0, (1, 4) = 0, (2, 1) = 2, (2, 2) = 2, (2, 3) = 0, (2, 4) = 0, (3, 1) = 0, (3, 2) = 6, (3, 3) = 4, (3, 4) = 0});
print(`output redirected...`); # input placeholder
Matrix(%id = 4436610626)
RowSpace(A);
print(`output redirected...`); # input placeholder
RowSpace(Matrix(%id = 4436610626))

 Why is this failing to give the row space? I copied and pasted it from the tutorial on this and it seems to run in the tutorial but not when I copy and paste it outside the tutorial? 

tutorial http://www.maplesoft.com/support/help/maple/view.aspx?path=LinearAlgebra%2FRowSpace

See picture

http://i.imgur.com/eQtEJU4.png?1

 

In general what is the best way to share maple code with users on this forum ? The above copy and paste looks messy. 

 

I'm trying to solve this system of ODEs by Laplace transform. 

> de1 := d^2*y(t)/dt^2 = y(t)+3*x(t)

> de2 := d^2*x(t)/dt^2 = 4*y(t)-4*exp(t)

with initial conditions 

> ICs := y(0) = 2, (D(y))(0) = 3, x(0) = 1, (D(x))(0) = 2

 

Using 

> deqns := de1, de2

and

> var := y(t), x(t)

 

I need to solve it for both y(t) and x(t), I have tried this by:

> dsolve({ICs, deqns}, var, method = laplace)

And

> dsolve({ICs, deqns}, y(t), method = laplace)

> dsolve({ICs, deqns}, x(t), method = laplace)

 

However I get this error message:

Error, (in dsolve/process_input) invalid initial condition

 

Any help is appreciated

Hi Maple friends.

Maple tends to spit out results(which comprise of variables) in very complicated forms, and I have to use the context menu to select 'simplify' to reduce them.

Is there a setting which will automatically simplify Maple's output?

Thanks in advance.

test.mw

can someone please look at this modified C.Love code, hopefully the problem i'm having is self explanatary.

In Maple 15 is it possible to differentiate with respect to say x - y instead of x

 

Thanks

Someone asked on math.stackexchange.com about plotting x*y*z=1 and, while it's easy enough to handle it with implicitplot3d it raised the question of how to get nice constained axes in the case that the x- or y-range is much less than the z-range.

Here's what WolframAlpha gives. (Mathematica handles it straight an an plot of the explict z=1/(x*y), which is interesting although I'm more interested here in axes scaling than in discontinuous 3D plots)

Here is the result of a call to implicitplot3d with default scaling=unconstrained. The axes appear like in a cube, each of equal "length".

 

Here is the same plot, with scaling=constrained. This is not pretty, because the x- and y-range are much smalled than the z-range.

 

How can we control the axes scaling? Resizing the inlined plot window with the mouse just affects the window. The plot itself remains  rendered in a cube. Using right-click menus to rescale just makes all axes grow or shrink together.

One unattractive approach it to force a small z-view on a plot of a much larger z-range, for a piecewise or procedure that is undefined outisde a specific range.

plots:-implicitplot3d(proc(x,y,z)
                        if abs(z)>200 then undefined;
                        else x*y*z-1; end if;
                      end proc,
                      -1..1, -1..1, -200..200, view=[-1..1,-1..1,-400..400],
                      style=surfacecontour, grid=[30,30,30]);

Another approach is to scale the x and y variables, scale their ranges, and then force scaled tickmark values. Here is a rough procedure to automate such a thing. The basic idea is for it to accept the same kinds of arguments are implicitplot3d does, with two extra options for scaling the axis x-relative-to-z, and axis y-relative-to-z.

implplot3d:=proc( expr,
                  rng1::name=range(numeric),
                  rng2::name=range(numeric),
                  rng3::name=range(numeric),
                  {scalex::numeric:=1, scaley::numeric:=1} )
   local d1, d2, dz, n1, n2, r1, r2, rngs, scx, scy;
   uses plotfn=plots:-implicitplot3d;
   (n1,n2) := lhs(rng1), lhs(rng2);
   dz := rhs(rhs(rng3))-lhs(rhs(rng3));
   (scx,scy) := scalex*dz/(rhs(rhs(rng1))-lhs(rhs(rng1))),
                scaley*dz/(rhs(rhs(rng2))-lhs(rhs(rng2)));
   (r1,r2) := map(`*`,rhs(rng1),scx), map(`*`,rhs(rng2),scy);
   (d1,d2) := rhs(r1)-lhs(r1), rhs(r1)-lhs(r1);
   plotfn( subs([n1=n1/scx, n2=n2/scy], expr),
           n1=r1, n2=r2, rng3, _rest[],
           ':-axis[1]'=[':-tickmarks'=[seq(i=evalf[3](i/scx),i=r1,d1/4)]],
           ':-axis[2]'=[':-tickmarks'=[seq(i=evalf[3](i/scy),i=r2,d2/4)]],
           ':-scaling'=':-constrained');
end proc:

The above could be better. It could also detect user-supplied custom x- or y-tickmarks and then scale those instead of forming new ones.

Here is an example of using it,

implplot3d( x*y*z=1, x=-1..1, y=-1..1, z=-200..200, grid=[30,30,30],
            style=surfacecontour, shading=xy, orientation=[-60,60,0],
            scalex=1.618, scaley=1.618 );

Here is another example

implplot3d( x*y*z=1, x=-5..13, y=-11..5, z=-200..200, grid=[30,30,30],
            style=surfacecontour, orientation=[-50,55,0],
            scaley=0.5 );

Ideally I would like to see the GUI handle all this, with say (two or three) additional (scalar) axis scaling properties in a PLOT3D structure. Barring that, one might ask whether a post-processing routine could use plots:-transform (or friend) and also force the tickmarks. For that I believe that picking off the effective x-, y-, and z-ranges is needed. That's not too hard for the result of a single call to the plot3d command. Where it could get difficult is in handling the result of plots:-display when fed a mix of several spacecurves, 3D implicit plots, and surfaces.

Have I overlooked something much easier?

acer

Hello,

 

I have a function, lets say g(x,y,...), that depends on many other functions. But I don´t want the results that are inside certain intervals, and I need to receive the results of those functions as something like NULL when asking for a result that is inside any of those intervals. That way, g(x,y,...) would also have to result in something like NULL if any of the lesser functions are NULL fora any given values.

 

I tryied using the piecewise command, and for the intervals that I wanted, it worked, but for those I wanted to be NULL, they were understood as 0, and so G(x,y,..) continued to exist but with a very different value.

 

To clarify what I need, I will try an exemple:

 

Imagine I have the function f(x)=x

 

I want to disconsider the results for x<2 and x>6, in a way that if I try the command 'f(1)', I will receive something like NULL and know that it is outside the range.

 

In the same way, I need the plot of this function f(X) to show the function only from 2 to 6, but not existing for the delimited intervals.

 

Ad if I continue and make g(x)=f(x)+10 , I don´t want g(x) to exist if f(x) doesn´t exist, and same for the g(x) plot, which shouldn´t be shown in the intervals where f(x) don´t exist.

 

 

Thank you very much for your atention!

 

 

Hello,

 

could you help me solve this error ? I don't understand what it means.

 


> eq3:=diff(x(t),t,t)+Gamma*diff(x(t),t)+omega[0]^2*(x(t)-(diff(x(t),t,t)+Gamma*diff(x(t),t)+omega[0]^2*x(t)+omega[0]^2*X[0])/omega[0]^2) = -omega[0]^2*X[0]:
> dsolve(eq3);
Warning, it is required that the numerator of the given ODE depends on the highest derivative. Returning NULL.

 

Thanks.

I remember to have seen and used a command to make the graph of the output of FeynmanDiagrams but I can not find more sample files. Someone can tell me how to do (plot a Feynman graph using the result of FeynmanDiagrams).

Thanks and sorry for my english.

Dear all,

I am trying to find the intial velocity of a ball that is shot under an angle while only the start and end coordinates are given. The air resistance should also be taken into account. 

In order to do that I have build the following Maple sheet:

Assignment_question_1.mw 

I have used two differential equations that both include the variables v0 and t, and then try to solve them. Only I receive an answer in the form of RootOf, which I cannot remove with for example allvalues. 

I have been working for quite a long time on this but I am not coming any further, so is there anyone who can find what I am doing wrong/what I should be doing else? Or maybe my whole approach is not right?

Even a small step in the right direction would be appreciated a lot!

Thanks in advance,

Elise

How do I permanently change the current directory on a MAC computer?

 

Kind regards

 

Per Kirkegaard

Hi all

The aim of following program is minimization but it is unable to produce it. where is the mistake?

Taylor2.mws

thanks a lot.

Mahmood   Dadkhah

Ph.D Candidate

Applied Mathematics Department

First 1307 1308 1309 1310 1311 1312 1313 Last Page 1309 of 2223