|
|
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.