MaplePrimes Questions

Hi all, I am new to MAPLE, I have been using Mathematica mostly. Here is what i am trying to do in MAPLE, 

Use the procedure plotmotion 2 on the plotmo worksheet to animate the motion of a marble in a bowl of the shape of the bottom half of the ellipsoid x^2/16 + y^2/12 +z^2/9 = 1

Can anyone help? It will be greatly appreciated.

Hello,

I'm trying to define a matrix that will have all the possible game results in a 9 team competition in one round. The results are a column in the form [game 1 home team, game 1 away team, game 2 home team, ... , game 9 away team]

This should give me a matrix with 18 rows, and 512 columns, where each column is a possible game result.

In essense, to continue the pattern of this matrix on for all 512 columns.

I've fiddled with for statements, but it ends up getting very messy given the number of nests that are required.

Thoughts? Any help would be much appreciated,

Thanks, Alex, Australia.

 

 

 

 

 

 

Hi Maple People,

I got stuck trying to write a procedue for my project today.  My topic is Prime Constellations.

If you want to help me, please read my primer on constellations

https://docs.google.com/viewer?a=v&pid=sites&srcid=ZGVmYXVsdGRvbWFpbnxwcmltZWNvbnN0ZWxsYXRpb25zfGd4OjNmOWZlNWQxZDY2ODY4YmQ

The procedure should refine an offset Vector and increase a multiplier scalar.

 

I want to delay putting the broken procedure until I can simplify it and try again.

Thanks,

Matt Anderson

hi friends

I encountered a problem and I can not draw the plot of this code

> sol := fsolve({diff(S, x) = 0, diff(S, y) = 0}, {x, y});



> with(VectorCalculus);
> with(linalg);
> s1 := evalf(subs(sol, linalg[grad](S, [x, y])));

> with(VectorCalculus);
> with(LinearAlgebra);
> s2 := evalf(subs(sol, linalg[hessian](S, [x, y]))); pmp0 := [x-subs(sol, x), y-subs(sol, y)]; sapprox := s0+evalm(`&*`(`&*`(transpose(pmp0), s2), pmp0));
> with(Statistics);
>
> with(stats); statevalf[icdf, chisquare[4]](.95);

> with(VectorCalculus);
> with(plottools);
> with(plots);
> with(linalg); ellips := {seq(stats*([statevalf[icdf, chisquare[4]]])(c) = sapprox, c = [.5, .95, .999])};
> plots(ellips(x, y), x = 950 .. 1000, y = 700 .. 750, grid = [50, 50], view = [950 .. 1000, 700 .. 750]);

 

 

can you helpe me?Thank you

Dear Friends,

I am solving 6 ODEs using maple15. then i got this error. anyone know abou this? thank you.

problem2.mw

 

 

restart:with (plots): B:=1:M:=1:Gr:=0.5:Pr:=3:w:=0.02:blt:=5:Bi:=10:

Eq1:=diff(f(eta),eta,eta,eta)-(diff(f(eta),eta))^(2)+f(eta)*diff(f(eta),eta,eta)+B*H(eta)*(F(eta)-diff(f(eta),eta))-M*diff(f(eta),eta)+Gr*theta(eta)=0;

diff(diff(diff(f(eta), eta), eta), eta)-(diff(f(eta), eta))^2+f(eta)*(diff(diff(f(eta), eta), eta))+H(eta)*(F(eta)-(diff(f(eta), eta)))-(diff(f(eta), eta))+.5*theta(eta) = 0

(1)

Eq2:=(1+Nr)*diff(theta(eta),eta,eta)+Pr*f(eta)*diff(theta(eta),eta)+(2/3)*H(eta)*B*(theta1(eta)-theta(eta))=0;

(1+Nr)*(diff(diff(theta(eta), eta), eta))+3*f(eta)*(diff(theta(eta), eta))+(2/3)*H(eta)*(theta1(eta)-theta(eta)) = 0

(2)

Eq3:=H(eta)*F(eta)+H(eta)*diff(G(eta),eta)+G(eta)*diff(H(eta),eta)=0;

H(eta)*F(eta)+H(eta)*(diff(G(eta), eta))+G(eta)*(diff(H(eta), eta)) = 0

(3)

Eq4:=F(eta)^2+G(eta)*diff(F(eta),eta)+B*(F(eta)-diff(f(eta),eta))=0;

F(eta)^2+G(eta)*(diff(F(eta), eta))+F(eta)-(diff(f(eta), eta)) = 0

(4)

Eq5:=G(eta)*diff(G(eta),eta)+B*(f(eta)+G(eta))=0;

G(eta)*(diff(G(eta), eta))+f(eta)+G(eta) = 0

(5)

Eq6:=G(eta)*diff(theta1(eta),eta)+l*B*(theta1(eta)-theta(eta))=0;

G(eta)*(diff(theta1(eta), eta))+l*(theta1(eta)-theta(eta)) = 0

(6)

bcs:=f(0)=0,(D(f))(0)=1,(D(theta))(0)=-Bi*(1-theta(0)),(D(f))(blt)=0,F(blt)=0,G(blt)=-f(blt),H(eta)=w,theta(blt)=0,theta1(blt)=0;

f(0) = 0, (D(f))(0) = 1, (D(theta))(0) = -10+10*theta(0), (D(f))(5) = 0, F(5) = 0, G(5) = -f(5), H(eta) = 0.2e-1, theta(5) = 0, theta1(5) = 0

(7)

L:=[0.5,1,1.5,2];

[.5, 1, 1.5, 2]

(8)

for k from 1 to 4 do p:=dsolve(eval({Eq1,Eq2,Eq3,Eq4,Eq5,Eq6,bcs},Nr=L[k]),[f(eta),F(eta),G(eta),H(eta),theta(eta),theta1(eta)],numeric,output=listprocedure);end do:

Error, (in dsolve/numeric/bvp) unevaluated names in system not allowed: {Y[9], Y[10]}

 

``

``

``

 

Download problem2.mw

 

When I take the derivative of this function wrt p, I am getting this:

 

Why the program gives , instead of only

Sorry for the format, I just copy and paste.

Thanks,

 

 

 

Hi,

I want to display several plots insequence and tried creating a list of polygonplot3d objects called "plotlist" and using them as an argument like this:

display(plotlist)

which triggered the error message:

Error, (in plots:-display) expecting plot structure but received: plotlist

Where's my mistake and how could I achieve the desired outcome? (displaying polygon3dplots in sequence)

 

Best Regards

 

 

I have generated a picture using pointplot and display that creates a grid of colors based on the angle between vectors and (1,0) in R^2. Each angle is given a unique color via the HSV color space. I would like to provide a legend explaining the colour, something along the lines of

red (brown) green (blue) red
0                Pi                 2*Pi

and possibly

red       (brown)       green 
angle1                     angle2   

Is there an easy way to do this, or will I have to mess around with a custom pointplot/textplot solution?

Good day everyone,

Please I do get numerical output/values in this solution U1.mw

Best regards

please help me a bout:


> u := am+bm*m0*x(t)+cm*n0*y(t);
> v := an+bn*m0*x(t)+cn*n0*y(t);
> eq1 := diff(x(t), t) = 2*A*(sinh(u)-x(t)*cosh(u));
> eq2 := diff(y(t), t) = 2*A*(sinh(v)-y(t)*cosh(v));
> init := x(0) = X, y(0) = Y;
> sol := dsolve({eq1, eq2, init}, {x(t), y(t)});
> param := A = 1/2, am = 0, bm = 1.2*exp(-4), cm = .5*exp(-5), m0 = 10000, an = 0, bn = -exp(-4), cn = 1.2*exp(-3), n0 = 1000;
> ;
> save u, v, param, " narm.sav";
Warning, unassigned variable `u` in save statement
Warning, unassigned variable `v` in save statement
Warning, unassigned variable `param` in save statement
> init := x(0) = 0, y(0) = 0.1e-1;
> eq1 := subs(param, eq1);
> eq2 := subs(param, eq2);
> sol := dsolve({eq1, eq2, init}, {x(t), y(t)}, numeric);
> with(plots);
> odeplot(sol, [x(t), y(t)], 0 .. 60, numpoints = 300, view = [-1 .. 1, -1 .. 1]);

How do I control colors in such Maple program using plot3d? Is it possible to define the conditions?

> restart: with(plots): mandelbrot := proc(x, y) local c, z, m; c := evalf(x+y*I); z := c; for m from 0 to 30 while abs(z) < 2 do z := z^2+c od; m
end: > plot3d(0, -2 .. 0.7, -1.2 .. 1.2, orientation=[-90,0], grid=[40, 40], style=point, scaling=constrained, color=mandelbrot);

I have a (251*1) vector for which I need to compute the cumulative sum. I'm currently using :

s:=Statistics[CumulativeSum](prices[1..251,1])

defined under the Statistics package but this method doesn't return an output, even though the argument passed is a vector.

Is there another method I can use?

In my  Standard GUI Maple 2015 (32 bit)  on Windows 8.1  plots[spacecurve]  command does not work:

plots[spacecurve]([cos(t), sin(t), t], t = 0 .. 2*Pi);

                         

 

Can someone confirm this bug?

dx/dt=2x(1-x/2)-xy, 

dy/dt=y(9/4 -y^2)-(x^2)y 

First 1291 1292 1293 1294 1295 1296 1297 Last Page 1293 of 2434