/*
 * Gear Curve
 * x	=	r*cos*theta	
 * y	=	r*sin*theta
 * r=a+(1/b)*tanh(b*sin(n*theta)), 	
 *where tanhx is the hyperbolic tangent. Plots above show gear curves for a=1, b=10, and n=1 to 12. 
*/

//g++ -Wall -o "%e" "%f" -lboost_iostreams -lboost_system -lboost_filesystem -lSDL -lSDL_image -lSDL_ttf -lSDL_mixer -lGLU -lglut -lXi -lXmu `pkg-config --libs opencv` 
#include <iostream>
#include <fstream>
#include <math.h>
#include "/usr/include/gnuplot-iostream.h"

using namespace std;

void init(){
	
	Gnuplot gp;
	
	}

void gear(double teeth){
	Gnuplot gp;
	ofstream plot;
    plot.open("gear.dat");
 // x	=	r*cos*theta	
 // y	=	r*sin*theta
 // r=a+(1/b)*tanh(b*sin(n*t))
	
	double a =1.5;
	double b =12;
	const double pi = 3.14159265359;
	double r =0;
	double cR=0;
	double n = 1;// teeth?
	n=teeth;
	
	double xx=0;
	double yy = 0;
	double cx = 0;
	double cy = 0;
	
		for(int t=0; t<=360; t++){
			double theta =(t*pi)/180;
	    	 
	        r=a+(1/b)*tanh(b*sin(n*theta));
	        xx = r * cos(theta);
	        yy = r * sin(theta);
	        cR=a+(1-b)*tanh(b*sin(0*theta));
	        cx=cR*cos(theta);
	        cy=cR*sin(theta);
	  plot << r <<" " << xx << " " << yy << " " << cR << " " <<cx  << " "<< cy << endl; 
	  //plot << r << " " << sin(theta) << endl;   
       
     }
     plot.close();
     gp << "unset polar\n";
     gp << "plot \"gear.dat\" using 2:3 with lines\n"; 
    // gp << "plot \"gear.dat\" using 5:6 with lines\n"; 
     plot.close();
	}

int main(int argc, char **argv){
	//init();
	cout << "Enter Number of teeth: ";
	double k = 5;
	cin >> k;
	
	gear(k);
	
	
	
	return 0;
	}