#include <stdio.h>
#include <math.h>

int main()
{
	float f, wd, dw, d1;
	float a, aa, ab, ac;
	float b, ba, bb, bc;
	float c, ca, cb, cc;
	float d, da, db, dc;
	float e;
	float wli, wf, wfi;
	int u = 0;

	printf( "Program to calculate the dimensions of a Moxon Rectangle.\n" );
	printf( "All equations correlated to NEC antenna modeling software for wire diameters" );
	printf( "     from 1E-5 to 1E-2 wavelengths.\n" );
	printf( "L. B. Cebik, W4RNL\n" );
	printf( "Enter Desired Frequency in MHz: " );
	scanf( "%f",&f ); 
	while( u > 3 || u < 1 ) {
		printf( "Select Units for Wire Diameter in\n\t1. Inches\n\t2. Millimeters\n\t3. Wavelengths\n" );
		printf( "Choose 1. or 2. or 3. " );
		scanf( "%d", &u );
	}
	printf( "Enter Wire Diameter in your Selected Units: " );
	scanf( "%g", &wd );
	if( u == 1 ) {
		wli = 11802.71 / f;
		dw = wd / wli;
	}
	if( u == 2 ) {
		wli = 299792.5 / f;
		dw = wd / wli;
	}
	if( u == 3 )
		dw = wd;
	printf( "Wire Diameter in Wavelengths: %g\n", dw );
	d1 = 0.4342945 * log( dw );
	if( d1 < -6 )
		printf( "Wire diameter less than 1e-6 wavelengths: results uncertain!\n" );
	if( d1 > -2 )
		printf( "Wire diameter greater than 1e-2 wavelengths: results uncertain!\n" );
	aa = -0.0008571428571;
	ab = -0.009571428571;
//	printf( "aa = %g\nab = %g\n", aa, ab );
	ac = 0.3398571429;
	a = ( aa * pow( d1, 2 )) + ( ab * d1 ) + ac;
	ba = -0.002142857143;
	bb = -0.02035714286;
	bc = 0.008285714286;
	b = ( ba * pow( d1, 2 )) + ( bb * d1 ) + bc;
	ca = 0.001809523381;
	cb = 0.01780952381;
	cc = 0.05164285714;
	c = ( ca * pow( d1, 2 )) + ( cb * d1 ) + cc;
	da = 0.001;
	db = 0.07178571429;
	d = ( da * d1 ) + db;
	e = b + c + d;

	printf( "Moxon Dimensions in Wavelengths:\n" );
	printf( "A = %.4g\n", a );
	printf( "B = %.4f\n", b );
	printf( "C = %.4f\n", c );
	printf( "D = %.4f\n", d );
	printf( "E = %.4f\n", e );

	wf = 983.5592 / f;
	wfi = wf * 12;
	printf( "Wavelength = %.4f Feet or %.4f Inches\n", wf, wfi );
	printf( "Dimensions in Feet and Inches\n" );

	printf( "A = %.4f Feet or %.4f Inches\n", a * wf, a * wfi );
	printf( "B = %.4f Feet or %.4f Inches\n", b * wf, b * wfi );
	printf( "C = %.4f Feet or %.4f Inches\n", c * wf, c * wfi );
	printf( "D = %.4f Feet or %.4f Inches\n", d * wf, d * wfi );
	printf( "E = %.4f Feet or %.4f Inches\n", e * wf, e * wfi );
}