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
|
clear
window 0
defaults
!
generate theta 0,,180 100 ! make a vector with 100 elements from 0 to 180
x = cosd(theta)
y = x^2
x2 = [0;45;90;135;180] ! the upper axis large tic mark locations
x3 = [5;10;15;20;25;30;35;40] ! the upper axis small tic mark increments
lx2 = len(x2) ! length of vector x2 = number of large tic marks
y2 = [1:lx2] ! y2 = [1;2;3;4;...]
xupper = cosd(x2) ! convert to functional values
lx3 = len(x3) !
y3 = [1:lx3] !
set
%yupperaxis 85 ! make room for a label
%xloweraxis 15
graphbox 0 ! turn off automatic tic marks on top axis
plotsymbol 0 ! no plotting symbol
autoscale on ! autoscale the plot
xlabel 'Cos(<theta>)'
ylabel 'Cos<^>2<_>(<theta>)'
xlabelon 1
ylabelon 1
curvecolor red
ydroplastnumber 1 ! do not draw ymax
graph x y ! draw the graph with the bottom x-axis
set
autoscale off
curvecolor black
xticson 0
graphbox 1
graph\axesonly
!
get
%xnumbersheight xnumsz ! number size on x-axis
%yloweraxis ylaxis ! lower end of y-axis (as % of window)
%yupperaxis yuaxis ! upper end of y-axis (as % of window)
%ximagticlength xiticl ! distance from x-axis to top of numbers (as %)
yloweraxis ylaxisw ! lower end of y-axis (world units)
yupperaxis yuaxisw ! upper end of y-axis (world units)
ximagticlength xiticlw ! distance from x-axis to top of numbers (world units)
ymin ymin ! number at bottom of y-axis
ymax ymax ! number at top of y-axis
ticlen = xiticlw/(yuaxisw-ylaxisw)*(ymax-ymin)
y1 = ymax-2*ticlen
y11 = ymax-ticlen
yloc = yuaxis+1
set
%textheight xnumsz ! height of text (% of window height)
textalign 2 ! centre justify text
textinteractive 0 !
%ytextlocation yloc ! vertical location of text
textangle 0 ! text angle (degrees)
!set clip 0
world\percent xupper y2 xp yp ! convert graph units to % of window
do j = [1:lx2] ! for each large tic mark ...
set %xtextlocation xp[j] ! set horizontal location of text
text rchar(x2[j]) ! draw the upper x-axis numbers
graph\overlay [xupper[j];xupper[j]] [ymax;y1] ! draw large tic marks
xtmp = cosd(x2[j]+x3) ! on top of box
world\percent xtmp y3 xpp ypp
do k = [1:lx3]
graph\overlay [xtmp[k];xtmp[k]] [ymax;y11] ! draw small tic marks
enddo ! on top of box
enddo
! find the mid point of the x-axis
get
%xupperaxis xux
%xloweraxis xlx
set
%xtextlocation (xux+xlx)/2 ! position a label at the mid point
%ytextlocation 90
textalign 2 ! centre justify
textinteractive 0
text '<theta><degree>'
|
