Items tagged with graph graph Tagged Items Feed

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 (`θIC`, `ωIC`, n) local a, Mh, Msat, G, e, beta, M, Eqns, ICs, soln; option remember; global `ω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) = `θIC`, omega(0) = `ω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)/`ωH`], 0 .. n*Th, labels = ["θ(t)","ω(t)/ωH"], axes = boxed, style = plottools:-point, size = [.25, .75]) end proc:

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

Download Poincare_section_Boyd_plot_fixing_theta.mw

bia Man

Hello maple experts,

please is there any other way of solving this problem Wht.mw

Some complained that the graphs do not satisfy the boundary conditions smoothly.

Thanks.

VIM.3.mwHi,How do I make a graph and error table for exact and by using this method i.e VIM 

I am trying to design a question, where studenst have to find the equation of a linear function given its graph. For this to work, I have to be able to draw gridlines, which for some reason is not possible in Maple TA with the ordinary plot function. I know it should be possible using an applet; given the problems this is going to create with all the different browsers, however, this is not a viable solution. I have found another solution which is unreasonably complex but should work. The problem is that sometimes it doesn't work. When I press Refresh algorithm preview in the Question Designer I get "Broken Maple plot. Verify your plot statement" roughly every third time I refresh the preview. Any help would be greatly appreciated.

 

The algorithm is here:

$a1=range(-4,4,1);
$a2=range(1,4,1);
condition:not(eq($a1,$a2));
$px=range(0,5,1);
$y0=range(0,5,1);
$x0=if(gt($a1,0),$px,-$px);
$f=$a1/$a2*(x-$x0)+$y0;
$G=plotmaple("
 plots[display](
  [plots[coordplot](cartesian, [-10 .. 10, -10 .. 10], grid = [11, 11], color = grey),
   plot($f, x = -10 .. 10, y = -10 .. 10, color = [green, blue, red], thickness = 2, tickmarks = [[1], [1]], labels = [``, ``]),
   plots[textplot]([0, 10, 'y'], align = {above, right}),
   plots[textplot]([10, 0, 'x'], align = {above, right})
  ],
  labels = [``, ``], axes = normal, view = [-10 .. 10, -10 .. 10]),
 plotoptions='width=350,height=350'");
$g=gcd($a1,$a2);
$at=$a1/$g;
$an=$a2/$g;
$a=if(eq($an,1),$at,$at/$an);
$ans=$a*x+($y0-$a*$x0);

Hi, Everyone.

I'm trying to figure out how to make a graph (actually a "line") change color based on certain criteria so I can create a Maple Application.

I've figured out the following, which indicates the sign of the graph by the color, but it's way more complicated than it ought to be.

It's also not robust as I need it to be, this happens to depend on the zeros occuring at integer values.  I could deal with that on a case by case basis, but not for something automatically.

I'm at a loss; maybe I'm not Googling the right things... maybe Maple can't do this (although I'd heard color tools were being added in Maple 16).

Does anyone have any suggestions?

Thanks In Advance!

Joseph

In the following worksheet, I can evaluate S(32), but not S(31). As a result, I can't use any value that depends on S(31) (or any i<31), for example, the graph is incomplete. How can I get Maple to evaluate S(31)?

restart

sigma := 167:

t := proc (i) options operator, arrow; exp(-(1/2)*i^2/sigma^2) end proc:

S := proc (l) options operator, arrow; 1+2*(sum(t(i), i = l .. infinity)) end proc:

evalf(S(32))

356.9784188

(1)

evalf(S(31))

1.+2.*(sum(exp(-(1/55778)*i^2), i = 31 .. infinity))

(2)

Expected result of evalf(S(31)):

evalf(S(32)+2*t(31))

358.9442559

(3)

plot(S(y), y = 20 .. 50)

 

``


All told, S(i) fails for i in [0..31], i = 115,116, and 221. S(i) succeeds for all other i < 3500.

Download test.mw

 


> restart;
> with*plots;
> Eq1 := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta, eta))-(diff(f(eta), eta))^2 = 0; 

> N := 1;

> blt := 10;
> Eq2 := (diff(theta(eta), eta, eta))/Pr+f(eta)*(diff(theta(eta), eta)) = 0; 
> bcs1 := f(0) = 0, (D(f))(0) = 1, (D(f))(blt) = 0;
> bcs2 := (D(theta))(0) = -N*(1+theta(0)), theta(blt) = 0;
> L := [2.5, 3, 5, 7, 10];
> for k to 5 do R := dsolve(eval({Eq1, Eq2, bcs1, bcs2}, Pr = L[k]), [f(eta), theta(eta)], numeric, output = listprocedure); X1 || k := rhs(R[3]); X2 || k := rhs(R[4]); Y1 || k := rhs(R[5]); Y2 || k := -rhs(R[6]) end do;

 

 how I will draw the graph for Pr against theta   for Pr=2.5 until 7  taking rest of the parameter fix

 

Dear  i want to plot the graphs for different values of R[2,5,6.5]   here is the codes: Either ido loop or array i have to use or other procedure 

The graph display should be at R=2, one graph, R=5 another and R=6.5 

onemore graph

 

restart:
with(DEtools):
eta:=1+k*x+epsilon*sin(2*Pi*x):
A1(x):=-exp(-alpha*x)*J^2*R/(2*eta^3+6*xi*J^2*eta^2):

psi0:=A1(x)*y^3:
psi:=delta*psi0:
V:=-diff(psi,x):
delta:=0.1:
epsilon:=0.01:
alpha:=1:
xi:=0.001:
k:=0.1:
As:= [0, 2, 4]:
x:=0.2
plot([seq]([eval(V, J= A), y, y= 0..eval(eta, J= A)], A= As),title=" velocity at R=2.0 ",labels= ["v", "y"],color= [green, red, blue], linestyle= [solid, dash, dot],legend= [seq](J = A, A= As), axes=boxed);

 

 

 

 

 

 

i want to plot the graph based on this equation to know the stability of this equation :

y := A*(1/x+x*exp(-2*sqrt(-1)*b))+4*sin(h)^2*(2*exp(-sqrt(-1)*b)-3*sin(h)^2*x^(-sin(h))*exp(sqrt(-1)*b*(-sin(h)-1))+3*x^(-sin(h))*exp(sqrt(-1)*b*(-sin(h)-1)))/(3*(1-r))-exp(-2*sqrt(-1)*b)/x-x;

with A = (1+r)/(1-r), r = (1/3)*sin(h)^2, b = m*Pi*h, m = 1, h=0.05

can somebody help me..??

Hi, I was told there are 3 errors about this code, could anyone please help ?

Thank you very much 


DiSpanningTree := proc(G,v)
local comp,edges,recurse,visited;
recurse := proc(w)
local d,z;
if assigned(visited[]) then return false; end if;
comp(numelems(comp)+1) := w;
visited[w] := true;
d := Departures(G,v);
for z in d do
if recurse(z) then edges(numelems(edges)+1) := [w,z]; end if;
end do;
return true;
end proc;
visited := table();
comp := Array([]);
edges := Array([]);
recurse(v);
return Graph(convert(comp,'list'),edges);
end proc:

When I print a Table with graphs in it, the graphs are compressed vertically to an extreme extent. There has to be some compression from the screen, of course, since I have a wide-screen (although the actual worksheet does not use the whole screen width because of the palette and the fact that I don't use full-screen mode. However, the vertical compression  (that is, in a direction where there is plenty of room on the page) is many times more than the horizontal compression. In order to get the printed graph to look normal, I have to stretch the screen graph vertically an extreme amount, which, of course, looks horribly distorted on the screen.

Also when printing a graph as part of a Table, the graph area is clipped so that the title and legend do not print.

I thought that the Table struction would be a good way to display related graphs. But this behavior is making the use of Tables tedious and much less satisfactory because of the missing titles and legends.

None of these things happen when I print a graph in a regular execution group outside a Table structure.

Has anybody else had this problem, and have you found a way to avoid the compression and get the legends?

Dear all,

Is there a Maple command to get ColumnGraph in 3D as bar3 in Matlab.

In 2D, I use:

Statistics[ColumnGraph](A);  # equivalent to bar(A) in Matlab

Thanks

how can you graph when piecewise function is not continuous?

 

example:

BB := piecewise (-1.57 < c and c < -1.56, h, -0.06< c and  c < -0.05, aa, -0.5< c and c < 0.04, bb);

 

 

how to graph in maple 

for example

 

-2 < x < -3, h

-1 < x < -2, b

 

why do I get the error Error, (in rtable/Sum) invalid arguments

In positive numbers, I get it ok

Hello,

What is the combination key for zooming in/out on a Maple 3D graphic?

It seems I have forgotten... (CTRL+left mouse button does not work)

Thankx

JMC

1 2 3 4 5 6 7 Last Page 1 of 48