#N canvas 384 121 780 602 10;
#X obj 257 361 +~;
#X obj 132 112 *~ 1024;
#X obj 257 160 tabread4~ \$0-cos;
#X obj 131 350 tabread4~ \$0-sin;
#X obj 132 421 /~;
#X obj 132 132 +~ 1;
#X obj 480 250 *~ 1024;
#X obj 480 278 +~ 1;
#X obj 410 311 tabread4~ \$0-1+a2;
#X obj 544 310 tabread4~ \$0-norm;
#X obj 133 444 *~;
#X obj 380 258 *~ -2;
#X obj 257 311 *~;
#X obj 134 473 outlet~;
#X obj 379 45 inlet~;
#X obj 489 66 table \$0-1+a2 1027;
#X obj 489 85 table \$0-norm 1027;
#X text 133 503 out = sqrt[(1-a^2)(1+a^2)] * [sin(x) / (1+a^2 - 2a*cos(x))]
;
#X obj 489 105 table \$0-cos 1027;
#X obj 489 127 table \$0-sin 1027;
#N canvas 436 142 729 607 init-tables 0;
#X obj 125 75 until;
#X obj 125 104 f 0;
#X obj 153 105 + 1;
#X obj 125 53 t a b;
#X obj 125 150 / 1024;
#X obj 125 127 t a a;
#X obj 343 417 tabwrite \$0-1+a2;
#X obj 125 12 loadbang;
#X obj 547 122 s \$0-cos;
#X obj 509 175 s \$0-sin;
#X msg 509 148 sinesum 1024 1;
#X msg 547 97 cosinesum 1024 0 1;
#X obj 509 67 loadbang;
#X obj 124 412 tabwrite \$0-norm;
#X msg 125 33 1027;
#X text 106 499 writes f(a) = 1+a^2 rsp. normalization into tables
;
#X obj 124 330 sqrt;
#X obj 125 214 *;
#X obj 179 272 + 1;
#X obj 125 241 t a a;
#X obj 125 267 swap 1;
#X obj 125 290 -;
#X obj 125 310 /;
#X obj 125 193 t a a;
#X text 107 434 sqrt((1-$f1*$f1)/(1+$f1*$f1));
#X text 341 434 1+$f1*$f1;
#X msg 164 75 0;
#X obj 452 361 + 0;
#X connect 0 0 1 0;
#X connect 1 0 2 0;
#X connect 1 0 5 0;
#X connect 2 0 1 1;
#X connect 3 0 0 0;
#X connect 3 1 26 0;
#X connect 4 0 23 0;
#X connect 5 0 4 0;
#X connect 5 1 27 0;
#X connect 7 0 14 0;
#X connect 10 0 9 0;
#X connect 11 0 8 0;
#X connect 12 0 11 0;
#X connect 12 0 10 0;
#X connect 14 0 3 0;
#X connect 16 0 13 0;
#X connect 17 0 19 0;
#X connect 18 0 22 1;
#X connect 18 0 6 0;
#X connect 19 0 20 0;
#X connect 19 1 18 0;
#X connect 20 0 21 0;
#X connect 20 1 21 1;
#X connect 21 0 22 0;
#X connect 22 0 16 0;
#X connect 23 0 17 0;
#X connect 23 1 17 1;
#X connect 26 0 1 1;
#X connect 27 0 6 1;
#X connect 27 0 13 1;
#X restore 489 45 pd init-tables;
#X text 131 550 cf: James A. Moorer: "Synthesis of Complex Audio Spectra"
and Dodge/Jerse: "Computer Music" page 163f;
#X obj 132 89 phasor~ \$1;
#X obj 132 46 inlet;
#N canvas 228 198 627 317 LICENSE-BSD 0;
#X text 121 56 This software is copyrighted by Miller Puckette \, Reality
Jockey Ltd. and others. The terms (the "Standard Improved BSD License")
apply to all files associated with the software unless explicitly disclaimed
in individual files.;
#X text 123 148 See the file LICENSE.txt for the full license text.
;
#X restore 633 30 pd LICENSE-BSD;
#X connect 0 0 4 1;
#X connect 1 0 5 0;
#X connect 2 0 12 0;
#X connect 3 0 4 0;
#X connect 4 0 10 0;
#X connect 5 0 2 0;
#X connect 5 0 3 0;
#X connect 6 0 7 0;
#X connect 7 0 8 0;
#X connect 7 0 9 0;
#X connect 8 0 0 1;
#X connect 9 0 10 1;
#X connect 10 0 13 0;
#X connect 11 0 12 1;
#X connect 12 0 0 0;
#X connect 14 0 6 0;
#X connect 14 0 11 0;
#X connect 22 0 1 0;
#X connect 23 0 22 0;
