Download (11Mb) - tuprints

A. Quellcode des verwendeten Skriptes
A.2. Quellcode
1 #################################################
2 #
3 # f u n c t i o n INTERPOLATE
4 #
5 # Do a l i n e a r i n t e r p o l a t i o n between two 2D p o i n t s
6 #
7 #################################################
8
9 function INTERPOLATE( x1 , y1 , x2 , y2 , x3 , y3 )
10 {
11 y3 = y1 + ( ( ( y2−y1 ) / ( x2−x1 ) ) ∗ ( x3−x1 ) ) ;
12 return y3 ;
13 }
14
15 ### END f u n c t i o n INTERPOLATE
16
17 ### BEGIN main program
18
19 {
20 # At f i r s t , f i l l up " readahead " b u f f e r
21 i f ( $1 < s t e p ) {
22
23 # Skip t e x t l i n e s at t h e b e g i n n i n g
24 while ( ! ( match( $1 , " ^[0 −9] " ) ) ) {
25 getline ;
26 }
27 i = 0 ;
28
29 # F i l l t h e b u f f e r , t a k e care t h a t i t w i l l contain one l i n e
30 # over s t e p width
31 do {
32 r [ i ] = $1 ;
33 pot [ i ] = $2 ;
34 i ++;
35 getline ;
36 } while ( r [ i −1] < r [ 1 ] + s t e p )
37
38 # s t o r e s i z e o f b u f f e r f o r l a t e r use
39 # ( Step width o f raw data i s assumed to be c o n s t a n t over
40 # t h e whole i n p u t data . )
41 b u f f e r s i z e = i ;
42
43 } ### END b u f f e r f i l l i n g
44
45 # Get second v a l u e ( ' current ' r−v a l u r e + s t e p width ) from
46 # readahead b u f f e r , c a l c u l a t e s t e p v o l t a g e and p r i n t i t
47 i = 0 ;
174