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