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
|
############################################################################
#
# File: complex.icn
#
# Subject: Procedures to perform complex arithmetic
#
# Author: Ralph E. Griswold
#
# Date: June 21, 2000
#
############################################################################
#
# This file is in the public domain.
#
############################################################################
#
# The following procedures perform operations on complex numbers.
#
# complex(r,i) create complex number with real part r and
# imaginary part i
#
# cpxadd(z1, z2) add complex numbers z1 and z2
#
# cpxdiv(z1, z2) divide complex number z1 by complex number z2
#
# cpxmul(z1, z2) multiply complex number z1 by complex number z2
#
# cpxsub(z1, z2) subtract complex number z2 from complex number z1
#
# cpxstr(z) convert complex number z to string representation
#
# strcpx(s) convert string representation s of complex
# number to complex number
#
############################################################################
record complex(rpart, ipart)
procedure strcpx(s)
s ? {
="(" | fail
return complex(numeric(upto('+-')),
2(move(1), numeric(upto(')')), tab(-1)))
}
end
procedure cpxstr(z)
if z.ipart < 0 then return "(" || z.rpart || z.ipart || "i)"
else return "(" || z.rpart || "+" || z.ipart || "i)"
end
procedure cpxadd(z1, z2)
return complex(z1.rpart + z2.rpart, z1.ipart + z2.ipart)
end
procedure cpxsub(z1, z2)
return complex(z1.rpart - z2.rpart, z1.ipart - z2.ipart)
end
procedure cpxmul(z1, z2)
return complex(z1.rpart * z2.rpart - z1.ipart * z2.ipart,
z1.rpart * z2.ipart + z1.ipart * z2.rpart)
end
procedure cpxdiv(z1, z2)
local denom
denom := z2.rpart ^ 2 + z2.ipart ^ 2
return complex((z1.rpart * z2.rpart + z1.ipart * z2.ipart) / denom,
(z1.ipart * z2.rpart - z1.rpart * z2.ipart) / denom)
end
procedure cpxconj(z)
return complex(z.rpart, -z.ipart)
end
procedure cpxabs(z)
return sqrt(z.rpart ^ 2 + z.ipart ^ 2)
end
|