-
Notifications
You must be signed in to change notification settings - Fork 0
/
2.11.scm
52 lines (50 loc) · 2.02 KB
/
2.11.scm
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
;; (define (mul-interval x y)
;; (let ((p1 (* (lower-bound x) (lower-bound y)))
;; (p2 (* (lower-bound x) (upper-bound y)))
;; (p3 (* (upper-bound x) (lower-bound y)))
;; (p4 (* (upper-bound x) (upper-bound y))))
;; (make-interval (min p1 p2 p3 p4)
;; (max p1 p2 p3 p4))))
(define (sign lower upper)
(cond ((and (< lower 0) (< upper 0)) -1)
((and (< lower 0) (> upper 0)) 0)
(else 1)))
(define (mul-interval x y)
(let ((xl (lower-bound x))
(xu (upper-bound x))
(yl (lower-bound y))
(yu (upper-bound y)))
(let ((xs (sign xl xu))
(ys (sign yl yu)))
(cond ((< xs 0)
(cond ((< ys 0) ; - - - -
(make-interval (* xu yu)
(* xl yl)))
((> ys 0) ; - - + +
(make-interval (* xl yu)
(* xu yl)))
(else ; - - - +
(make-interval (* xl yu)
(* xl yl)))))
((> xs 0)
(cond ((< ys 0) ; + + - -
(make-interval (* xu yl)
(* xl yu)))
((> ys 0) ; + + + +
(make-interval (* xl yl)
(* xu yu)))
(else ; + + - +
(make-interval (* xu yl)
(* xu yu)))))
(else
(cond ((< ys 0) ; - + - -
(make-interval (* xu yl)
(* xl yl)))
((> ys 0) ; - + + +
(make-interval (* xl yu)
(* xu yu)))
(else ; - + - +
(make-interval (min (* xl yu)
(* xu yl))
(max (* xl yl)
(* xu yu))))))))))