Function: quickplot3d - fast 3-d curve or point plotting of a list or hfarray

Calling sequence:

quickplot3d(L,options);

Parameters:

L - a list, hfarray or Matrix of the points to be plotted. In the case of a list, each element is a point (represented as a list with three elements). An hfarray or Matrix should have two indices: the first goes from 1 to the number of points, the second goes from 1 to 3. In the case of a Matrix the datatype should be float[8] and the order should be C_order .

options - most of the usual three-dimensional plot options.

Description:

Examples:

> H:= hfarray(1..3000,1..3,
(i,j)->if j=1 then 0.001*i elif j=3 then cos(.1*i) else sin(.1*i) fi);

                          [ 1..3000 x 1..3 2-D Array ]

                     H := [ Data Type: float[8]      ]

                          [ Storage: rectangular     ]

                          [ Order: C_order           ]

> quickplot3d(H,style=line,axes=box,scaling=constrained);

> quickplot3d(H,axes=box,coords=spherical,scaling=constrained);

The next example represents the attractor of the Rössler system of differential equations for , and .

> de:= {diff(x(t),t) = -y(t)-z(t),diff(y(t),t)=x(t)+y(t)/5,
diff(z(t),t)=1/5+z(t)*(x(t)-5),x(0)=0,y(0)=0,z(0)=0}:

> sol:= dsolve(de,[x(t),y(t),z(t)],numeric,value=
array([seq(100+0.05*j,j=0..4000)]),maxfun=50000):

> pts:= hfarray(linalg[submatrix](sol[2,1],1..4001,2..4));

                           [ 1..4001 x 1..3 2-D Array ]

                    pts := [ Data Type: float[8]      ]

                           [ Storage: rectangular     ]

                           [ Order: C_order           ]

> quickplot3d(pts,axes=frame,title="Rössler system",shading=XY,labels=[x,y,z]);

See also:

Maple Advisor Database R. Israel 2000