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
|
#
# dreieck.m
#
# (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule
#
w = 12
N = 10
rand("seed", 1);
angles = 80 * rand(1,N)
radii = 2 * rand(1,N) + 0.4
angle = 20
radii = radii * w / (cosd(angle) * sum(radii))
radius = sum(radii)
radius * cosd(angle)
points = zeros(2,N);
ray = zeros(2,N);
p = [ 0; 0 ];
for i = (1:N)
p = p + radii(1,i) * [ cosd(angles(1,i)); sind(angles(1,i)) ];
points(:, i) = p;
ray(:, i) = sum(radii(1,1:i)) * [ cosd(angle); sind(angle) ];
end
points
ray
fn = fopen("drei.inc", "w");
for i = (1:N)
fprintf(fn, "\\coordinate (A%d) at (%.4f,%.4f);\n", i,
points(1,i), points(2,i));
fprintf(fn, "\\coordinate (B%d) at (%.4f,%.4f);\n", i,
ray(1,i), ray(2,i));
end
fprintf(fn, "\\def\\r{%.4f}\n", radius);
fclose(fn);
|
