gauche で動作する超循環評価器 問題を解く前に、先に超循環評価器を作ってプログラムを実行できる環境を作ってみることにする。 最初、以下のようなエラーが出て、動作しなくて悩んでいたが。 ;;; M-Eval input: (cons ‘1 ‘2) *** ERROR: invalid application: ((primitive #<subr cons>) 1 2) eval …続きを読む

問題4.1 まず、gauche でそれぞれの評価順の場合の手続きを作って確かめてみる。 (define val 10) (define expression ‘((set! val (+ val 2)) (set! val (* val 2)))) 左から右に評価する場合 (define (list-of-values-left-to-right exps) (if (null? exps) ‘(…続きを読む

問題3.82 (define (estimate-integral p x1 x2 y1 y2) (define monte-carlo-stream (stream-map (lambda (m) (* (- x2 x1) (- y2 y1) m)) (monte-carlo random-in-range-experiment-stream 0.0 0.0))) (define random-…続きを読む

問題3.81 問題文が解りにくいが、 generate または reset の要求を表す値（’generate と 整数 を使う）から成るストリームを引数にとって、それぞれの要求に対応して作られた乱数からなるストリームを作る。 (define (rand-update x) (remainder (+ x 1812433253) 4294967296)) (define (rand input-s…続きを読む

問題3.80 (define (RLC R L C dt) (define (rlc vC0 iL0) (define dvC (scale-stream iL (/ -1 C))) (define diL (add-streams (scale-stream vC (/ 1 L)) (scale-stream iL (- (/ R L))))) (define iL (integral (del…続きを読む

問題3.79 (define (solve-2nd f dt y0 dy0) (define ddy (stream-map f dy y)) (define dy (integral (delay ddy) dy0 dt)) (define y (integral (delay dy) y0 dt)) y) 計算機プログラムの構造と解釈 posted with amazlet at 08.11.…続きを読む

問題3.78 図3.35 をそのまま手続きに書き起こす。 (define (solve-2nd a b dt y0 dy0) (define ddy (add-streams (scale-stream dy a) (scale-stream y b))) (define dy (integral (delay ddy) dy0 dt)) (define y (integral (delay dy…続きを読む

問題3.77 被積分値を先に評価する force の部分はすぐにわかったのだけれど、 integral の再帰部分の第1引数を delay しなければいけないところがわからなかった。 solve 手続きは、本文に記載されているものは Gauche では動作しないので修正が必要となった。 参考 : SICP Exercise 3.77. – 理工系学生の演習日記 (define (int…続きを読む

問題3.76 (define (smooth st) (define (iter st lv) (let ((average (/ (+ (stream-car st) lv) 2.0))) (cons-stream average (iter (stream-cdr st) (stream-car st))))) (iter st (stream-car st))) (define zero-c…続きを読む

問題3.75 Louis の変更した手続きでは平滑化された値（avpt）を次の last-value として渡している。 さらに境界のチェックは平均値同士から算出しないとだめらしい。 したがって修正した make-zero-crossings 手続きは以下のようになる。 (define (make-zero-crossings input-stream last-value last-avpt) …続きを読む