Items tagged with plotting


I have an two arrays of elements Op1[i,j] with -5 <= i,j <= 5 and Op2[i,j] with 5<=i,j<= 35 . When I use surfdata like


it produces a plot with both the range of Op1 and Op2 between 0 and 1 (see Fig1).

I tried surfdata({0p1,Op2},-5..35,-5..35) but this only gave the range of the figure from -5 to 35, -5 to 35 but with both plots scaled so that they cover the whole range(see Fig2).

How do i make surfdata use different ranges on the same plot for both arrays?

I found from this forum that to plot a 2D array of points use can be made of the Maple procedure surfdata.

Does anyone have suggestions on how to plot contours in (preferrably) Maple 16 or Maple 17?

I tried the following

Output := Array(-10 .. 10, -10 .. 10, proc (i, j) options operator, arrow; i^2+j^2 end proc):
F := proc (x, y) -> x^2+y^2 end proc:
surfdata(Output, color = F, dimension = 2);

but "the option dimension = 2" is a Maple 18 addition.

Ideally, I would like also to be able to plot contours with options found in the procedure

Hello every body i need help with my program. how can i plot r&s versus each other?




with(PDEtools, casesplit, declare)

[casesplit, declare]


with(DEtools, gensys)



declare((Omega, H)(z));

Omega(z)*`will now be displayed as`*Omega


H(z)*`will now be displayed as`*H


`&epsilon;` := -1



h := 70



r := 0.3e-3



n := 42





ode1 := -(diff(Omega(z), z))+Omega(z)*((3*(1-2*`&epsilon;`*h*sqrt(r)/H(z)-Omega(z)+2*Omega(z)^(3/2)/(3*n)))*(1+h*`&epsilon;`*sqrt(r)/sqrt(h^2*r+H(z)^2*(1-2*`&epsilon;`*h*sqrt(r)/H(z))))-2*sqrt(Omega(z))/n)/(1+z) = 0:

ode2 := -(diff(H(z), z))-(3*H(z)*(1/2))*(1-2*`&epsilon;`*h*sqrt(r)/H(z)-Omega(z)+2*Omega(z)^(3/2)/(3*n))*(1+h*`&epsilon;`*sqrt(r)/sqrt(h^2*r+H(z)^2*(1-2*`&epsilon;`*h*sqrt(r)/H(z))))/(1+z) = 0:

q := -1+(3/2*(1-2*`&epsilon;`*h*sqrt(r)/H(z)-Omega(z)+2*Omega(z)^(3/2)/(3*n)))*(1+h*`&epsilon;`*sqrt(r)/sqrt(h^2*r+H(z)^2*(1-2*`&epsilon;`*h*sqrt(r)/H(z)))):


sys := {ode1, ode2}

ics := {H(0) = h, Omega(0) = .77}:
sol := dsolve(`union`(sys, ics), numeric, output = listprocedure, stiff = true):``






r2 := q+2*q^2+(1+z)*(diff(q, z)):



s := (r2-1)/(3*(q-1/2)):












Hi everyone,


I am trying to plot the standard map equations which is defined by

p_{n+1}=p_{n}+k sin x_{n}  (mod 2pi)
x_{n+1}=x_{n}+p_{n}+k sin x_{n} (mod 2pi).


Is there is any simple code for plotting this map in Maple? 

Thank you in advance.

Let A and B be regions in space defined through f(x,y,z) > 0 and g(x,y,z) > 0, respectively.  How does one plot the surface of the object A ∩ B in Maple?

In Mathematica there is a command RegionPlot3D for that.  See sample here.  Perhaps there is an equivalent one in Maple that I don't know of.

How can I ask Maple to plot intersection of two implicitplot3d? It is explained how to draw the union in the Maple help by simply entering a list to combine plots, but I didn't see anything about intersection.

I don't recall ever changing any settings in my copy of Maple 16. But for some reason, I tried to make a plot and it spit out line-printed output rather than a regular plot.

How do I change this back? I have no use for line-printed output.

plot([E1(x), E2(x), E3(x)], x = 0 .. 5, title = "Nitrogen Nucleus Second-Order Energy Shifts")

                  Nitrogen Nucleus Second-Order Energy Shifts                 
  +                                                                     HHHHH
6e-08                                                               HHHHH     
  +                                                            HHHHHH         
4e-08                                                     HHHHHH              
  +                                                  HHHHHH                   
  +                                           HHHHHHH                         
2e-08                                  HHHHHHHH                               
  +                           HHHHHHHHHH                                      
  +             HHHHHHHHHHHHHHH                                               
0 +              1    HHHHHHHHHH2HH           3              4              5
  +                               HHHHHHHHH                                   
-2e-08                                    HHHHHHHH                            
  +                                              HHHHHHH                      
-4e-08                                                 HHHHHH                 
  +                                                         HHHHHH            
  +                                                              HHHHH        
-6e-08                                                                HHHHH   
  +                                                                       HHH







I have a second problem which is related to converting horizontal axis to degrees ( or in terms of Pi :  0 - Pi/2 ).

Thank you for your helps


s1 := (1/3*(4*(1.66+10^(-1))^2-x^2))^(3/2)/(x^2*(1+.66+10^(-1)))-sin(2*a);
s2 := (1/3*(4*(1.66+10^(-2))^2-x^2))^(3/2)/(x^2*(1+.66+10^(-2)))-sin(2*a);
s3 := (1/3*(4*(1.66+10^(-3))^2-x^2))^(3/2)/(x^2*(1+.66+10^(-3)))-sin(2*a);
k1 := series(s1, x);
       4.769067180 x     + (-0.5773502688 - 1. sin(2. a))

                           2    / 4\
          + 0.01164914507 x  + O\x /
k2 := series(s2, x);
       4.293792441 x     + (-0.5773502692 - 1. sin(2. a))

                           2    / 4\
          + 0.01293857500 x  + O\x /
k3 := series(s3, x);
       4.247636750 x     + (-0.5773502693 - 1. sin(2. a))

                           2    / 4\
          + 0.01307916821 x  + O\x /
p1 := convert(k1, polynom);
  4.769067180                                                2
  ----------- - 0.5773502688 - 1. sin(2. a) + 0.01164914507 x 
p2 := convert(k2, polynom);
  4.293792441                                                2
  ----------- - 0.5773502692 - 1. sin(2. a) + 0.01293857500 x 
p3 := convert(k3, polynom);
  4.247636750                                                2
  ----------- - 0.5773502693 - 1. sin(2. a) + 0.01307916821 x 
with(plots, implicitplot);
implicitplot([p1, p2, p3], a = 0 .. (1/2)*Pi, x = 0 .. 3, style = point, symbol = [box, cross, circle], labeldirections = [horizontal, vertical], labels = ["&theta; (radians)", "Switching  field (T)"], color = [blue, red, black], symbolsize = 13);

Dear sir, 


I request your kind help with regard to the above cited subject.



With warm regards.



Associate Professor in Mathematics.

Aurora's Technological & Research Institute, 


What coding will display a spherical cap of an arbitrary sphere (an arbitrarily located centre and arbitrary radius) where the cap covers an arbitrary solid angle and the radius to the cap's centre is any arbitrary radius of the sphere?

Assume the sphere is defined by its radius and its azimuth and polar angles.

hello experts,

I was using maple for a physical  problem,

and things turned very complicated with a equation with bessel function in it,

like this BesselI(1, (0.9067480359e-2+0.9067480359e-2*I)*sqrt(f))andBesselI(0., (0.9067480359e-2+0.9067480359e-2*I)*sqrt(f)),

which include complex,

the whole equation is as followed:

the variable is f and RV,,dependent variable is RV.

how am I supposed to plot RV when f=100..4000?

it is certain RV has real part and imaginary part,maybe i need a 3Dplot?

please let me know if you have any idea.

best regards,


-7.873519774*10^18*RV^4+(2.676513624*10^12-3.842712573*10^15*(-1)^(3/4)*BesselI(1, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))/(Pi*sqrt(f)*BesselI(0, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))*(4+623.8617593*(-1)^(3/4)*BesselI(1, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))/(BesselI(0, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))*sqrt(f)))))*RV^2+80864.83845+1.440831316*10^9*(-1)^(3/4)*BesselI(1, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))/(Pi*sqrt(f)*BesselI(0, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))*(4+623.8617593*(-1)^(3/4)*BesselI(1, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))/(BesselI(0, 0.1282335370e-1*(-1)^(1/4)*sqrt(f))*sqrt(f)))) = 0

want to plot the relationship between RV and f,but how?


complexplot(RV, f = 100 .. 4000, labels = ["Re", "Im"])

complexplot(RV, f = 100 .. 4000, labels = ["Re", "Im"])





I'm having problems plotting these three functions on the same graph. I can get two of them, but as soon as I try to add the third, it throws a couple errors no matter what I do. It doesn't help that the docs are confusing. The first two functions are approximating diff(y(x),x) = -2xy(x)+1 with the Euler method (step size of 0.05 and 0.1) that I implemented using a for loop. I gather the data into an array and then plot those. That works. But once I try to set y(x) equal to a Taylor approximation, I can plot that graph but the values for the other functions are strangely absent. I've attached my work. What is wrong with my syntax and how can I plot all three functions at once?

I want to plot a function f(x,y)=x^3-2x^2y. But I only want to plot values where $x^2+y^2 <= 1$. So far I specified the plotting range in the standard way, but I really need to specify it like I described. I want to use plot3d.

I am trying to plot some graphs for a differential equations class. I need to plot two equations (soln1a, soln1b1) in one graph, with the two curves overlayed on each other. I am trying to plot them on a graph named "gr1c". The error I get says "Error, (in plots/multiple) empty plot". Clicking on the error results in a webpage saying that "There is no help page available for this error".

When I plot each equation separately, there are no errors. The graphing command that works for soln1a is 



and the graphing command that works for soln1b1 is


My code is as follows. I would greatly appreciate any help ASAP to remove the error. Thank you!



deG:=diff(theta(t),t,t) + mu*diff(theta(t),t)+K*sin(theta(t))= 0;
deL:=diff(theta(t),t,t) + mu*diff(theta(t),t)+K*theta(t)= 0;
Iv:=theta(0)=0.75, D(theta)(0)=2.0;



dear all,

I have a problem when I try to find points on boundary of the surface.

the surface ploted from a matrix as follow:


please help me to generate a bondary curve of this surface.





1 2 Page 1 of 2