/*Translated by FOR_C, v3.4.2 (-), on 07/09/115 at 08:33:08 */
/*FOR_C Options SET: ftn=u io=c no=p op=aimnv pf=,p_cwofz s=dbov str=l x=f - prototypes */
#include <math.h>
#include "fcrt.h"
#include <stdio.h>
#include <stdlib.h>
#include "p_cwofz.h"
/*     program DRCWOFZ
 *>> 1996-06-10 DRCWOFZ Krogh Moved formats up, M77CON changes for C.
 *>> 1992-03-13 DRCWOFZ WV Snyder Created from DRCWOFZ
 *>> 1991-11-20 DRCWOFZ WV Snyder Original Code
 * Conversion should only be done from "Z" to "C" for processing to C.
 *--C replaces "?": DR?WOFZ, ?WOFZ
 *
 *     Demonstration driver for CWOFZ.
 * */
 
int main( )
{
	long int _n, iflag, ix, iy;
	float w[2], z[2];
		/* OFFSET Vectors w/subscript range: 1 to dimension */
	float *const W = &w[0] - 1;
	float *const Z = &z[0] - 1;
		/* end of OFFSET VECTORS */
 
 
	/*     Evaluate w(z) at 25 points [0..4] x [0..4] in the complex plane.
	 * */
	printf("  x  y             Re w                      Im w\n");
	for (ix = 0; ix <= 4; ix++)
	{
		for (iy = 0; iy <= 4; iy++)
		{
			Z[1] = (float)( ix );
			Z[2] = (float)( iy );
			cwofz( z, w, &iflag );
			printf("%3ld%3ld", ix, iy);
			for(_n=0L; _n < sizeof(w)/sizeof(float); _n++)
				printf("%27.18g", w[_n]);
			printf("\n");
		}
	}
	exit(0);
} /* end of function */