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
|
# ex_triangle.ppl
#
# The code in this file is part of Pyxplot
# <http://www.pyxplot.org.uk>
#
# Copyright (C) 2006-2012 Dominic Ford <coders@pyxplot.org.uk>
# 2008-2012 Ross Church
#
# $Id: ex_triangle.ppl 1261 2012-07-11 21:38:05Z dcf21 $
#
# Pyxplot is free software; you can redistribute it and/or modify it under the
# terms of the GNU General Public License as published by the Free Software
# Foundation; either version 2 of the License, or (at your option) any later
# version.
#
# You should have received a copy of the GNU General Public License along with
# Pyxplot; if not, write to the Free Software Foundation, Inc., 51 Franklin
# Street, Fifth Floor, Boston, MA 02110-1301, USA
# ----------------------------------------------------------------------------
# This script produces a figure for the Pyxplot Users' Guide
# Call common initialisation script
reset
title = "ex_triangle"
load "examples/fig_init.ppl"
# BEGIN
set unit angle nodimensionless
set unit of length cm # Display lengths in cm
set unit of angle degree # Display angles in degrees
set numeric sigfig 3 display latex # Correct to 3 significant figure
cm = unit(cm) # Shorthand to save space
deg = unit(deg)
turn(a) = matrix( [cos(a),-sin(a)], \
[sin(a), cos(a)] )
# Define subroutine for drawing triangles
subroutine triangleDraw(B,AB,AC,BC)
{
# Use cosine rule to find interior angles
ABC = acos((AB**2 + BC**2 - AC**2) / (2*AB*BC))
BCA = acos((BC**2 + AC**2 - AB**2) / (2*BC*AC))
CAB = acos((AC**2 + AB**2 - BC**2) / (2*AC*AB))
# Positions of three corners of triangle
C = B + vector(BC,0*cm)
A = B + turn(ABC)*vector(AB,0*cm)
# Draw triangle
polygon [A,B,C]
# Draw angle symbols
arcRad = 0.4*cm # Radius of angle arcs
arc at B radius arcRad from 90*deg-ABC to 90*deg
arc at C radius arcRad from -90*deg to -90*deg+BCA
arc at A radius arcRad from 90*deg+BCA to 270*deg-ABC
# Label lengths of sides
textGap = 0.1*cm
text "%s"%(BC) at (B+C)/2 gap textGap hal c val t
text "%s"%(AB) at (A+B)/2 gap textGap rot ABC hal c val b
text "%s"%(AC) at (A+C)/2 gap textGap rot -BCA hal c val b
# Label angles
arcRad2 = vector(1.4*arcRad , 0*cm)
text "%s"%CAB at A+turn(-90*deg+ABC-BCA)*arcRad2 hal c val t
text "%s"%ABC at B+turn( ABC/2)*arcRad2 hal l val c
text "%s"%BCA at C+turn(180*deg - BCA/2)*arcRad2 hal r val c
# Label points ABC
text "A" at A gap textGap hal c val b
text "B" at B gap textGap hal r val c
text "C" at C gap textGap hal l val c
}
# Display diagram with three triangles
set multiplot ; set nodisplay
call triangleDraw(vector([2.8,3.2])*cm, 3*cm, 4*cm, 4*cm)
call triangleDraw(vector([0.0,0.0])*cm, 3*cm, 4*cm, 5*cm)
call triangleDraw(vector([6.5,0.0])*cm, 3*cm, 3*cm, 3*cm)
set display ; refresh
# END
# Call common cleanup script
load "examples/fig_end.ppl"
|
