summaryrefslogtreecommitdiff
path: root/graphics/circuit_macros/examples/svg/LoglogSVG.m4
blob: d6b29d0fa0297e8e366a19b675687e0e060d831a (plain)
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
.PS
# LoglogSVG.m4
gen_init(svg_font(sans-serif,11bp__))

define(`rint',`sign(`$1')*int(abs(`$1')+.5)')
S:[
  horiz = 4
  vert = horiz*2/3
  tmax = 10
  tmin = 1
  vmax = 100
  vmin = 0.1
  tic = 0.08
  denv = log(vmax) - log(vmin)
  denh = log(tmax) - log(tmin)

  Origin: Here
  define(`coord',
    `Origin+((log(`$1')-log(tmin))/denh*horiz,(log(`$2')-log(vmin))/denv*vert)')

thinlines_
# Left axes and labels
  psset_(linecolor=gray)
  for i = rint(log(vmin)) to rint(log(vmax)) do {
    move to coord(tmin,exp(i))
    line right horiz
    sprintf("%g",exp(i)) wid 0.2 rjust at last line.start
    if i < rint(log(vmax)) then {
      for j = 2 to 9 do { move to coord(tmin,j*exp(i))
        line right horiz
        if (j==2) || (j==5) then {
          sprintf("%g",j*exp(i)) rjust at last line.start }
        } } }

#Horizontal axis and labels
  for t = rint(log(tmin)) to rint(log(tmax)) do {
    move to coord(exp(t),vmin)
    line up vert
    sprintf("%g",exp(t)) wid 0.15 below at last line.start
    if t < rint(log(tmax)) then {
      for j = 2 to 9 do { move to coord(j*exp(t),vmin)
        line up vert
        if (j==2) || (j==5) then {
          sprintf("%g",j*exp(t)) below at last line.start }
        } } }
  psset_(linecolor=black)

# arrow from Origin up vert chop 0 chop -0.1
  line from Origin up vert
  "svg_it(y(t))" wid 0.3 at Origin+(-0.1,0.5*vert) rjust

# arrow from Origin right horiz chop 0 chop -0.1
  line from Origin right horiz
  "svg_it(t)" at Origin+(.5*horiz,-0.1) below

thicklines_
  alpha = 2.5
  beta = 3.0
  tm = max(tmin,(vmin/alpha)^(1/beta))
  tM = min(tmax,(vmax/alpha)^(1/beta))
  line from coord(tm,alpha*tm^beta) to coord(tM,alpha*tM^beta)
P: box invis fill_(1) wid 0.75 ht 0.18 \
   "svg_it(y(t) = 1.5 t)" with .w at coord(2.005,14)
   "svg_it(svg_small(3))" with .w at P.ne #ljust

  alpha = 0.05
  beta = 2.0
  tm = max(tmin,(vmin/alpha)^(1/beta))
  tM = min(tmax,(vmax/alpha)^(1/beta))
  line from coord(tm,alpha*tm^beta) to coord(tM,alpha*tM^beta)
Q: box invis fill_(1) wid 0.80 ht 0.18 \
   "svg_it(y(t) = 0.05 t)" with .e at coord(4.6,1.4)
   "svg_it(svg_small(2))" with .w at Q.ne #ljust

]

# box wid S.wid ht S.ht at S
# print (S.wid,S.ht)/(1 pc__)

 command "</g>" # end font
.PE