// generate a .fig file containing a Dirichlet function plot from 0 to 1

#include <cstdio>
#include <cln/integer.h>

const int scale = 18000;

void printline(double x, double y)
{
	printf("2 1 0 1 0 7 50 -1 -1 0.000 0 0 -1 0 0 2\n"
		"\t%d %d %d %d\n", int(scale*x), scale/2, int(scale*x), scale/2-int(scale*y)
	); 
}

int main(int argc, char* argv[])
{
	printf("#FIG 3.2 Generated by make_diri\n"
	"Landscape\n"
	"Center\n"
	"Metric\n"
	"A4\n"
	"100.00\n"
	"Single\n"
	"-3\n"
	"1200 2\n"
	);

	for(unsigned n=2; n<750; ++n)
	{
		for(unsigned z=1; z<n; ++z)
		{
			if( cln::gcd(z,n) > 1 ) continue;

			const double r = 1/double(n);
			printline(z*r, r);
		}
	}
	fprintf(stderr, "Okay.\n");
}