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           ]```