(include "common-stalingrad") (letrec ((loop (lambda (i r) (if (zero? i) r (loop (- i 1) (let* ((start (list (real 1) (real 1))) (f (lambda (x1 y1 x2 y2) (- (+ (sqr x1) (sqr y1)) (+ (sqr x2) (sqr y2))))) ((list x1* y1*) (multivariate-argmin-F (lambda ((list x1 y1)) (multivariate-max-R (lambda ((list x2 y2)) (f x1 y1 x2 y2)) start)) start)) ((list x2* y2*) (multivariate-argmax-R (lambda ((list x2 y2)) (f x1* y1* x2 y2)) start))) (list (list (write-real x1*) (write-real y1*)) (list (write-real x2*) (write-real y2*))))))))) (loop (real 1000) (list (list (real 0) (real 0)) (list (real 0) (real 0)))))