#include <stdio.h>
#include <math.h>
#include <dos.h>
#include "tools.h"

float Xmax = 2.8,Xmin = -2.8,Ymax = 2,Ymin = -2, X = 0, Y = 0, Z = 0;
float deltaX,deltaY,Xtemp,Ytemp,Ztemp;
int col,row,j,max_row = 349, max_col = 639,color;
float a = 2.24, b = .43, c = -.65, d = -2.43, e = 1;
long int max_iterations=50000,i;
int OPERATOR = 0;
char ch;
main()
{
	setMode(16);
	deltaX = max_col/(Xmax - Xmin);
	deltaY = max_row/(Ymax - Ymin);
	for (j=0; j<2; j++)
	{
		cls(0);
		for (i=0; i<max_iterations; i++)
		{
			Xtemp = sin(a*Y) - Z*cos(b*X);
			Ytemp = Z*sin(c*X) - cos(d*Y);
			Z = e*sin(X);
			X = Xtemp;
			Y = Ytemp;
			if (j==0)
			{
				col = (X - Xmin)*deltaX;
				row = (Y - Ymin)*deltaY;
			}
			else
			{
				col = (Y - Xmin)*deltaX;
				row = (Z - Ymin)*deltaY;
			}
			if ((col>0) && (col<=max_col) &&
				(row>0) && (row<=max_row))
			{
				color = readPixel(col,row);
				color = (++color)%15 + 1;
				plot(col,row,color);
			}
		}
	getch();
	}
}
               /*ÚÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ¿
                 ³                                                       ³
                 ³      readPixel = Read a Pixel from the Screen         ³
                 ³                                                       ³
                 ³                                                       ³
                 ÀÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÙ*/



int readPixel(int x, int y)
{

	union REGS reg;

	reg.h.ah = 0x0D;
	reg.x.cx = x;
	reg.x.dx = y;
        int86 (0x10,&reg,&reg);
	return (reg.h.al);
}