Items tagged with plot3d plot3d Tagged Items Feed

Hello everyone!

I'm working with Maple for 8 years and this is the first time I encounter such a problem:

I solve a linear ODE system for my quantum Mechanics research and i want to calculate a quantity formed by the solutions of the system. Everything is going fine with laplace transformation and the symbolic answer is very fast although extremely big (it uses RootOf simplification)

When i use evalf to take numerical approximations of the quantity, the answer is very fast. OK

When i leave unspecified two parameters in order to get a 3d plot the answer takes about 40mins to appear. I think it is too slow but anyway, i can live with it. OK

I want this 3d plot for different parameters of the ODE so i have to do the above process many times. But for some parameter values of the ODE, i get a 3dplot for which:

1) The display of numpoints (surface and line) is incorrect.

2) Plot gets deformed as long as i rotate it!!! The orientation angle changes the appearance of the plot, something obviously unacceptable. I have access to maple 17, 15, 14, and 12. Maple 17 and 15 exhibit this behavior (64 bit Maple editions). But when i tried with Maple 12-standard interface(64 bit) and Maple 14 -Classical interface (32bit) everything went fine, but at the cost of  300 extra secs. I got a solid 3d-plot that is invariant under rotations.

My video card is Ati radeon 6800 HD and i installed the latest drivers. I don't know why is this happening. I thought of putting the blame on me, but in maple 12 and maple 14 everything is fine. It is a pity, because Maple 17 gives very nice and smooth 3d plots.

The first plot is the correct one with maple 12. The second is produced with maple 14. All good.

The third and the second are produced with different orientations of the same plot (Maple 17). Wrong display of data and change of image after rotation. The final one is again from Maple 17 but with style=points. It is in is in full agreement with the first two plots. 

Conclusion: Style=surface and line doesnt produce the correct plot, while style=point does. As for hardware acceleration, i unticked it with no results.

Extra information: When i try to open the maple 12 worksheet which is correct, with Maple 17, again i take the wrong displayed and unstable picture. Moreover only maple 14 classical edition works well. Maple 14 standard interface is also pathological.

Heeeeelp guys!!!



Suppose i have an mx3 matrix M containing a 3D data where 1,2,3 columns are x,y,z coordinates.

Here is what i tried,

Statistics:-SurfacePlot(M[..,1],M[..,2],M[..,3]) does not produce what i want.

plots:-surfdata(M) plots what i want but only works when i use subset of the data for all the data it just hangs.

Statistics:-ScatterPlot3D(M, lowess, fitorder = 2) produces a smoothened version but i want the interpolated plot.

Please help. Here is one sample of the data i want to plot, Mdata.xlsxMdata.xlsx


So if you are given a 3x3 matrix i was just wondering how to use plot3d and plot the 3 column vectors given by it?  Thanks.

So im trying to write a maple script that computes the Jordan form of a given (3x3)- matrix
A. If {a,b,c} is a basis with respect to which A is in Jordan form, then I'm trying to make it
plot the three lines spanned by a, b and c, in the standard coordinate system. I was hinted to use plot3d here.

sidenote: I know how to compute the jordan matrix of A, such by find the eigen vectors and generalised eigen vectors and putting them in as columns in a 3x3 matrix say S,   where S is invertible    then  (S^-1)*(A)*(S) = (J).

Thanks in advance. <3

Hi guys!


I want to do a fractalzoom into the mandelbrot fractal. The basic idea was to produce all single plot pictures, and then display/animate the zoom by iterate through the list. But I cant find a solution that really works. Any help would be very appreciated :)


So where am I? A single picture is made like this: (Sorry math tool didnt take it)


plot3d(0, -2 .. 1, -1.5 .. 1.5, orientation = [-90, 0], grid = [250, 250], style = patchnogrid, color = mandelbrot, lightmodel = None):

As u can see, I use the colorfunc property to color the plane. The mandelbrot procedure evalueates how many iterations are needed to go beyond 2 (> 2).

I cant post pictures, so heres the code if u want to reproduce some:

mandelbrot := proc (x, y)

local c, z, iterationCount, position, iterations;

z := 0; iterations := 0; iterationCount := 35;

position := Complex(x, y);

c := evalf(position);

while iterations < iterationCount do

    if 2 <= abs(z) then return iterations;

    else z := z^2+c; iterations := iterations+1;

    end if


return iterations



As you can see i get a list of PLOT3D structures, each with different axis sizes or views. The problem now is when i use the straightforward way:

display(plots, insequence=true);

the animated frames get smaller and smaller, because the frame is set to the same size as the first frame. 

Do you have any idea how I can make the axis somewhat dynamic, or do you see a better approach?


Here's the rest of the code to produce the effect. Its nearly hardcoded, because its still in experimental phase, sorry for that :P

The test does the following: Start with window -2 .. 1, -1.5 .. 1.5 and with linear interpolation zoom in to the window -0.82 .. -0.7, -0.2 .. -0.08.

 v is used to store the view explicitely, is explained afterwards


steps := 3:

d1 := (2-0.82)/steps:

d2 := (1+0.7)/steps:

d3 := (1.5-0.2)/steps:

d4 := (1.5+0.08)/steps:

l := Array(1 .. 4):

v := Array(1 .. 4):

for i from 0 to steps do

    v(i+1) := [d1*i-2 .. -d2*i+1, -1.5+d3*i .. 1.5-d4*i, 0 .. 0];

    l(i+1) := plot3d(0, d1*i-2 .. -d2*i+1, -1.5+d3*i .. 1.5-d4*i, orientation = [-90, 0], grid = [100, 100], style = patchnogrid, color = mandelbrot, lightmodel = None);



display a single picture:



create the animation:

ll := convert(l, list):

display(ll, insequence = true);


Another approach:

I tried to explicitely give the view as argument:

display(l(i), view = v(i)), v(i) stored in previous render loop. This works, but then the question to resolve is how to draw each frame into the same plot figure and not creating a new plot window for each one?


Now you should have all informations, any help is very welcome!

Thanks in advance,




thank for any  help.

I need you help to make a movies ( animation)  according to the index k=1..10.  k is used in Vect_T[k], vector of time to make an animation of the true solution, and it's also used in Matrix(M[..,..,k]), h

True solution:ow can I make a seq of plots for different value of k. 


plot3d(True_solution(x,y,Vect_T[k]),x=0..a,y=0..b,axes=boxed,grid= [11,11]);

Analytic solution

plots:-matrixplot(Matrix(M[..,..,k]),axes=boxed,tickmarks= [[1="0", 2=" ", 3=" ", 4=" ", 5=".4", 6=" ", 7=" 0.6", 8=" ", 9=".8", 10=" ", 11="1"],[1="0", 2=" ", 3=" ", 4=" ", 5=".4", 6=" ", 7=" 0.6", 8=" ", 9=".8", 10=" ", 11="1"],default],labels=[x,y,u],axes = boxed);

I tried, this:

S:=seq(plots:-matrixplot(Matrix(M[..,..,i]),labels=[x,y,u],axes = boxed),i=1..10):
 But I get only one figure.


I tried

for k from 1 to 10 do
BoxTrue:=BoxTrue, plot3d(True_solution(x,y,Vector_time[k]),x=0..a,y=0..b,axes=boxed,grid= [11,11]): od:
display(BoxTrue, insequence=true);
also, no display


y = -x+sqrt(-5*x)

how to plot this in 2d or 3d graph


TSol:= f(x,y,Q); # The actual funtion is in the attached maple sheet because of its size.

animate( plot3d, [TSol,y=-1..1,x=0..1], Q=0..2,contours=65,style=patchnogrid, grid=[80,80], orientation=[-360,-180], lightmodel=light4, shading=zhue, transparency=0.3,style=contour);



plot3d(subs(Q=2,TSol),y=-1..1,x=0..1,contours=65,style=patchnogrid, grid=[80,80], orientation=[-360,-180], lightmodel=light4, shading=zhue, transparency=0.3,style=contour);


For Q = 2, why there is a big difference between these two ouputs? (If the maple sheet is not accessible then let me know, I will copy past the function (f(x,y,Q))?)

Hello my friends

I have a problem in displaying the mesh grid in my plot3d graphs in maple 17. For instance, I when I type:

plot3d(x^2-y^2, x = -1 .. 1, y = -1 .. 1)

The graph has no mesh grid on it. I mean I am not seeing the black lines on the 3d surface. Can anyone help me to display the mesh grid please? Thanks. 

I'm trying to animate transparency, here's are the plots:

A :=tranz-> plot3d(-2/sqrt(x^2+y^2)+5/sqrt((x-1.6)^2+y^2), x = -5 .. 5, y = -5 .. 5, view = [-5 .. 5, -5 .. 5, -5 .. 5], transparency = tranz) end proc;
B := sphere([0, 0, 0], 2*(1/10), color = magenta, style = patchnogrid):
C := sphere([1.6, 0, 0], 5*(1/10), color = green, style = patchnogrid);
E := plot3d(0, x = -5 .. 5, y = -5 .. 5, view = [-5 .. 5, -5 .. 5, -5 .. 5], style = wireframe, shading = zgrayscale):

here I am displaying them without a problem:

display(A(.5), B, C, E, scaling = constrained, view = [-5 .. 5, -5 .. 5, -5 .. 5], axes = normal)

here I am trying to animate the transparency in plot A to no avail :( :

animate(display, [A(tranz), B, C, E, scaling = constrained, view = [-5 .. 5, -5 .. 5, -5 .. 5], axes = normal], tranz = 0.1 .. 0.9)
Error, (in plot3d) expecting option transparency to be of type {"default", realcons} but received tranz

if you guys could help me find where I'm going wrong I'd be super grateful :D

thanks everyone.



Hi everyone,




plot3d({lambda1,lambda2}, S=-10..10, alpha=-5..5, contours=50,grid=[100,100]);








contourplot({lambda1,lambda2}, S=-10..10, alpha=-5..5, contours=50,grid=[100,100]);



How to minimize the disconnect in the 3d plot and countourplot to get a smooth plot?


Plot3d color oddity...

December 28 2013 Julesp 35

Hello, I was trying to control color of a plot3d. 
I find this answer :
And this post of @Carl Love : 

Here's how to do it with a continuous transformation to your existing color function, which is presumed to return a value between 0 and 1 (the HUE color scale). Keeping it continuous is very very nice when you want colors to represent  numeric values. Let's say your existing color function is C, and your coordinate functions for a parametrized surface are FxFyFz.

Gamma:= 1.15:
     [Fx, Fy, Fz],  a..b, c..d,
     color=  [
          (x,y)-> (1-C(x,y))^Gamma/3, #Hue
          (x,y)-> 1-C(x,y)/4,         #Saturation
          (x,y)-> 1-C(x,y)/7,         #Value
          colortype= HSV
     lightmodel= NONE,
     style= patchnogrid     

There are several parameters that can be adjusted; I've chosen some of them by my personal taste for color .

  • Gamma controls the evenness of the distribution between red and green. I gave this one a name because this is a well-known concept (see the Wikipedia article "Gamma correction").
  • The 3 in the Hue selects the fraction (1/3 in this case) of the full color spectrum that you want. If you want green to red, it will need to be pretty close to 3.
  • The Hue value is subtracted from 1 to make the scale go green to red rather than red to green.
  • The 4 in the Saturation controls (to some extent) how "light" the light-green is.
  • The 7 in the Value controls (to some extent) how dark the dark-red is (lower values will make it darker).
  • lightmodel= NONE is used so that the colors will not change due to shadows when the plot is rotated. "

I made some test to see the impact of the Gamma parameter. 
And with Gamma = 1, it's odd. 



C := proc (x, y) x end proc;

proc (x, y) x end proc


Gamma := 1.15:


a := 1:



It looks like with gamma = 1, plot3d makes an automatic scaling of the colors.
But I don't understand why.
Does anyone know ?


How to draw an object in 3d which is like a spring

3 layers of circle already enough

Is there some sort of analog of  option  discont=true  in  plot3d?  Can I remove these vertical cylindrical surfaces in the following example? They correspond to the zeros of the denominator:

plot3d((x^2+sin(x))/(y^2+cos(x)), x = -7 .. 7, y = -5 .. 5, view = [-7 .. 7, -6 .. 6, -1 .. 7], numpoints = 50000, scaling = constrained, axes = normal);

plots[implicitplot](y^2+cos(x), x = -7 .. 7, y = -5 .. 5, scaling = constrained, numpoints = 5000);



how to draw this...

November 06 2013 rit 200


n := 31;

h := 0;

c1 := Array(1 .. 100);

c2 := Array(1 .. 100);

c3 := Array(1 .. 100);

for tt from 0 to n-4 do

c1[tt+1] := ...

c2[tt+1] := ...

c3[tt+1] := ...


plot3d([c1[x],c2[x],c3[x]], x=1..27, y=1..27);

1 2 3 4 5 Page 1 of 5