RESIDUE.BAK

/*
void plot(double x,double y,int c,double sx,double sy)
{

	if(fabs(x)<getmaxy())
	{
		if(y>0)
		{
			if(sx+1>fabs(y))
			{
				putpixel(sy+x,sx-y,c);
				lineto(sy+x,sx-(y));

			}
		}
		if(y<0)
		{
			if(sx+1>fabs(y))
			{
				putpixel(sy+x,sx+fabs(y),c);
				lineto(sy+x,sx+fabs(y));
			}
		}
	}
}
*/
/*
int get_n(int n)
{
	cout<<"\n\t Enter the value of n";
	cin>>n;
	return(n);
}
*/
double f=0.05;
	double T=1/f;
	double dx=1;
	//cout<<sx<<"       "<<sy;getch();

	for(double x=-sy;x<sy;x+=dx)
	{
		double y=0;
		/*
		for(int n=0;n<3;n++)
		{
			//y+=coeff(n,f,T,x)*sin(2*M_PI*f*n*x);
			y+=sin(2*M_PI*f*n*x);
		}
		*/
		//term=sin(2*M_PI*f*n*x);
		/*
		int a=100;
		for(int j=0; j<a; j++)
		{
		*/
		//plot(x*j,1000,4,sx,sy);
		//delay(500);
		//plot (x,y*100,4,sx,sy);
		//cou t<<"("<<x<<","<<y<<")"<<endl;
		//ge tch();
		//p lot(x*j,1000,2,sx,sy);
		// delay(500);


	}
	/*
	for(x=-sy;x<sy;x+=dx)
	{
		if(fabs(x*scx)<sy)
		{       y=f(x);
			if(errno!=EDOM)
			{
			y*=scy;
			x*=scx;
			if(x>=l&&x<=u)
			{
				s+=y*dx;
				line(sy+x,sx,sy+x,sx-y);
			}
			plotr(x,y,c,0,len,sx,sy,m);
			asdf++;

			delay(1);


			m=(f(x+dx)-f(x))/dx;
			m=-1/m;
			y/=scy;
			x/=scx;

		}
		}
	} /* */