Items tagged with plot3d plot3d Tagged Items Feed

In a 3d coordinate system I have a circular spacecurve with z-minimum -4 and z-maximum +4. In the same 3d coordinate system I have a 3d surface plot with z-minimum -0.5 and z-maximum +1.3 . When I choose the color option "Z(Hue)" in order to color-code the z-values on the 3d surface and make the topography more clear, I mostly get a totally green 3d surface. It seems that the color scaling is coupled with the spacecurve with z-values of +-4 . How can I uncouple the color scaling from the spacecurve and couple it with the z-range of the surface, while the color-limits shall be at +-1.3 ?

Hi everyone,

 

I am using the following code (by dr. Corless) to animate a Riemann surface:

 

restart;
with(plots);
B := 1.5;
u2 := r*cos(th); v2 := r*sin(th);
w1 := u1+I*v1; w2 := u2+I*v2;
z1 := evalc(w1^2); z2 := evalc(w2^2);
x1 := evalc(Re(z1)); x2 := evalc(Re(z2));
y1 := evalc(Im(z1)); y2 := evalc(Im(z2));
f1 := proc (theta) options operator, arrow; plot3d([-x1, -y1, v1], u1 = -6 .. 1, v1 = -B .. B, grid = [50, 50], orientation = [theta, 80], color = u1) end proc;
f2 := proc (theta) options operator, arrow; plot3d([-x2, -y2, v2], r = 1 .. 1, th = -Pi .. Pi, grid = [50, 50], orientation = [theta, 80], color = black) end proc;
display(seq([f1(10*k), f2(10*k)], k = -17 .. 18), insequence = true, axes = box, view = [-1 .. 1, -1 .. 1, -B .. B]);

f1 is the Riemann surface and f2 is the winding curve.

The animation works fine on the individual plots, i.e. when I do:

display(seq(f1(10*k), k = -17 .. 18), insequence = true, axes = box, view = [-1 .. 1, -1 .. 1, -B .. B]);

or

display(seq(f2(10*k), k = -17 .. 18), insequence = true, axes = box, view = [-1 .. 1, -1 .. 1, -B .. B]);

but with both f's (as in my modification)  the animation does not combine correctly the two plots and shows them separately. Is there any way I can combine the plots so that display produces a smooth animation with each rotated frame containing its winding curve?

 

Many thanks,

Yiannis

The regular plot3d command takes in the range of x and y but if I have been asked to "choose a range of z to get a nice picture" what command can I use apart from implicitplot3d, which in this case doesn't seem appropriate?

For my research at my school, I need to plot the Schrodinger equation using my given wave function. I am having issues though as I am receiving an error telling me that I need to submit my argument in standard or parametric form.

 

restart;

schro := {[4*(x^2-1)*(diff(psi(x, y), x, x))/(x^2-y^2)+2*x*(diff(psi(x, y), x))-(y^2-1)*(diff(psi(x, y), y, y))-2*y*(diff(psi(x, y), y))+(-2/(x+y)-2/(x-y)+1)*psi(x, y)] = energy};

psi := proc (x, y) options operator, arrow; e^((-1)*.5*[2*x+2*y])+e^((-1)*.5*[2*x-2*y]);

plot3d(schro, x= 1..100, y = -1 ..1);

 

Unfortunately, I am unable to plot my function. I was able to substitute my psi(x,y) function into my schro equation, but I am unable to do much else. I want to plot my function such that my energy is represented in the z-axis as a relationship between x and y which are a specific coordinate system (conical coordintes). 

Hi,

How do you plot a list of values/an array into a surface?

E.g.

t = [seq(0+i*(2*evalf(Pi)*(1/10)), i = 0 .. 10)]

x:=[3,4,2*evalf(Pi), 7.83]:

y:=[2.5,4.3,6,2*evalf(Pi)+2]:

z:=[-2,0.3,1.5,evalf(Pi)]:

plot3d([x[i]*cos(t), x[i]*sin(t), y[i]*cos(t-z[i])],-10..10,-10..10);

 

 

I would like to look at a 3d plot including an condition about the two variables of the plot

For this simple case, I only want to see the plot with the condition that x>0. Is that possible?

plot3d(x*y,[x,0,1],[y,0,1]);

Hi,

 

For my thesis I would like to illustrate the inclination of the solar system objects in a plane. Imagine the Solar system as a circular plane which is the average of the motions of all objects orbiting the Sun. Each individual planet/object is inclined towards this averaged plane - some more, some less. E.g. this image. Behind Neptune there lies the so-called Kuiper belt with many thousands of dwarf planets (Pluto is one of them and there are over 1000 objects known out there already). Now imagine that the big heavy Jupiter and other big planets perturb those small objects out there: so they are also inclined towards the average plane, see this image. Their inclination depends on their radial distance to the Sun (measured usually in Astronomical Units...1 unit is the distance Sun-Earth...we are talking about 40 to 50 units here). For one, I want to visualise this: imagine a circular plane and each orbit out there has another angle to the average plane. This is the first. But now: imagine the circle with it's 360°. Each object reaches it's highest point on its orbit around the Sun on another angle on this 360° circle. Neptune e.b. at 170°, Pluto at 250° etc. You get the picture. So not only are the objects in the Kuiper belt inclined differently, but their maximum orbit positions are also scattered across an imaginary 360° circle. I want to show this with a 3d-plane like this image but not with two peaks for one orbit. Is there a way to do/plot/visualise this in Maple? I am just interested in a visualisation of the principle without any empirical data behind this.

Thanks.

I am trying to make a plot3d on my Surface Pro  3 running Windows8.1.

The plot looks like a ghost.  Pale grey.  It looked nice and colourful when I executed the same command on an older computer a few hours ago.  (I know it was exactly the same command because I saved the file on my old computer and then opened in on the Surface Pro 3)

Is there a way to fix this?

 

plot3d(-x^2+y^2-3*x+2*y+100, x = -5 .. 8, y = -5 .. 5, transparency = 0, lightmodel = light1, glossiness = 1, filled = true, style = wireframeopaque, shading = zhue, view = [-5 .. 10, -5 .. 8, 0 .. 150]);

 OOPS!  Now the problem has corrected itself.  

I'm getting worried about all the weird things Maple is doing on this Surface Pro 3.  I'm trying to make up a final exam tonight (the exam is tomorrow).

As stated, contourplot3d is not displaying contours. Even in the Maple Help examples the contours are not displayed in the help screen.

In this example:

contourplot3d(-5*x/(x^2+y^2+1), x = -3 .. 3, y = -3 .. 3, filledregions = true);

which is directly from Maple Help, the surface is plotted (with boxed axes), but no contours appear.

When I try

contourplot3d(-5*x/(x^2+y^2+1), x = -3 .. 3, y = -3 .. 3);

the plot shows nothing except boxed axes.

I've also tried

plot3d(-5*x/(x^2+y^2+1), x = -3 .. 3, y = -3 .. 3, contours=6);

and the surface shows but without countours. Furthermore, the Maple scripts that worked previously in earlier versions (showing contours via plot3d) now do not show contours in Maple 18. This occurs on both Win 7 and Win 8 laptops.

Thanks!

 

I'm importing a STL geometry with plottools[importplot].

The result is a PLOT3D data structure. Applying "op" to the structure reveals that it is formed by POLYGONS structures.

How can I import the same STL geometry, and, instead of getting a POLYGONS structure, getting a POINTS structure?

I think it should be direct for Maple, because in the plot view configs one can select the style to be surface, contour and even point.

If the conversion has to be by hand, does someone know a not too painful method? For me it's being painful, as I'm trying a solution all the day.

So this is the errorcode I get and my code is

> r := (x1, y1, x2, y2) options operator, arrow; (sqrt((x2-x1)^2+(y2-y1)^2))^2
> hreg := (A, C, x) options operator, arrow; A*x+C
> hpconf := (Q, x1, y1, x2, y2) options operator, arrow; Q*(1/4)*ln(2.25*b*K*t/(S*r(x1, y1, x2, y2)))/(3.14*K*b)


> he1 := (x, y, C) options operator, arrow; hreg(A1, C, x)+hpconf(Q1, x, y, xp1, yp1)+hpconf(Q2, x, y, xp2, yp2)+hpconf(Q3, x, y, xp3, yp3)+hpconf(Q4, x, y, xp4, yp4)
> C1 := solve(he1(-75, 0, C) = 20, C);


> plot3d(he1(x, y, C1), x = -75 .. 75, y = -75 .. 75, style = surfacecontour, axes = box);

I have defined all the constants, and squareboxes shouldn't be needed for the plot, right?

 

I currently have:

plot3d(78.9*x+49.1*y+15.7*(1-x-y), x = 0 .. 1, y = 0 .. 1)

I have ratios x, y and z and the ratio of z is implicit from the ratios of x and y (z=1-x-y), therefore whenever x+y>1 the surface is meaningless as it gives z as a negative ratio.

I want the surface to be plotted only for x+y<1, so something like:

plot3d(78.9*x+49.1*y+15.7*(1-x-y), x = 0 .. 1, y = 0 .. 1, x+y<1)

But I can't find a suitable command.

Thanks in advance,

Tom.

 

I would like to plot this :

yz(x+y+2z)^3-8=0 such as x^2+x(y+2z)-1<0 and y^2+y(x+2z)-3<0

The difficulty is that bounds depend on other variables. I was trying to use Maple implicitplot3d function, but I don't see a way.

Is it possible to use Maple for such plot ? If so, how ?

 

Thank you for your answer

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!!!

 

 

1 2 3 4 5 6 Page 1 of 6