Items tagged with plots plots Tagged Items Feed

A quick question, Can I add legends to 3D plots? 

 

I have:

r:=t-><4cos2t+4sint,5cos3t+4cost,(sin3t+2sint)^2>

with(plots):

 

then I tried:

spacecurve(r(t), t = 0 .. 2*Pi, thickness = 2, color = black, axes = normal, labels = [x, y, z], numpoints = 150)

But I keep getting the error: Warning, unable to evalute the function to numeric values in the region. 

But I thought <4cos2t+4sint, 5cos3t+4cost, (sin3t+2sint)^2> was defined on the region t=0 to 2Pi??? 

I'm using Maple 2015 if that helps. 

Hi all,

I've got a parametric curve of the form f(x),h(x) x=0..1, and want to make an animation of this curve when changing an exogenous variable h0.

My attempt was to do it like this:

plots[animate](plot, [f(x), h(x), x = 0 .. 1], h0 = 0 .. 1)

but it doesn't work.

Does anyone know if/how this can be done?

Thanks,

Paul

would like to draw the graph of x^2-1

with below sys instead of x^2-1 

sys := {-(1/2)*(-x-1+sqrt(-3*x^2-2*x-3))*(diff(y(a, b), b))/x^2+(diff(y(a, b), a))*(diff(y(a, b), b)), (1/2)*(x+1+sqrt(-3*x^2-2*x-3))*(diff(y(a, b), a))/x^2+(x^2+x+1)/x^4, (1/2)*(x+1+sqrt(-3*x^2-2*x-3))*(diff(y(a, b), a))/x^2-(1/2)*(-x-1+sqrt(-3*x^2-2*x-3))*(diff(y(a, b), b))/x^2+(x^2+x+1)/x^4};
IBC:={x=1,y(a,3)=(3-1)*(3+1),y(0,b)=0,y(1,b)=0}:
pds:=pdsolve(sys,IBC,numeric,spacestep=140);
plots[display]([seq(pds:−plot(v,t=i10),i=0..5)]):

 

I never expected that the reflected light direction of sun from moon in the sky would be so dificult to imagine ...

at the following article mentioned :

we derive an equation for the magnitude of the moon tilt illusion that can be applied to all con gurations of sun and moon in the sky.

THE MOON TILT ILLUSION

 

since the calculations contains many steps and high level mathematical formula , there is no way rather to recourse to maple (powerful math assistant )

I hope there was adaptations between a lots of functions and predefined schema of maple and this problem so that the calculations and visualization facilitated several times ?

your effort will be a graet present for all the people of the world that look to the moon crescent everytime !

 

Hello,

I would like to plot an non coupled non linear oscillator.

The equations are the following:

K:=Matrix([<0, -1, 1, -1>,<-1, 0, -1, 1>,<-1, 1, 0,-1>,<1, -1, -1,0>]);
omega[sw]:=beta/(1-beta)*omega[s];
for i to 4
do
r[i]:=sqrt((u[i](t)/(L/2))^2+(v[i](t)/H)^2):
omega[i]:=omega[st]/(1+exp(b*v[i](t)))+omega[sw]/(1+exp(-b*v[i](t))):
Equ[i]:=diff(u[i](t),t)=Au*(1-r[i]^2)*u[i](t)+omega[i]*(L/2)/H*v[i](t):
Eqv[i]:=diff(v[i](t),t)=Av*(1-r[i]^2)*v[i]+omega[i]*(L/2)/H*v[i](t)+MatrixVectorMultiply(K,<seq(v[i](t),i=1..4)>)[i]:
EqSys[i]:=[Equ[i],Eqv[i]]:
end do:

My parameters are the following:

paramsGeo:=L=0.015,H=0.015,beta=0.5,Vf=0.3;
omegaS:=eval(Pi*Vf/L, [paramsGeo]);
paramsCycle:=omega[s]=omegaS,Au=1,Av=1,b=100;
params:=paramsGeo,paramsCycle;

I'm not sure with my initial equations. But, may be it is possible to start with:

ic:=[u[1](0)=0.8, v[1](0)=0,u[2](0)=0.8, v[2](0)=0,u[3](0)=0.8, v[3](0)=0,u[4](0)=0.8, v[4](0)=0];

For these equations, I would like to obtain the following plots:
- plot 1: horizontal axis : u[1](t) vertical axis : v[1](t).
- plot 2: horizontal axis : u[2](t) vertical axis : v[2](t).
- plot 3: horizontal axis : u[3](t) vertical axis : v[3](t).
- plot 4: horizontal axis : u[4](t) vertical axis : v[4](t).
- plot 5: horizontal axis : t, vertical axis : v[1](t), v[2](t), v[3](t), v[4](t).

For this last plot, I would like to obtain this kind of curve:

I image that since my equations are coupled i can not use directly use Deplot function but Dsolve.

May you help me for defining a good syntax for solving my system and then deducing the following plots?

Thanks a lot for your help

Dears,

I would like to draw a spectrum of plots. For example sin(mx) when m ∈ [1,2] and x ∈ [0,2] that is a spectrum of sinous plots. Obviously, if I write

plot(sin(m*x), m = 1 .. 2, x = 0 .. 2)

I will receive the error

Warning, expecting only range variable m in expression sin(m*x) to be plotted but found name x

I expect that there exist an alternative plot command in Maple, able to do that. I checked some members of Plots family that their name conducted me to check them. Moreover, I searched the web: but the most common retrieved results are related to drawing the area between two plots or under one plot.

A real example of what I am looking for is (not exaclty but something like) the below picure

A plot of Gaussian plots

Could you please let me know the alternative plot command, satisfying my requirement?

Hello everyone, I'm a new one to Maple, I've just learnt some basic tools.

I want to creat a command that can animate the graph of line y=ax+b by the parameter a, and b will be subscribe later. For example, I can plot y=x+b by:

F:=b->plot(x+b,x=-1..1):

display(F(1));

It did work.

However, applying this with animation didn't seem to work. 

F:=b->animate(plot,[a*x+b,x=-1..1],a=0..10,frames=5):

display(F(1));

It did not create an animation, instead 5 frames of this graph for a=0, 2.5, 5, 7.5, 10

Please show me a solution for this problem, thank you

 

hi all

i have tow plots and i want to display them together, but i cant, plz help me

restart;

I5:=Int((x^2)/((x^4-2*x0*x^2+x0^2+1)^(3/2)),x=0..infinity);
I6:=(1/2)*Int(1/(x^4-2*x0*x^2+x0^2+1)^(1/2),x=0..infinity);
f:=((-4/Pi)*((x0*I6-I5)/(x0*I5+I6)^(1/3)));
B:=x0/(x0*I5+I6)^(2/3): #B:=mu/eF
C:=B/((16/Pi^2)*((x0*I6-I5)^2/(x0*I5+I6)^(2/3))): #C:=mu/epsilon0/2

p1:=plot(B,x0=0..4):
p2:=plot(C,x0=-3..0):
display({p1,p2},x0=-3..4);

 

Hi all,

Thanks to help me to plot my result of dsolve for different values of a parameter.

With plots:-animate, I get one curve according to one value of the parameter. But, I need all curves in the same figure.

ws1.mw

 

Can anyone offer an explanation why the following error occurs, and how it can be avoided? Thanks!

 

plots:-shadebetween(2, 1/sqrt(-x^2+1), x = 0 .. (1/2)*sqrt(3))

Error, (in t) invalid input: nops expects 1 argument, but received 0

 

``

 

Download error.mw

  restart: with(plots):
  H := a -> piecewise(a>=0,1):
  f1 := y->(H(y-1*Pi)-H(y-2*Pi))*sin(y)^2:

  g1:=(f1(y)/sqrt(4*Pi*t))*(sin((x-y)^2/4/t+Pi/4)-sin((x+y)^2/4/t+Pi/4)):
  g2:= int(g1, y= 0..100):

  g3:= diff(g2,t):
  g4:= diff(g2,x$2):
  g5:= (g3^2+g4^2)/2:
  E2:= unapply(Int(g5, x= 0..100, epsilon= 1e-4, digits= 7), t):
  CodeTools:-Usage( plot(E2, 0..20, numpoints= 50, labels= [t, E]));

The above mentioned code should give constant figure, but it takes a lot of time and not accurate.

If you can help me to improve these codes, I would be pleased.

 

Thanks!