#!perl -w
# CABebez.pl
# (c) Copyright 2006 by H. M�ller (mollerh@math.uni-muenster.de).
# Version 1.2 for Cabri-g�om�tre II with MacOS 9.x, Virtual Printer as PostScript driver, MacPerl 5.6, and LATEX-packages 'pict2e' and 'ebezier'.
# This program may be distributed and/or modified under the conditions of the LaTeX Project Public License, either version 1.3 of this license or (at your option) any later version.
# The latest version of this license is in http://www.latex-project.org/lppl.txt.
# This program has the LPPL maintenance status "maintained". The Current Maintainer of this program is H. M�ller.
#
use POSIX('ceil','floor');
#________________________________________________________
# Definable by user:
# Unitlength in pt:
$ul = 1.0;
# Fill factor (for filling with magnification up to 500%):
$fillf = 5;
# B�zier factor:
$bezf = 2.0;
# Point factor:
$pointf = 0.3;
#________________________________________________________
# Constants:
# Color names:
$yellow = "0.9843900.9511410.020249";
$orange = "1.0000000.3927370.009949";
$red = "0.8649270.0342110.025910";
$purple = "0.9486080.0325630.519234";
$violet = "0.2769050.0000000.645487";
$navy = "0.0000000.0000000.828138";
$blue = "0.0088040.6692610.917967";
$green = "0.1215990.7170980.078874";
$darkgreen = "0.0000000.3933010.069093";
$darkbrown = "0.3359430.1742730.020081";
$brown = "0.5657890.4428780.227359";
#________________________________________________________
# Further Constants:
# Pi:
$Pi = "3.14159265358979";
# Constant in cubic B�zier curves for quarters of a circle:
$l90 = "0.552284749830794";
# Constants in dotted figures:
$uli = sp(4 / $ul);
$ule = sp(0.8 / $ul);
#________________________________________________________
@lines = <>;
do {
	$_ = $lines[$i++];
	if (/ setrgbcolor \s/o) {
		s/ //go;
		s/setrgbcolor\s/ /o;
		$c = $_;
		$_ = $lines[$i++];
		s/ moveto//o;
		s/lineto stroke/stroke/o;
		s/curveto stroke/curveto/o;
		s/ setlinewidth stroke//o;
		s/ lineto//go;
		if (/stroke/o) {
			$line[++$#line] = $c.$_;
		}
		elsif (/closepath fill/o) {
			$vector[++$#vector] = $_;
		}
		elsif (/arc /o) {
			$circle[++$#circle] = $c.$_;
		}
		elsif (/arcn/o) {
			$arc[++$#arc] = $c.$_;
		}
		elsif (/curveto/o) {
			do {
				$conic[++$#conic] = $c.$_;
				$_ = $lines[$i++];
				s/ moveto//o;
				s/curveto stroke/curveto/o;
			}
			until $_ !~ /curveto/o;
		}
	}
}
until $i == $#lines;
#
$pflag = 1;
$sflag = 1;
$thicknessflag = 1;
$coun = 0;
$bmax = 500;
$xtex = "";
$mtex = "";
$btex = "\\documentclass{article}\n\\usepackage{ebezier}\n\n";
#________________________________________________________
# Lines and polygons
$cflag = 1;
foreach (@line) {
	@coo = split;
	$co0 = $coo[0];
	$co2 = (-1) * $coo[2];
	$co4 = (-1) * $coo[4];
	if ($co0 ne $yellow) {
		if ($cflag) {
			$xtex .= "%Lines, arrows, polygons and B�zier curves\n";
			$cflag = 0;
		}
		if ($co0 ne $blue) {
			bound($coo[1],$co2);
			bound($coo[3],$co4);
		}
		if (($co0 ne $red) and ($co0 ne $blue) and ($co0 ne $brown)
			and ($co0 ne $darkbrown) and ($co0 ne $navy)) {
			lin($co0,$coo[1] * $ul,$co2 * $ul,$coo[3] * $ul,$co4 * $ul);
		}
		if (($co0 ne $green) and ($co0 ne $darkgreen) and ($co0 ne $violet)) {
			if ($pflag) {
				$cb1 = $coo[1];
				$cb2 = $co2;
				$pol = $co0." ".$cb1." ".$cb2;
				$pflag = 0;
			}
			else {
				$pol .= " ".$coo[1]." ".$co2;
				if (abs($coo[3] - $cb1) + abs($co4 - $cb2) < 2.0E-6) {
					$poly[++$#poly] = $pol;
					$pflag = 1;
				}
			}
		}
	}
}
#
foreach (@poly) {
	@po = split;
	$p0 = $po[0];
	$pon = $#po;
	if (($p0 eq $red) or ($p0 eq $purple) or ($p0 eq $darkbrown)
		or ($p0 eq $orange) or ($p0 eq $brown)) {
		if ($pon == 6) {
			tri($p0,$po[1],$po[2],$po[3],$po[4],$po[5],$po[6]);
		}		
		elsif ($pon == 8) {
			($p0,$u1,$v1,$u2,$v2,$u3,$v3,$u4,$v4) = @po;
			$s1 =	abs($u1 - $u4) + abs($u2 - $u3) + abs($v1 - $v2) + abs($v3 -$v4);
			$s2 =	abs($u1 - $u2) + abs($u3 - $u4) + abs($v1 - $v4) + abs($v2 -$v3);
			if (($s1 < 4.0E-6) or ($s2 < 4.0E-6)) {
				bound($u1,$v1);
				bound($u3,$v3);
				rect($p0,$u1,$v1,$u2,$v2,$u3,$v3,$u4,$v4);
			}
			else {
				bound($u1,$v1);
				bound($u2,$v2);
				bound($u3,$v3);
				bound($u4,$v4);
				tri($p0,$u1,$v1,$u2,$v2,$u3,$v3);
				tri($p0,$u1,$v1,$u3,$v3,$u4,$v4);
			}
		}
		elsif ($pon > 8) {
			bound($po[1],$po[2]);
			for (my $j = 3; $j <= $pon - 3; $j += 2) {
				bound($po[$j],$po[$j + 1]);
				tri($p0,$po[1],$po[2],$po[$j],$po[$j + 1],$po[$j + 2],$po[$j + 3]);
			}
			bound($po[$pon - 1],$po[$pon]);
		}
	}
	elsif ($p0 eq $navy) {
		for (my $j = 1; $j <= $pon - 3; $j +=2) {
			lin($p0,$po[$j] * $ul,$po[$j + 1] * $ul,$po[$j + 2] * $ul,$po[$j + 3] * $ul);
		}
	}
	elsif ($p0 eq $blue) {
		if ($pon == 4) {
#Text marker and B�zier curves
		$coun++;
		$bo3 = $po[1] + ($po[3] - $po[1]) / $ul;
		$bo4 = $po[2] + ($po[4] - $po[2]) / $ul;
		bound($po[1],$po[2]);
		bound($bo3,$bo4);
		$po1 = sp($po[1]);
		$po2 = sp($po[2]);
		$mtex .= "\\put(".$po1.",".$po2."){".$coun."}\n";
			}	
		elsif ($pon == 6) {
			$xtex .= "%Quadratic B�zier curve\n";
			qbez($po[1],$po[2],$po[3],$po[4],$po[5],$po[6]);
		}
		elsif ($pon == 8) {
			$xtex .= "%Cubic B�zier curve\n";
			cbez($po[1],$po[2],$po[3],$po[4],$po[5],$po[6],$po[7],$po[8]);
		}
	}
}
#________________________________________________________
# Arrows
foreach (@vector) {
	@ve = split;
	if ($cflag) {
		$xtex .= "%Arrows\n";
		$cflag = 0;
	}	
	$ve[1] = (-1) * $ve[1];
	$ve[3] = (-1) * $ve[3];
	$ve[5] = (-1) * $ve[5];
	$ve[7] = (-1) * $ve[7];
	bound($ve[0],$ve[1]);
	bound($ve[2],$ve[3]);
	bound($ve[6],$ve[7]);
	$vu0 = $ve[4] + ($ve[0] - $ve[4]) / $ul;
	$vu1 = $ve[5] + ($ve[1] - $ve[5]) / $ul;
	$vu2 = $ve[4] + ($ve[2] - $ve[4]) / $ul;
	$vu3 = $ve[5] + ($ve[3] - $ve[5]) / $ul;
	$vu6 = $ve[4] + ($ve[6] - $ve[4]) / $ul;
	$vu7 = $ve[5] + ($ve[7] - $ve[5]) / $ul;
	tri($red,$vu0,$vu1,$vu2,$vu3,$ve[4],$ve[5]);
	tri($red,$vu0,$vu1,$vu6,$vu7,$ve[4],$ve[5]);
}
#________________________________________________________
# Conics
$cflag = 1;
foreach (@conic) {
	@po = split;
	$p0 = $po[0];
	if ($p0 ne $yellow) {
		if ($cflag) {
			$xtex .= "%Conics\n";
			$cflag = 0;
		}
		$po[2] = (-1) * $po[2];
		$po[4] = (-1) * $po[4];
		$po[6] = (-1) * $po[6];
		$po[8] = (-1) * $po[8];
		cbez($po[1],$po[2],$po[3],$po[4],$po[5],$po[6],$po[7],$po[8]);
	}
}
#________________________________________________________
# Circles, halves and quarters of circles
$cflag = 1;
$aflag = 1;
foreach (@circle) {
	@po = split;
	$p0 = $po[0];
	if ($p0 ne $yellow) {
		$po[2] = (-1) * $po[2];
		if ($po[4] > 1.0E-3 or abs($po[5] - 360) > 1.0E-3) {
			if ($aflag) {
				$xtex .= "%Arcs\n";
				$aflag = 0;
			}
			$arce = ($po[4] > 0) ? 360 - $po[4] : 0;
			$arcb = ($po[5] > 0) ? 360 - $po[5] : 0;
			$darc = $arce - $arcb;
			if ($darc < 0) {$darc += 360}
			$quar = int($darc / 90);
			$lq = ($p0 eq $darkgreen) ? $pointf : $bezf;
			if ($quar > 0) {
				for (my $k = 1; $k <= $quar; $k++) {
					quart($po[1],$po[2],$po[3],$arcb,$lq);
					$arcb += 90;
					if ($arcb > 360) {$arcb -= 360}
				}
			}
			if ($darc > $quar * 90) {
				arc($po[1],$po[2],$po[3],$arcb,$arce,$lq);
			}
		}
		else {
			if ($cflag) {
			$xtex .= "%Circles, halves and quarters of circles\n";
			$cflag = 0
			}
			if ($p0 eq $navy) {
				quart($po[1],$po[2],$po[3],0,$bezf);
				quart($po[1],$po[2],$po[3],90,$bezf);
				quart($po[1],$po[2],$po[3],180,$bezf);
				quart($po[1],$po[2],$po[3],270,$bezf);
			}
			elsif ($p0 eq $purple) {
				quart($po[1],$po[2],$po[3],90,$bezf);
				quart($po[1],$po[2],$po[3],180,$bezf);
			}
			elsif ($p0 eq $red) {
				quart($po[1],$po[2],$po[3],0,$bezf);
				quart($po[1],$po[2],$po[3],270,$bezf);
			}
			elsif ($p0 eq $orange) {
				quart($po[1],$po[2],$po[3],180,$bezf);
				quart($po[1],$po[2],$po[3],270,$bezf);
			}
			elsif ($p0 eq $darkbrown) {
				quart($po[1],$po[2],$po[3],0,$bezf);
				quart($po[1],$po[2],$po[3],90,$bezf);
			}
			elsif ($p0 eq $blue) {
				quart($po[1],$po[2],$po[3],180,$bezf);
			}
			elsif ($p0 eq $green) {
				quart($po[1],$po[2],$po[3],90,$bezf);
			}
			elsif ($p0 eq $brown) {
				quart($po[1],$po[2],$po[3],270,$bezf);
			}
			elsif ($p0 eq $violet) {
				quart($po[1],$po[2],$po[3],0,$bezf);
			}
			elsif ($p0 eq $darkgreen) {
				quart($po[1],$po[2],$po[3],0,$pointf);
				quart($po[1],$po[2],$po[3],90,$pointf);
				quart($po[1],$po[2],$po[3],180,$pointf);
				quart($po[1],$po[2],$po[3],270,$pointf);
			}
		}
	}
}
#________________________________________________________
# Arcs
$aflag = 1;
foreach (@arc) {
	@po = split;
	if ($aflag) {
		$xtex .= "%Arcs\n";
		$aflag = 0;
	}
	$po[2] = (-1) * $po[2];
	$arcb = ($po[4] > 0) ? 360 - $po[4] : 0;
	$arce = ($po[5] > 0) ? 360 - $po[5] : 0;
	$darc = $arce - $arcb;
	if ($darc < 0) {$darc += 360}
	$quar = int($darc / 90);
	$lq = ($po[0] eq $darkgreen) ? $pointf : $bezf;
	if ($quar > 0) {
		for (my $k = 1; $k <= $quar; $k++) {
			quart($po[1],$po[2],$po[3],$arcb,$lq);
			$arcb += 90;
			if ($arcb > 360) {$arcb -= 360}
		}
	}
	if ($darc > $quar * 90) {
		arc($po[1],$po[2],$po[3],$arcb,$arce,$lq);
	}
}
#________________________________________________________
# Frame
if ($bmax > 500) {
	$bmax = 100 * ceil($bmax / 100);
	$btex .= "\\renewcommand{\\qbeziermax}{".$bmax."}\n";
}
if ($xtex . $mtex  ne "") {
$xtex = $btex."\\begin{document}\n\n\\setlength{\\unitlength}{".
$ul."pt}\n"."\\begin{picture}(".ceil(($xmax - $xmin)).",".
ceil(($ymax - $ymin)).")(".floor($xmin).",".floor($ymin).")\n".
"\\linethickness{0.8pt}\n"."\\thicklines\n".$xtex;
$xtex .= $mtex."\\end{picture}\n\n\\end{document}";
}
print $xtex."\n";
#
#________________________________________________________
# Lines
sub lin {
	my ($c,$xb,$yb,$xe,$ye) = @_;
	if (($c eq $green) or ($c eq $violet) or ($c eq $orange) or
			($c eq $purple) or ($c eq $navy)) {
		$dx = $xe - $xb;
		$dy = $ye - $yb;
		$xbu = sp($xb / $ul);
		$ybu = sp($yb / $ul);
		if (abs($dx) < 1.0E-6) {
			$sy = $dy <=> 0;
			$leu = sp(abs($dy / $ul));
			$xtex .= "\\put(".$xbu.",".$ybu."){\\line(0,".$sy."){".$leu."}}\n";
		}
		elsif (abs($dy) < 1.0E-6) {
			$sx = $dx <=> 0;
			$leu = sp(abs($dx / $ul));
			$xtex .= "\\put(".$xbu.",".$ybu."){\\line(".$sx.",0){".$leu."}}\n";
		}
		else {
			$flag1 = 1;
			foreach $j (0,1,2,3,4,5,6,8,10,12,13,15,16,18,20,22,24,25,26,27,29,30,34) {
				$num = int($j / 6) + 1;
				$den = ($j % 6) + 1;
				if ($den < $num) {
					$lem = (int(10 * $den / $num) + 1) / $ul;
				}
				else {
					$lem = 10 / $ul;
				}
				if (abs((abs($dy / $dx)) - $num / $den) < 5.0E-4) {
					if ($dy < 0) {
						$num = -$num;
					}
					if ($dx < 0) {
						$den = -$den;
					}
					$le = abs($dx);
					$leu = sp(abs($dx) / $ul - 0.1);	
					if ($leu > $lem) {
						$xtex .= "\\put(".$xbu.",".$ybu."){\\line(".$den.",".$num."){".$leu."}}\n";
						$flag1 = 0;
						last;
					}
					else {
						last;
					}
				}
			}
			if ($flag1) {
				if (not $thicknessflag) {
					$xtex .= "\\linethickness{0.8pt}\n";
					$thicknessflag = 1;
				}
				$len = floor($bezf * (sqrt(($xe - $xb)**2 + ($ye - $yb)**2)));
				if ($len > 0) {
					if ($len > $bmax) {$bmax = $len}
					$xeu = sp($xe / $ul);
					$yeu = sp($ye / $ul);
					$xtex .= "\\Lbezier[".$len."](".$xbu.",".$ybu.")(".$xeu.",".$yeu.")\n";
				}
			}
		}
	}
	elsif ($c eq $darkgreen) {
	# Dotted line
		if ($thicknessflag) {
			$xtex .= "%Dotted line\n\\linethickness{1.2pt}\n";
			$thicknessflag = 0;
		}
		$len = floor($pointf * (sqrt(($xe - $xb)**2 + ($ye - $yb)**2)));
		if ($len > 0) {
			if ($len > $bmax) {$bmax = $len}
				$xbu = sp($xb / $ul);
				$ybu = sp($yb / $ul);
				$xeu = sp($xe / $ul);
				$yeu = sp($ye / $ul);
			$xtex .= "\\Lbezier[".$len."](".$xbu.",".$ybu.")(".$xeu.",".$yeu.")\n";
		}
	}
}
#________________________________________________________
# Triangles
sub tri {
my ($q0,$qx1,$qy1,$qx2,$qy2,$qx3,$qy3) = @_;
	$qx1 = $qx1 * $ul;
	$qy1 = $qy1 * $ul;
	$qx2 = $qx2 * $ul;
	$qy2 = $qy2 * $ul;
	$qx3 = $qx3 * $ul;
	$qy3 = $qy3 * $ul;
	if ($q0 eq $red) {
# Filled triangle
		%ha = ($qx1,$qy1,$qx2+1e-07,$qy2+1e-07,$qx3+2e-07,$qy3+2e-07);
		@hb = ();
		@hc = ();
		foreach (sort { $ha{$a} <=> $ha{$b} } keys %ha) {
			$hb[++$#hb] = $_;
			$hc[++$#hc] = $ha{$_};
		}
		($qx1,$qx2,$qx3) = @hb;
		($qy1,$qy2,$qy3) = @hc;
		$lin = int(($qy3 - $qy1) * $fillf);
		$xtex .= "%Filled triangle\n\\linethickness{0.1pt}\n";
		$si = ($qy3 - $qy1) * $qx2 - ($qx3 - $qx1) * $qy2 - $qx1 * $qy3 + $qx3 * $qy1 <=> 0;
		$dex = ($qx3 - $qx1) / ($qy3 - $qy1) / $fillf;
		$d1 = $qy2 - $qy1;
		$d2 = $qy3 - $qy2;
		if ($d1 >= 1.0E-3) {
			$fx1 = ($qx2 - $qx1) / $d1;
			$sx1 = $qx1 - $qy1 * $fx1;
		}
		if ($d2 >= 1.0E-3) {
			$fx2 = ($qx3 - $qx2) / $d2;
			$sx2 = $qx2 - $qy2 * $fx2;
		}
		for ($k = 1; $k <= $lin; $k++) {
			$xb = $qx1 + $k * $dex;
			$yb = $qy1 + $k / $fillf;
			if ($yb <= $qy2) {
				if ($d1 >= 1.0E-3) {
					$leu = sp((abs($sx1 + $yb * $fx1 - $xb) + 0.5) / $ul);
				}
				else {
					$leu = sp((abs($qx2 - $qx1) + 0.5)/ $ul);
				}
			}
			else {
				if ($d2 >= 1.0E-3) {
					$leu = sp((abs($sx2 + $yb * $fx2 - $xb) + 0.5)/ $ul);
				}
				else {
					$leu = sp((abs($qx3 - $qx2) + 0.5)/ $ul);
				}
			}
			if ($si > 0) {
				$xbu = sp($xb / $ul);
			}
			else {
				$xbu = sp($xb / $ul);
			}
			$ybu = sp($yb / $ul);
			$xtex .= "\\put(".$xbu.",".$ybu."){\\line(".$si.",0){".$leu."}}\n";
		}
		$xtex .= "\\linethickness{0.8pt}\n";
	}
	elsif (($q0 eq $purple) or ($q0 eq $darkbrown)) {
# Dotted triangle
		%ha = ($qx1,$qy1,$qx2+1e-07,$qy2+1e-07,$qx3+2e-07,$qy3+2e-07);
		@hb = ();
		@hc = ();
		foreach (sort { $ha{$a} <=> $ha{$b} } keys %ha) {
			$hb[++$#hb] = $_;
			$hc[++$#hc] = $ha{$_};
		}
		($qx1,$qx2,$qx3) = @hb;
		($qy1,$qy2,$qy3) = @hc;
		$xtex .= "%Dotted triangle\n";
		$si = ($qy3 - $qy1) * $qx2 - ($qx3 - $qx1) * $qy2 - $qx1 * $qy3 + $qx3 * $qy1 <=> 0;
		$dy1 = 2 * ceil($qy1 / 2);
		$dy3 = 2 * floor($qy3 / 2);
		$lin = $dy3 - $dy1;
		$dex = ($qx3 - $qx1) / ($qy3 - $qy1);
		$xbh = $qx1 + ($dy1 - $qy1 - 2.0) * $dex;
		$dex = 2 * $dex;
		$d1 = $qy2 - $qy1;
		$d2 = $qy3 - $qy2;
		if ($d1 >= 1.0E-3) {
			$fx1 = ($qx2 - $qx1) / $d1;
			$sx1 = $qx1 + ($dy1 - $qy1) * $fx1;
		}
		if ($d2 >= 1.0E-3) {
			$fx2 = ($qx3 - $qx2) / $d2;
			$sx2 = $qx2 + ($dy1 - $qy2) * $fx2;
		}
		for ($k = 0; $k <= $lin; $k += 2) {
			$qy = $dy1 + $k;
			$xbh = $xbh + $dex;
			($si > 0) ? ($xb = $xbh) : ($xe = $xbh);
			if ($qy <= $qy2) {
				($d1 >= 1.0E-3) ? ($xeh = $sx1 + $k * $fx1) : ($xeh = $qx1);
			}
			else {
				($d2 >= 1.0E-3) ? ($xeh = $sx2 + $k * $fx2) : ($xeh = $qx2);
			}
			($si > 0) ? ($xe = $xeh) : ($xb = $xeh);
			$xb = 2 * ceil($xb / 2);
			$xbd = $xb + (($xb + $qy) % 4);
			($xe >= $xbd) ? ($num = floor(($xe - $xbd) / 4) + 1) : ($num = 0);
			if (not $thicknessflag) {
				$xtex .= "\\linethickness{0.8pt}\n";
				$thicknessflag = 1;
			}
			$xbu = sp($xbd / $ul);
			$qyu = sp($qy / $ul);
			$xtex .= "\\multiput(".$xbu.",".$qyu.")(".$uli.",0){".$num;
			$xtex .= "}{\\line(1,0){".$ule."}}\n";
		}
	}
	elsif (($q0 eq $orange) or ($q0 eq $brown)) {
#Hatched triangle
		$xtex .= "%Hatched triangle\n";
		$d1 = $qx1 - $qy1;
		$d2 = $qx2 - $qy2;
		$d3 = $qx3 - $qy3;
		$qd1 = $d1;
		$qd2 = $d2;
		$qd3 = $d3;
		%ha = ($qx1,$qd1,$qx2+1e-07,$qd2+1e-07,$qx3+2e-07,$qd3+2e-07);
		@hb = ();
		foreach (sort { $ha{$a} <=> $ha{$b} } keys %ha) {
			$hb[++$#hb] = $_;
		}
		($qx1,$qx2,$qx3) = @hb;
		%ha = ($qy1,$qd1,$qy2+1e-07,$qd2+1e-07,$qy3+2e-07,$qd3+2e-07);
		@hb = ();
		@hc = ();
		foreach (sort { $ha{$a} <=> $ha{$b} } keys %ha) {
			$hb[++$#hb] = $_;
			$hc[++$#hc] = $ha{$_};
		}
		($qy1,$qy2,$qy3) = @hb;	
		($d1,$d2,$d3) = @hc;
		$si = (-$qy3 + $qy1) * $qx2 + ($qx3 - $qx1) * $qy2 + $qx1 * $qy3 - $qx3 * $qy1 <=> 0;
		$p1 = 4 * ceil($d1 / 4);
		$p2 = 4 * floor($d2 / 4);
		$p3 = 4 * floor($d3 / 4);
		$fx1 = ($qx1 - $qx3) / ($d3 - $d1);
		$sx1 = $qx3 + $fx1 * $d3;
		$fy1 = ($qy1 - $qy3) / ($d3 - $d1);
		$sy1 = $qy3 + $fy1 * $d3;
		$d21 = $d2 - $d1;
		$d32 = $d3 - $d2;
		if ($d21 >= 1.0E-3) {
			$fx2 = ($qx1 - $qx2) / $d21;
			$sx2 = $qx2 + $fx2 * $d2;
		}
		if ($d32 >= 1.0E-3) {
			$fx3 = ($qx2 - $qx3) / $d32;
			$sx3 = $qx3 + $fx3 * $d3;
		}
		for ($k = $p1; $k <= $p3; $k += 4) {
			$xbk = $sx1 - $k * $fx1;
			$ybk = $sy1 - $k * $fy1;
			if ($k <= $p2) {
				($d21 < 1.0E-3) ? ($le = abs($qx2 - $qx1)) : ($le = abs($sx2 - $k * $fx2 - $xbk));
			}
			else {
				($d32 < 1.0E-3) ? ($le = abs($qx3 - $qx2)) : ($le = abs($sx3 - $k * $fx3 - $xbk));
			}
			lin($green,$xbk,$ybk,$xbk + $si * $le,$ybk + $si * $le);
		}
	}
}
#________________________________________________________
# Rectangles
sub rect {
	my ($q0,$x1,$y1,$x2,$y2,$x3,$y3,$x4,$y4) = @_;
	$dx = abs($x2 - $x1);
	$dx = ($dx < 1.0E-6) ? abs($x3 -$x2) : $dx;
	$dx = $dx * $ul;
	$xb = ($x1 < $x2) ? (($x1 < $x3) ? $x1 : $x3) : (($x2 < $x3) ? $x2 : $x3);
	$xb = $xb * $ul;
	$dy = abs($y2 - $y1);
	$dy = ($dy < 1.0E-6) ? abs($y3 -$y2) : $dy;
	$dy = $dy * $ul;
	$yb = ($y1 < $y2) ? (($y1 < $y3) ? $y1 : $y3) : (($y2 < $y3) ? $y2 : $y3);
	$yb = $yb * $ul;
	$xe = $xb + $dx;
	$ye = $yb + $dy;
# Filled rectangle
	if ($q0 eq $red) {
		$xtex .= "%Filled rectangle\n\\linethickness{0.1pt}\n";
		$lin = 5 * $dx;
		if ($dy <= $dx) {
			$ybf = $yb - 0.2;
			$xbu = sp($xb / $ul);
			$dxu = sp($dx / $ul);
			for ($k = 0; $k <= $lin; $k++) {
				$ybf += 0.2;
				$ybu = sp($ybf / $ul);
				$xtex .= "\\put(".$xbu.",".$ybu."){\\line(1,0){".$dxu."}}\n";
			}
		}
		else {
			$xbf = $xb - 0.2;
			$ybu = sp($yb / $ul);
			$dyu = sp($dy / $ul);
			for ($k = 0; $k <= $lin; $k++) {
				$xbf += 0.2;
				$xbu = sp($xbf / $ul);
			$xtex .= "\\put(".$xbu.",".$ybu."){\\line(0,1){".$dyu."}}\n";
			}
		}
		$xtex .= "\\linethickness{0.8pt}\n";
	}
# Dotted rectangle
	elsif (($q0 eq $purple) or ($q0 eq $darkbrown)) {
		$xtex .= "%Dotted rectangle\n";
		$xbb = 2 * ceil($xb / 2);
		$ybb = 2 * ceil($yb / 2);
		for ($k = 0; $k <= 2; $k += 2) {
			$ybd = $ybb + $k;
			$xbd = $xbb + (($xbb + $ybd) % 4);
			$numx = floor(($xe - 0.13 - $xbd) / 4) + 1;
			$numy = floor(($ye - 0.13 - $ybd) / 4) + 1;
			$xbu = sp($xbd / $ul);
			$ybu = sp($ybd / $ul);
			$xtex .= "\\multiput(".$xbu.",".$ybu.")(".$uli.",0){".$numx;
			$xtex .= "}{\\begin{picture}(0,0)\\multiput(0,0)(0,".$uli."){".$numy;
			$xtex .= "}{\\line(1,0){".$ule."}}\\end{picture}}\n"
		}
	}
# Hatched rectangle
	elsif (($q0 eq $orange) or ($q0 eq $brown)) {
		$xtex .= "%Hatched rectangle\n";
		$p1 = 4 * ceil(($xb - $ye) / 4);
		if ($dx >= $dy) {
			$p2 = 4 * floor(($xb - $yb) / 4) + 4;
			$p3 = 4 * floor(($xe - $ye) / 4) + 4;
			$xp = $yb + $p2;
			$yp = $yb;
			$lp = $dy;
			$ip1 = 4;
			$ip2 = 0;
		}
		else {
			$p3 = 4 * floor(($xb - $yb) / 4) + 4;
			$p2 = 4 * floor(($xe - $ye) / 4) + 4;
			$xp = $xb;
			$yp = $xb - $p2;
			$lp = $dx;
			$ip1 = 0;
			$ip2 = -4;
		}
		$p4 = 4 * floor(($xe - $yb) / 4);
		for ($k = $p1; $k <= $p2 - 4; $k += 4) {
			lin($green,$xb,$xb - $k,$ye + $k,$ye);
		}
		$np = ($p3 - $p2) / 4;
		if ($p3 > $p2) {
			if ($lp >= 10 / $ul) {
				$xpu = sp($xp / $ul);
				$ypu = sp($yp / $ul);
				$ip1u = sp($ip1 / $ul);
				$ip2u = sp($ip2 / $ul);
				$lpu = sp($lp / $ul - 0.1);
				$xtex .= "\\multiput(".$xpu.",".$ypu.")(".$ip1u.",".$ip2u;
				$xtex .= "){".$np."}{\\line(1,1){".$lpu."}}\n";
			}
			else {
				for ($k = 0; $k <=$np -1; $k++) {
					$u = $xp + $k * $ip1;
					$v = $yp + $k * $ip2;
					lin($green,$u,$v,$u + $lp,$v + $lp);
				}
			}
		}
		for ($k = $p3; $k <= $p4; $k += 4) {
			lin($green,$yb + $k,$yb,$xe,$xe - $k);
		}
	}
}
#________________________________________________________
# Quadratic B�zier curve
sub qbez {
	my ($x1,$y1,$x2,$y2,$x3,$y3) = @_;
	$xb = $x1;
	$yb = $y1;
	$len = 0.0;
	for ($t = 0.02; $t <= 1.0; $t += 0.02) {
		bound($xb,$yb);
		$s = 1.0 - $t;
		$xe = $s * ($s * $x1 + $t * $x2) + $t * ($s * $x2 + $t * $x3);
		$ye = $s * ($s * $y1 + $t * $y2) + $t * ($s * $y2 + $t * $y3);
		$len += sqrt(($xe - $xb)**2 + ($ye - $yb)**2);
		$xb = $xe;
		$yb = $ye;
	}
	bound($x3,$y3);
	$le = int($bezf * $len * $ul);
	if ($le > $bmax) {$bmax = $le}
	$x1 = sp($x1);
	$y1 = sp($y1);
	$x2 = sp($x2);
	$y2 = sp($y2);
	$x3 = sp($x3);
	$y3 = sp($y3);
	$xtex .= "\\Qbezier[".$le."](".$x1.",".$y1.")(".$x2.",".$y2;
	$xtex .= ")(".$x3.",".$y3.")\n";
}
#________________________________________________________
# Cubic B�zier curve
#
sub cbez {
	my ($x1,$y1,$x2,$y2,$x3,$y3,$x4,$y4) = @_;
	$xb = $x1;
	$yb = $y1;
	$len = 0.0;
	for ($t = 0.02; $t <= 1.0; $t += 0.02) {
		bound($xb,$yb);
		$s = 1.0 - $t;
		$u1 = $s * $x1 + $t * $x2;
		$v1 = $s * $y1 + $t * $y2;
		$u2 = $s * $x2 + $t * $x3;
		$v2 = $s * $y2 + $t * $y3;
		$u3 = $s * $x3 + $t * $x4;
		$v3 = $s * $y3 + $t * $y4;
		$xe = $s * ($s * $u1 + $t * $u2) + $t * ($s * $u2 + $t * $u3);
		$ye = $s * ($s * $v1 + $t * $v2) + $t * ($s * $v2 + $t * $v3);
		$len += sqrt(($xe - $xb)**2 + ($ye - $yb)**2);
		$xb = $xe;
		$yb = $ye;
	}
	bound($x4,$y4);
	$le = int($bezf * $len * $ul);
	if ($le > 0) {
		if ($le > $bmax) {$bmax = $le}
		$x1 = sp($x1);
		$y1 = sp($y1);
		$x2 = sp($x2);
		$y2 = sp($y2);
		$x3 = sp($x3);
		$y3 = sp($y3);
		$x4 = sp($x4);
		$y4 = sp($y4);
		$xtex .= "\\cbezier[".$le."](".$x1.",".$y1.")(".$x2.",";
		$xtex .= $y2.")(".$x3.",".$y3.")(".$x4.",".$y4.")\n";
	}
}
#________________________________________________________
# Quarters of circles
#
sub quart {
	my ($xm,$ym,$r,$al,$bf) = @_;
	if ($al == 0) {
		$dx1 = $r;
		$dy1 = 0;
	}
	elsif ($al == 90) {
		$dx1 = 0;
		$dy1 = $r;
	}
	elsif ($al == 180) {
		$dx1 = (-1) * $r;
		$dy1 = 0;
	}
	elsif ($al == 270) {
		$dx1 = 0;
		$dy1 = (-1) * $r;
	}
	else {
		$al1 = $al * $Pi / 180;
		$dx1 = $r * cos($al1);
		$dy1 = $r * sin($al1);
		$gal = 90 * ceil($al / 90);
		if (($gal == 0) and  ($al != 0)) {
			$gx1 = $r;
			$gy1 = 0;
		}
		elsif (($gal == 90) and ($al != 90)) {
			$gx1 = 0;
			$gy1 = $r;
		}
		elsif (($gal == 180) and ($al != 180)) {
			$gx1 = (-1) * $r;
			$gy1 = 0;
		}
		elsif (($gal == 270) and ($al != 270)) {
			$gx1 = 0;
			$gy1 = (-1) * $r;
		}
		bound($xm + $gx1,$ym + $gy1);
	}
	$dx4 = (-1) * $dy1;
	$dy4 = $dx1;
	$x1 = $xm + $dx1;
	$y1 = $ym + $dy1;
	$x4 = $xm + $dx4;
	$y4 = $ym + $dy4;
	bound($x1,$y1);
	bound($x4,$y4);
	$x2 = $x1 + $l90 * $dx4;
	$y2 = $y1 + $l90 * $dy4;
	$x3 = $x4 + $l90 * $dx1;
	$y3 = $y4 + $l90 * $dy1;
	$le = int($bf * $Pi * $r * $ul/ 2);
	if ($le > $bmax) {$bmax = $le}
	$x1 = sp($x1);
	$y1 = sp($y1);
	$x2 = sp($x2);
	$y2 = sp($y2);
	$x3 = sp($x3);
	$y3 = sp($y3);
	$x4 = sp($x4);
	$y4 = sp($y4);
	$xtex .= "\\cbezier[".$le."](".$x1.",".$y1.")(".$x2;
	$xtex .= ",".$y2.")(".$x3.",".$y3.")(".$x4.",".$y4.")\n";
}
#________________________________________________________
# Short arcs
sub arc {
	my ($xm,$ym,$r,$al,$be,$bf) = @_;
	$al1 = $al * $Pi / 180;
	$be1 = $be * $Pi / 180;
	$dx1 = $r * cos($al1);
	$dy1 = $r * sin($al1);
	$dx4 = $r * cos($be1);
	$dy4 = $r * sin($be1);
	$x1 = $xm + $dx1;
	$y1 = $ym + $dy1;
	$x4 = $xm + $dx4;
	$y4 = $ym + $dy4;
	bound($x1,$y1);
	bound($x4,$y4);
	$gal = 90 * int($al / 90);
	$gbe = 90 * int($be / 90);
	if ($gal != $gbe) {
		if ($gbe == 0) {
			$gx1 = $r;
			$gy1 = 0;
		}
		elsif ($gbe == 90) {
			$gx1 = 0;
			$gy1 = $r;
		}
		elsif ($gbe == 180) {
			$gx1 = (-1) * $r;
			$gy1 = 0;
		}
		elsif ($gbe == 270) {
			$gx1 = 0;
			$gy1 = (-1) * $r;
		}
		bound($xm + $gx1,$ym + $gy1);
	}
	$aux = 4 * (2 * $r - sqrt(($dx1 + $dx4)**2 + ($dy1 + $dy4)**2)) / 3;
	$lam = $aux / sqrt(($dx1 - $dx4)**2 + ($dy1 - $dy4)**2);
	$x2 = $x1 - $lam * $dy1;
	$y2 = $y1 + $lam * $dx1;
	$x3 = $x4 + $lam * $dy4;
	$y3 = $y4 - $lam * $dx4;
	$d1 = abs($al1 - $be1);
	if ($d1 > $Pi / 2) {$d1 = 2 * $Pi - $d1}
	$le = int($bf * $d1 * $r * $ul);
	if ($le > $bmax) {$bmax = $le}
	$x1 = sp($x1);
	$y1 = sp($y1);
	$x2 = sp($x2);
	$y2 = sp($y2);
	$x3 = sp($x3);
	$y3 = sp($y3);
	$x4 = sp($x4);
	$y4 = sp($y4);
	$xtex .= "\\cbezier[".$le."](".$x1.",".$y1.")(".$x2.",";
	$xtex .= $y2.")(".$x3.",".$y3.")(".$x4.",".$y4.")\n";
}
#________________________________________________________
# Bounding box
sub bound {
	my ($xn,$yn) = @_;
	if ($sflag) {
		$xmin = $xn;
		$xmax = $xn;
		$ymin = $yn;
		$ymax = $yn;
		$sflag = 0;
	}
	else {
		if ($xn < $xmin) {$xmin = $xn} elsif ($xmax < $xn) {$xmax = $xn}
		if ($yn < $ymin) {$ymin = $yn} elsif ($ymax < $yn) {$ymax = $yn}
	}
}
#________________________________________________________
# Sprintf
sub sp {
	my $x = shift;
	return sprintf("%.3f",$x) + 0;
}
#________________________________________________________