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
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
|
proc x01 {{w loopback}} {
global xscale yscale xoff yoff
# For starting from scratch this call to pladv increments cursub, but
# then the following plssub sets it to zero so the whole thing is
# essentially a nop. However, for the case when other examples are run
# first, this call to pladv is absolutely essential to finish the
# preceding page.
$w cmd pladv 0
$w cmd plssub 2 2
# First plot
set xscale 6.0
set yscale 1.0
set xoff 0.0
set yoff 0.0
plot1 $w
# Second
set xscale 1.0
set yscale 0.0014
set yoff 0.0185
$w cmd plsyax 5
plot1 $w
# Third
plot2 $w
# Fourth
plot3 $w
# Restore defaults
# $w cmd plcol0 1
$w cmd plssub 1 1
$w cmd pleop
}
# This is supposed to work just like the plot1() in x01c.c
proc plot1 {w} {
global xscale yscale xoff yoff
set npts 60
matrix x f $npts
matrix y f $npts
for {set i 0} {$i < $npts} {incr i} {
x $i = [expr $xoff + ($xscale * ($i + 1)) / $npts]
y $i = [expr $yoff + $yscale * pow([x $i],2)]
}
set xmin [x [expr 0]]
set xmax [x [expr $npts-1]]
set ymin [y [expr 0]]
set ymax [y [expr $npts-1]]
matrix x1 f 6
matrix y1 f 6
for {set i 0} {$i < 6} {incr i} {
set j [expr $i*10+3]
x1 $i = [x $j]
y1 $i = [y $j]
}
$w cmd plcol0 1
$w cmd plenv $xmin $xmax $ymin $ymax 0 0
$w cmd plcol0 2
$w cmd pllab "(x)" "(y)" "#frPLplot Example 1 - y=x#u2"
# plot the data points
$w cmd plcol0 4
$w cmd plpoin x1 y1 9
# draw the line through the data
$w cmd plcol0 3
$w cmd plline x y
}
# This is supposed to work just like the plot2() in x01c.c
proc plot2 {w} {
$w cmd plcol0 1
$w cmd plenv -2 10 -.4 1.2 0 1
$w cmd plcol0 2
$w cmd pllab "(x)" "sin(x)/x" "#frPLplot Example 1 - Sinc Function"
# Fill up the array
matrix x1 f 100
matrix y1 f 100
for {set i 0} {$i < 100} {incr i} {
set x [expr ($i - 19.)/6.]
x1 $i = $x
y1 $i = 1
if {$x != 0} { y1 $i = [expr sin($x)/$x] }
}
$w cmd plcol0 3
$w cmd plwidth 2
$w cmd plline x1 y1
$w cmd plwidth 1
}
# This is supposed to work just like the plot3() in x01c.c
proc plot3 {w} {
$w cmd pladv 0
$w cmd plvsta
$w cmd plwind 0.0 360.0 -1.2 1.2
# Draw a box with ticks spaced 60 degrees apart in X, and 0.2 in Y.
$w cmd plcol0 1
$w cmd plbox "bcnst" 60.0 2 "bcnstv" 0.2 2
# Superimpose a dashed line grid, with 1.5 mm marks and spaces.
# plstyl expects two integer matrices for mark and space!
matrix mark i 1
matrix space i 1
matrix none i 0
mark 0 = 1500
space 0 = 1500
$w cmd plstyl mark space
$w cmd plcol0 2
$w cmd plbox "g" 30.0 0 "g" 0.2 0
mark 0 = 0
space 0 = 0
$w cmd plstyl none none
$w cmd plcol0 3
$w cmd pllab "Angle (degrees)" "sine" "#frPLplot Example 1 - Sine function"
matrix x f 101
matrix y f 101
for {set i 0} {$i < 101} {incr i} {
x $i = [expr 3.6 * $i]
y $i = [expr sin([x $i] * $::PLPLOT::PL_PI / 180.0)]
}
$w cmd plcol0 4
$w cmd plline x y
}
|
