Items tagged with plot

Feed
sol_L := dsolve({de_L, ic});

    {x(t) = (-y0 - x0) exp(-2 t) + (y0 + 2 x0) exp(-t),
      y(t) = -2 (-y0 - x0) exp(-2 t) - (y0 + 2 x0) exp(-t)}
How i can plot this?thanks

O(t) = 9.2-[(8.93*35.3)/(11.17-8.93)]*(exp(-8.93*t)-exp(-11.17*t))-8.3*exp(-11.17*t)

t betwen 0 and 1

How to plot f(x)=x in the interval -2<x<2 with a period 4?

The plot should look like this

Aslam o Alikum!

Hope everyone going fine with sound health I want to make the graph of following system of ODEs

Eq1 := diff(F(eta), eta, eta, eta, eta)-M*(eta*(diff(F(eta), eta, eta, eta))+3*(diff(F(eta), eta, eta))+(diff(F(eta), eta))*(diff(F(eta), eta, eta))-F(eta)*(diff(F(eta), eta, eta, eta)))-Ha^2*(diff(F(eta), eta, eta));
Eq2 := diff(G(eta), eta, eta)+Pr*M*(F(eta)*(diff(G(eta), eta))-eta*(diff(G(eta), eta)))+Pr*Ec*(diff(F(eta), eta, eta))^2+Nb*(diff(H(eta), eta))+diff(G(eta), eta)+Nt*(diff(G(eta), eta))^2;
Eq3 := diff(H(eta), eta, eta)+M*Sc*(F(eta)*(diff(H(eta), eta))-eta*(diff(H(eta), eta)))+Nt*(diff(F(eta), eta, eta))/Nb;
 

using the BC's

IC1 := F(0) = 0, ((D@@2)(F))(0) = 0, (D(G))(0) = 0, (D(H))(0) = 0;
     F(0) = 0, @@(D, 2)(F)(0) = 0, D(G)(0) = 0, D(H)(0) = 0
IC2 := F(1) = 1, (D(F))(1) = 0, G(1) = 1, H(1) = 1;
 

for Ec := .1; Nt := .1; Nb := .1; Sc := .5; Pr := 10; M := .5; and different values of Ha like Ha:=0, Ha:=2, Ha:=4, Ha:=6 and Ha:=8 in one coordinate. 

I am waiting your positive response.

Hi

Here is my question: 

I have been trying to plot this on maple but i only get an ampry box

-arctan((2.m.x)/(1-x^2))

Simple as that... M is between 0 to 10

And x and y are supposed to fmbe from 0 to infinity(i dont know what to write for infinity so i give a large number like 2milions)

> restart;

with(plots);

pr := .72; p := 0; n := [.5, 1, 1.5]; s := 0; a := .2; b := 0; L := [red, blue, green]; l := 0; k := 1;

for j to nops(n) do R1 := 2*n[j]/(1+n[j]); R2 := 2*p/(1+n);

sol1 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R1*(1-(diff(f(eta), eta))^2) = 0, diff(diff(theta(eta), eta), eta)+pr*k*f(eta)*(diff(theta(eta), eta))+R2*pr*k*(diff(f(eta), eta))*theta(eta)+(2*(a*(diff(f(eta), eta))+b*theta(eta)))/(1+n[j]) = 0, f(0) = 1, (D(f))(0) = b*((D@@2)(f))(0), (D(f))(1.8) = 0, theta(0) = 1+s*(D(theta))(0), theta(1.8) = 1], numeric, method = bvp);

fplt[j] := plots[odeplot](sol1, [eta, diff(diff(f(eta), eta), eta)], color = L[j], axes = boxed); tplt[j] := plots[odeplot](sol1, [[eta, theta(eta)]], color = L[j], axes = normal) end do; plots:-display([seq(fplt[j], j = 1 .. nops(n))]);

plots:-display([seq(tplt[j], j = 1 .. nops(n))]);

 

staganation_point1.mw
 

can we chage the axis sir ?? like  f'' vs eta to f'' vs lambda.

``

restart

l := 1:

1

 

1.5

 

.5

 

[blue, green, red, yellow]

(1)

``

for j to nops(p) do R1 := 2*n/(n+1); R2 := 2*p[j]/(n+1); R3 := 2/(n+1); sol1 := dsolve([diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))+R1*(1-(diff(f(eta), eta))^2)-M*(diff(f(eta), eta)) = 0, diff(diff(theta(eta), eta), eta)+pr*f(eta)*(diff(theta(eta), eta))-R2*pr*(diff(f(eta), eta))*theta(eta)+R3*(A*(diff(f(eta), eta))+B*theta(eta)) = 0, f(0) = 1, (D(f))(0) = L+b*((D@@2)(f))(0), (D(f))(7) = 1, theta(0) = 1+s*(D(theta))(0), theta(7) = 0], numeric, method = bvp); plots[odeplot](sol1, [eta, ((D@@2)(f))(eta)], color = red); fplt[j] := plots[odeplot](sol1, [eta, f(eta)], color = K[j], axes = boxed); tplt[j] := plots[odeplot](sol1, [[eta, theta(eta)]], color = K[j], axes = normal); fplt[j] := plots[odeplot](sol1, [eta, diff(f(eta), eta)], color = K[j], axes = boxed) end do:

 

 

plots:-display([seq(fplt[j], j = 1 .. nops(n))]);

 

sol1(0)

sol1(0)

(2)

sol1(.1)

[eta = .1, f(eta) = 1.05958091104306206, diff(f(eta), eta) = .643210624614908300, diff(diff(f(eta), eta), eta) = .881482678165403044, theta(eta) = .623284688471349546, diff(theta(eta), eta) = -.578039450700496560]

(3)

sol1(.2)

[eta = .2, f(eta) = 1.12800452943200891, diff(f(eta), eta) = .722346769554029544, diff(diff(f(eta), eta), eta) = .706526135439307756, theta(eta) = .568123251856343492, diff(theta(eta), eta) = -.525530979400813946]

(4)

sol1(.3)

[eta = .3, f(eta) = 1.20351830506746449, diff(f(eta), eta) = .785511903074783246, diff(diff(f(eta), eta), eta) = .561442941644520022, theta(eta) = .518103974464032668, diff(theta(eta), eta) = -.475257424178228970]

(5)

sol1(.4)

[eta = .4, f(eta) = 1.28466826824405134, diff(f(eta), eta) = .835505660630676662, diff(diff(f(eta), eta), eta) = .442470716586289281, theta(eta) = .472985640642506311, diff(theta(eta), eta) = -.427567049032814172]

(6)

sol1(.5)

[eta = .5, f(eta) = 1.37026161183094430, diff(f(eta), eta) = .874752886901313142, diff(diff(f(eta), eta), eta) = .345911467377074400, theta(eta) = .432494259338694842, diff(theta(eta), eta) = -.382764248064397461]

(7)

sol1(.6)

[eta = .6, f(eta) = 1.36678221814533528, diff(f(eta), eta) = .771028661281065508, diff(diff(f(eta), eta), eta) = .407805382194403932, theta(eta) = .876413930517023876, diff(theta(eta), eta) = -.197648778495384870]

(8)

sol1(2)

[eta = 2., f(eta) = 2.66120522956795602, diff(f(eta), eta) = .991532161353848585, diff(diff(f(eta), eta), eta) = 0.251405465681268682e-1, theta(eta) = .635967939441598018, diff(theta(eta), eta) = -.144641270049362308]

(9)

``

``

``

 

``


 

Download staganation_point1.mw

 

 

 

 

 

dear all,

how can I plot

 u:=(3*x-1)*(1-x):
v:=2*x*(1-2*x):

in (u,v)-plane

where x from 0..1

 

Hi Maple expert,

May I ask your time little bit, I am struggling with the 3D plot in maple, for basic 3D plot in maple is no problem for me, but this is a new case for me actually I want to perform a 3D plot based on these 3 different type of plots (or equations), here they are:

1. plot(2*(exp(-0.15*(u)^1.5))*(u)^0.5,u=0..2);

2. plot([0.2*(x-sin(x)),0.5*(1-cos(x)),x=0..10*Pi]);

3. polarplot(4+1*cos(5*x),x=0..2*Pi);

Can someone help me for this, many thanks for your help and time.

i would like to plot a cube with the color given by a 3 variables functions (a way to plot in R^4)

in this way the cube is the domain of the functions, its color is the value of the function, i would like to put in evidence the max min of such function

but the command that i tried to use does not work:

plots[display](plotools([cuboid]([0,0,0],[1,1,1]),trasparency=.7,colorscheme["xyzcoloring",proc(x,y,z) option operator , arrow; x^2 + y^2 - z^2 end(proc)])

thanks for your help, anna rita

Hi, I have been trying to solve the Schrodinger equation for harmonic oscillators using dsolve and plot the the wavefunctions for the different energy levels. However I am struggling to plot all the different wavefuntions on the same plot. I also want to normalize the wavefunctions to help compare their shapes and values. Here's my code:- schro := {diff(psi(x), x, x)-(alpha*x^4+x^2-energy)*psi(x) = 0}; // d / d \\ / 4 2 \ \ { |--- |--- psi(x)|| - \alpha x + x - energy/ psi(x) = 0 } \\ dx \ dx // / ic := {psi(3) = 0, (D(psi))(3) = 1}; {psi(3) = 0, D(psi)(3) = 1} schro1 := subs(energy = 3.30687, alpha = .1, schro); soln1 := dsolve(schro1 union ic, {psi(x)}, type = numeric); // d / d \\ / 4 2 \ \ { |--- |--- psi(x)|| - \0.1 x + x - 3.30687/ psi(x) = 0 } \\ dx \ dx // / proc(x_rkf45) ... end; with(plots); [animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, densityplot, display, dualaxisplot, fieldplot, fieldplot3d, gradplot, gradplot3d, implicitplot, implicitplot3d, inequal, interactive, interactiveparams, intersectplot, listcontplot, listcontplot3d, listdensityplot, listplot, listplot3d, loglogplot, logplot, matrixplot, multiple, odeplot, pareto, plotcompare, pointplot, pointplot3d, polarplot, polygonplot, polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus, semilogplot, setcolors, setoptions, setoptions3d, shadebetween, spacecurve, sparsematrixplot, surfdata, textplot, textplot3d, tubeplot] odeplot(soln1, [x, psi(x)], -3 .. 3); Thank in advance

Hi,

Is there any way to generate a subplot in Maple2016?

I want to show many plots in one table at once.

I appreciate any idea you may have.

 

plz help! How to draw the three-dimensional graphics (p(x,t))? When I run,it can not run.I do not know where the problem lies.The code is as follows:

restart:
with(PDEtools); with(student); KN := 3;
C2 := 1/.3; C1 := 0.6e-2/(.3); C3 := 4.3/(.3); beta := 0.43e-1;
ADM1 := proc (n) options operator, arrow; convert(subs(lambda = 0, value((Diff(F(Sum(lambda^i*U[i], i = 0 .. n)), `$`(lambda, n)))/factorial(n))), diff) end proc; A0[0] := F(U[0]);
for n to KN do A0[n] := ADM1(n) end do;
for n from 0 to KN do A[n] := unapply(simplify(convert(C1*(diff(subs({seq(U[i] = p[i](x, t), i = 0 .. KN)}, expand(subs(F(U[0]) = U[0]*exp(U[0]), A0[n]))), t)), diff)+C2*convert(subs({seq(U[i] = Diff(p[i](x, t), x), i = 0 .. KN)}, expand(subs(F(U[0]) = U[0]^2, A0[n]))), diff)), x, t) end do;
p[0] := proc (x, t) options operator, arrow; .2*sin((1/2)*Pi*x)*exp((-1/4)*t*Pi*Pi/C1) end proc;
p[1] := proc (x, t) options operator, arrow; (-(int(subs(t = s, C3*(diff(p[0](x, t), x))), s = 0 .. t))+int(subs(t = s, A[0](x, t)), s = 0 .. t))/beta end proc;
p[2] := proc (x, t) options operator, arrow; -(int(subs(t = s, C3*(diff(p[1](x, t), x))), s = 0 .. t))+int(subs(t = s, A[1](x, t)), s = 0 .. t) end proc;
p := unapply(subs(x = Zeta, t = tau, p[0](x, t)+p[1](x, t)+p[2](x, t)), Zeta, tau);
with(plots);
plot3d(p(x,t),x=0..1,t=0..1)




Hello out there,

i just started using maple for some kinetik calculations.

it is a combination of the following formulas:

 

Lo := piecewise(sin(t) < 0, 0, sin(t) >= 0, sin(t))

L := Lo*(int(2*exp(-b*x(t)*(j*cos(phi)+sqrt(J^2-j^2*sin(phi)^2))), phi = 0 .. Pi))/Pi

r:= (int(u*L*j/(L+k), j = 0 .. J))/J^2

diff(x(t), t) = -B*x(t)+x(t)*r

 

The values would be:

 

b:=1
u:=1.2
J:=2.5
k:=3
B:=0.8

 

 

I want to plot x over t and get the values for it...
I tried several aproaches, but nothing realy worked...
Thanks a lot

I tried to generate a simple plot based on an expression using units. Please see the screenshot of the worksheet below

The y-axis should be scaled in m3 not m2 because a, b and the variable r are in m. I know that it is possible to enforce the axis scaling in units using useunits but since the mathematical definition of the expression is clear this should not be necessary.

The other question is why is the y-axis missing the variable name? On this page on the Online Help http://www.maplesoft.com/support/help/Maple/view.aspx?path=MaplePortal/UnitsInPlots it is mentioned "When you plot an expression, the axes labeling shows the variable names and units." I think what I am plotting is an expression.

 

Hi everyone, 

I am presently writting an alternative to the procedure GraphTheory:-DrawGraph for a better rendering.

I come across the following problem : the edges (CURVES) are in the foreground and the polygons (POLYGONS) which contain the vertex names are in the background.

I often face this kind of problem as I use PLOT either in the form PLOT( MyPlot1, MyPlot2, ...) or in the form 
plots:-display(PLOT(Myplot1), PLOT(Myplot2), ...)
How does Maple decide which structure (CURVES, POLYGONS, ..) is in front of or behind another one ?

I suppose there must be a tip to manage this ?

Thank you for your answers

5 6 7 8 9 10 11 Last Page 7 of 184