00001 program make_data
00002 c
00003 c Make sample data for mcdraw.
00004 c
00005 parameter (N = 250, ND = 5, PI = 3.14159)
00006 c
00007 open(10,file='DATA.1')
00008 rewind 10
00009 do i = 1,N
00010 xx = 3.*i/float(N)
00011 write(10,*)i,xx,sin(2.*PI*xx),cos(2.*PI*xx),
00012 $ exp(-xx)*sin(2.*PI*xx)
00013 end do
00014 close(10)
00015 c
00016 open(10,file='DATA.2')
00017 rewind 10
00018 do j = 1,N/2
00019 yy = 3.*j/float(N/2)
00020 do i = 1,N/2
00021 xx = 2.*i/float(N/2)
00022 write(10,*)sin(2.*PI*xx)*cos(2.*PI*yy)
00023 end do
00024 end do
00025 close(10)
00026 c
00027 end
