3577l1.orig100644 1061 330 542 6746151742 11113 0ustar mlrljeditset k 0.05 set B 7.5 proc NewValue {Name Value} { puts "$Name [expr $Value]" } scale .k -variable k -label Damping \ -from 1.0 -to 0.0 -resolution 0.01 \ -command "NewValue k" scale .b -variable B -label Forcing \ -from 30.0 -to 0.0 -resolution 0.10 \ -command "NewValue B" pack .k .b -side left 3577l2.orig100644 1061 330 2002 6746152026 11121 0ustar mlrljeditcanvas .duffing -width 210 -height 210 .duffing create oval 50 190 60 200 \ -fill blue -tag movable .duffing create line 5 205 205 205 .duffing create line 5 5 5 205 .duffing create text 20 10 -text "k" .duffing create text 190 190 -text "B" .duffing bind movable { CanvasMark %x %y %W } .duffing bind movable { CanvasDrag %x %y %W } proc CanvasMark { x y w } { global canvas set canvas($w,obj) [$w find closest $x $y] set canvas($w,x) $x set canvas($w,y) $y } proc CanvasDrag { x y w } { global canvas if {$x < 5} { set x [expr 5] } if {$x > 205} { set x [expr 205] } if {$y < 5} { set y [expr 5] } if {$y > 205} { set y [expr 205] } $w move $canvas($w,obj) [expr $x-$canvas($w,x)]\ [expr $y-$canvas($w,y)] set canvas($w,x) $x set canvas($w,y) $y puts stdout "B [expr 0.150 * $x - 0.750]" puts stdout "k [expr -0.005 * $y + 1.025]" flush stdout } pack .duffing 3577l3.orig100664 1007 330 5103 6746152216 11442 0ustar darcyljeditBEGIN { if (MyHost == "") { "uname -n" | getline MyHost; close("uname -n") } if (MyPort == 0) MyPort = 8080 HttpService = "/inet/tcp/" MyPort "/0/0" MyPrefix = "http://" MyHost ":" MyPort SetUpServer() while ("GAWK" < "Perl") { RS = ORS = "\r\n" # header lines are terminated this way Status = 200 # this means OK Reason = "OK" Header = TopHeader Document = TopDoc Footer = TopFooter if (GETARG["Method"] == "GET") { HandleGET() } else if (GETARG["Method"] == "HEAD") { # not yet implemented } else if (GETARG["Method"]!=""){print "bad method", GETARG["Method"] } Prompt = Header Document Footer print "HTTP/1.0", Status, Reason |& HttpService print "Content-length:", length(Prompt) + length(ORS) |& HttpService print ORS Prompt |& HttpService while ((HttpService |& getline) > 0) ; # ignore all the header lines close(HttpService) # stop talking to this client HttpService |& getline # wait for new client request print systime(), strftime(), $0 # do some logging delete GETARG; delete MENUE; delete PARAM GETARG["Method"]=$1; GETARG["URI"]=$2; GETARG["Version"]=$3 for (i=length($2); (substr($2,i,1)!="?") && (i>0); i--) ; if (i > 0) { # is there a "?" indicating a CGI request ? split(substr($2, 1, i-1), MENUE, "[/:]") split(substr($2, i+1), PARAM, "&") for (i in PARAM) { j = index(PARAM[i], "=") GETARG[substr(PARAM[i], 1, j-1)] = substr(PARAM[i], j+1) } } else { # there is no "?", no need for splitting PARAMs split($2, MENUE, "[/:]") } } } ######################################################################## # The site independent part ends here. ######################################################################## function Doc(k, B) { return "
\ \ \ \
Damping k Forcing B

" } function SetUpServer() { k=0.05 ; B = 7.5 TopHeader = "Statistics with GAWK" TopDoc = Doc(k, B) TopFooter = "" } function HandleGET() { if ("k" in GETARG) print "k", k=GETARG["k"] if ("B" in GETARG) print "B", B=GETARG["B"] fflush() Document = Doc(k, B) } 3577l4.orig100640 1061 330 11364 6746123301 11144 0ustar mlrljedit#include #include #include #include #include #include #include #define SAMPLE_RATE 44100 #define BLOCKS_PER_SECOND 50 #if defined(LINUX) #include int SoundFD; #elif defined(IRIX) #include ALconfig Irix_AL_config; ALport SoundFD; #elif defined(CYGWIN) #include #define NONAMELESSUNION #define _stdcall STDCALL typedef HANDLE HTASK; #include #include static HWAVEOUT WaveDevice = NULL; static sem_t WavBufSema; static void CALLBACK WaveCallBack(HWAVE hWave, UINT uMsg, DWORD dwInstance, DWORD dwParam1, DWORD dwParam2) { if (uMsg == WOM_DONE) { WAVEHDR * wh = (WAVEHDR *) dwParam1; HGLOBAL hg; waveOutUnprepareHeader(WaveDevice, wh, sizeof(WAVEHDR)); hg = GlobalHandle(wh->lpData); GlobalUnlock(hg); GlobalFree(hg); hg = GlobalHandle(hg); GlobalUnlock(hg); GlobalFree(hg); sem_post (& WavBufSema); } } #endif int InitSoundSystem(int SampleRate) { int OK = -1; #if defined(LINUX) if ((SoundFD=open("/dev/dsp", O_WRONLY)) >=0) { int OSS_format = AFMT_S16_LE; int OSS_speed = SampleRate; ioctl(SoundFD, SNDCTL_DSP_RESET, 0); if ((ioctl(SoundFD, SNDCTL_DSP_SETFMT, &OSS_format) != -1) && (OSS_format == AFMT_S16_LE) && (ioctl(SoundFD, SNDCTL_DSP_SPEED, &OSS_speed ) != -1) && (OSS_speed == SampleRate)) { OK = SoundFD; } } #elif defined(IRIX) long PV [10] = { AL_LEFT_SPEAKER_GAIN, 10, AL_RIGHT_SPEAKER_GAIN, 10, AL_OUTPUT_RATE, }; PV [5] = AL_RATE_44100; if ((ALsetparams(AL_DEFAULT_DEVICE,PV,6)==0) && ((Irix_AL_config = ALnewconfig()) != NULL)){ ALsetchannels (Irix_AL_config, AL_STEREO); ALsetqueuesize (Irix_AL_config, SAMPLE_RATE); ALsetwidth (Irix_AL_config, AL_SAMPLE_16); if (((SoundFD=ALopenport("Speaker","w", Irix_AL_config)) != NULL)) OK = 0; } #elif defined(CYGWIN) WaveDevice = NULL; sem_init (& WavBufSema, 0, BLOCKS_PER_SECOND); if (waveOutGetNumDevs() > 0) { WAVEFORMATEX outFormatex; outFormatex.wFormatTag = WAVE_FORMAT_PCM; outFormatex.wBitsPerSample = 16; outFormatex.nChannels = 1; outFormatex.nSamplesPerSec = SAMPLE_RATE; outFormatex.nBlockAlign = sizeof(short int); outFormatex.nAvgBytesPerSec = outFormatex.nSamplesPerSec * sizeof(short int); if (waveOutOpen(&WaveDevice, WAVE_MAPPER, &outFormatex, (DWORD) WaveCallBack, 0, CALLBACK_FUNCTION) == MMSYSERR_NOERROR) OK = 0; } #endif return OK; } void SoundOutAndWait(short int * Sample, int HowMany) { int HowManyBytes = HowMany * sizeof(short int); #if defined(LINUX) write(SoundFD, Sample, HowManyBytes); #elif defined(IRIX) ALwritesamps (SoundFD, Sample, HowMany); #elif defined(CYGWIN) HGLOBAL HGheader = GlobalAlloc (GMEM_MOVEABLE | GMEM_ZEROINIT, sizeof(WAVEHDR)); HGLOBAL HGdata = GlobalAlloc (GMEM_MOVEABLE, HowManyBytes); LPWAVEHDR WaveHeaderPtr = GlobalLock(HGheader); WaveHeaderPtr->dwBufferLength = HowManyBytes; WaveHeaderPtr->lpData = GlobalLock(HGdata); CopyMemory(WaveHeaderPtr->lpData, Sample, HowManyBytes); waveOutPrepareHeader(WaveDevice, WaveHeaderPtr, sizeof(WAVEHDR)); waveOutWrite (WaveDevice, WaveHeaderPtr, sizeof(WAVEHDR)); sem_wait (& WavBufSema); #endif } int main (int argc, char *argv[]) { double k=0.05, B=7.5, x1=3.0, dx=4.0, dt=M_PI/100.0, ddt, x0, x2; int i, j = 0; if (InitSoundSystem(SAMPLE_RATE) < 0) exit(0); ddt = dt * dt; x0 = x1 - dx*dt; while (1) { short int SampleBuffer[SAMPLE_RATE / BLOCKS_PER_SECOND]; fd_set stdin_set; struct timeval timeout = { 0, 0 }; FD_ZERO(& stdin_set); /* non-blocking read */ FD_SET(fileno(stdin), & stdin_set); if (select(fileno(stdin)+1, & stdin_set, NULL, NULL, & timeout) > 0) { char line[128], Name[128]; float Value; fgets(line, sizeof(line), stdin); sscanf (line, "%s %f", Name, &Value); if (strcmp(Name,"k")==0) k=(double) Value; if (strcmp(Name,"B")==0) B=(double) Value; } for (i=0; i<(SAMPLE_RATE/BLOCKS_PER_SECOND); i++, j++) { SampleBuffer[i] = 4096 * x1; x2 = (ddt*B * cos(j*dt) + (2+dt*k)*x1 - x0 - ddt*x1*x1*x1) / (1 + dt * k); if ((j%1000) == 0) /* Poincare section */ printf("%f %f\n", x2, (x2-x1)/dt); x0=x1; x1=x2; } SoundOutAndWait(SampleBuffer, SAMPLE_RATE / BLOCKS_PER_SECOND); } return 0; } 3577l5.orig100644 1061 330 654 6746152100 11110 0ustar mlrljeditcanvas .p -width 210 -height 210 .p create line 5 105 205 105 .p create line 105 5 105 205 .p create text 180 90 -text "x" .p create text 120 10 -text "x'" pack .p proc draw { } { if { ([gets stdin line] >= 0) && \ ([scan $line "%f %f" x y] >= 2) } { set x [expr 25 * $x + 105] set y [expr -12.5 * $y + 105] .p create oval $x $y $x $y after 1 { draw } } } after 0 { draw }