|
|
1.1 ! root 1: define(N,50) # size of sample grid in rho-sigma space. ! 2: define(M,11) # number of samples in theta ! 3: define(RMAX,3.) # rho search limit ! 4: define(SMAX,3.) # sigma search limit ! 5: ! 6: procedure ! 7: integer i, j, k ! 8: real theta, sigma, rho, fmin, fmax, outsid, f(N,N) ! 9: real fun ! 10: ! 11: do k=1,M{ ! 12: theta = (k-1)*50./max(1,M-1) ! 13: do j=1,N{ ! 14: sigma = j*SMAX/N ! 15: do i=1,N{ ! 16: rho = i*RMAX/N ! 17: f(i,j) = fun(rho,sigma,theta) ! 18: } ! 19: } ! 20: fmin = 0 ! 21: fmax = 0 ! 22: outsid = -1.e20 ! 23: mvefil(30,theta,N,N,fmin,fmax,outsid,f) ! 24: } ! 25: ! 26: end ! 27: ! 28: real procedure fun ( rho, sigma, theta ) ! 29: real rho, sigma, theta ! 30: real t, f, deg ! 31: deg = 4 * atan(1.) / 180 ! 32: f = sin(rho*10*deg) * cos(sigma*10*deg) ! 33: return( f ) ! 34: end
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.