Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
#! /usr/bin/perl
$NSTEPS = 32;
$RMIN = 0.1;
$RMAX = 0.5;
$NDOTS = 15;
$THETAMAX = 3.14159;
printf "#! /usr/bin/env partiview
# Generated by mksnake.pl
datavar 0 lum
datavar 1 color
eval lum lum %d 0
eval color color
eval slum 1000
eval speed 8
", $NDOTS+2;
for($time = 0; $time <= $NSTEPS; $time++) {
$tfrac = $time / $NSTEPS;
printf "\n";
printf "datatime %d\n", $time;
for($i = 0; $i < $NDOTS; $i++) {
$r = $RMIN + ($RMAX-$RMIN) * $i/$NDOTS;
$theta = $tfrac * $THETAMAX * sqrt($RMAX / $r);
printf "%g %g %g %d %g\n", $r*cos($theta), $r*sin($theta), 0, $i, $theta;
}
}