#define XMin -2.20 /*定义坐标及矩阵参数*/
#define XMax 0.60
#define YMin -1.20
#define YMax 1.20
#define MaxIter 30.0
#define Res 160.0

char source[160];
int srcoff,srcseg;

int Iterate(float cx, float cy) /*Julia集向Mandelbrot集的变换*/
{
int Iters;
float x, y;
float x2, y2;
float temp;

x=cx;
x2=SqrFP(x);
y=cy;
y2=SqrFP(y);
Iters=0;
while((Iters<MaxIter) && (x2+y2<4))
{
temp=cx+x2-y2;
y=cy+2*x*y;
y2=SqrFP(y);
x=temp;
x2=SqrFP(x);
++Iters;
}
return(Iters);
}

void Pix(int x, int y, int col) /*在图上绘制一象素点*/
{
int a;

Plot(x, y, col);
Plot(x, Res-y-1, col);
}

void CalcMSet() /*根据迭代次数决定颜色并绘制象素点*/
{
int ix, iy;
int Iters;
float cx, cy;
float dx, dy;
int L1, L2;

Draw(0, 0, Res-1, 0, 35);
Draw(Res-1, 0, Res-1, Res-1, 35);
Draw(Res-1, Res, 0, Res, 35);
Draw(0, Res, 0, 0, 35);
Draw(Res, 0, 2*Res-1, 0, 35);
Draw(2*Res-1, 0, 2*Res-1, Res-1, 35);
Draw(2*Res-1, Res, Res, Res, 35);
L1=(MaxIter-(MaxIter*2)/3);
L2=(MaxIter-(MaxIter*5)/6);
dx=(XMax-XMin)/(Res-1);
dy=(YMax-YMin)/(Res-1);
for(iy=2; iy<=((Res-1)/2); iy++)
{
cy=YMin+iy*dy;
for(ix=2; ix<=Res-3; ix++)
{
cx=XMin+ix*dx;
Iters=Iterate(cx, cy);
if(Iters==MaxIter)
Pix(ix, iy, 143);
else
{
if(Iters>L1)
Pix(ix, iy, 215);
else
{
if(Iters>L2)
Pix(ix, iy, 35);
}
}
}
}
}

void CalcJSet(float cx, float cy) /*利用反转技术计算Julia集*/
{
int xp, yp;
float dx, dy;
float r;
float Theta;
float x, y;

x=0;
y=0;
do
{
dx=x-cx;
dy=y-cy;
if(dx>0)
Theta=atan(dy/dx)*0.5;
else
{