mathspic script

git-svn-id: svn://tug.org/texlive/trunk@18868 c570f23f-e606-0410-a88d-b1316a301751

3 files changed, 3331 insertions, 0 deletions

diff --git a/Build/source/texk/texlive/linked_scripts/Makefile.am b/Build/source/texk/texlive/linked_scripts/Makefile.am
index 84c7eae8b36..0fc7821a3e7 100644
--- a/Build/source/texk/texlive/linked_scripts/Makefile.am
+++ b/Build/source/texk/texlive/linked_scripts/Makefile.am
@@ -58,6 +58,7 @@ nobase_dist_script_SCRIPTS = \
 	latexmk/latexmk.pl \
 	luaotfload/mkluatexfontdb.lua \
 	listings-ext/listings-ext.sh \
+	mathspic/mathspic.pl \
 	mkjobtexmf/mkjobtexmf.pl \
 	mkgrkindex/mkgrkindex \
 	oberdiek/pdfatfi.pl \
diff --git a/Build/source/texk/texlive/linked_scripts/Makefile.in b/Build/source/texk/texlive/linked_scripts/Makefile.in
index e3afcbaa939..c8e8a72ddc7 100644
--- a/Build/source/texk/texlive/linked_scripts/Makefile.in
+++ b/Build/source/texk/texlive/linked_scripts/Makefile.in
@@ -207,6 +207,7 @@ nobase_dist_script_SCRIPTS = \
 	latexmk/latexmk.pl \
 	luaotfload/mkluatexfontdb.lua \
 	listings-ext/listings-ext.sh \
+	mathspic/mathspic.pl \
 	mkjobtexmf/mkjobtexmf.pl \
 	mkgrkindex/mkgrkindex \
 	oberdiek/pdfatfi.pl \
diff --git a/Build/source/texk/texlive/linked_scripts/mathspic/mathspic.pl b/Build/source/texk/texlive/linked_scripts/mathspic/mathspic.pl
new file mode 100644
index 00000000000..132517b0b17
--- /dev/null
+++ b/Build/source/texk/texlive/linked_scripts/mathspic/mathspic.pl
@@ -0,0 +1,3329 @@
+#!/usr/bin/perl
+#
+#(c) Copyright 2005  Apostolos Syropoulos       & R.W.D. Nickalls
+#                    apostolo@ocean1.ee.duth.gr   dicknickalls@compuserve.com
+#                    apostolo@obelix.ee.duth.gr
+#
+# This program can be redistributed and/or modified under the terms
+# of the LaTeX Project Public License Distributed from CTAN
+# archives in directory macros/latex/base/lppl.txt; either
+# version 1 of the License, or any later version.
+#
+package DummyFH;
+my $index = 0;
+sub TIEHANDLE {
+    my $class = shift;
+    my $self= shift;
+    bless $self, $class;
+}
+sub READLINE {
+    my $self = shift;
+    #shift @$self;
+    if ($index > $#$self) {
+      $index = 0;
+      return undef;
+    }
+    else {
+      return $self->[$index++];
+    }
+}
+
+package main;
+use Math::Trig;
+our $version_number = "1.00 Feb 14, 2005";
+our $commandLineArgs = join(" ", @ARGV);
+our $command  = "";
+our $curr_in_file = "";
+our %PointTable = ();
+our %VarTable = ();
+our %ConstTable = ();
+our $no_errors = 0;
+our $xunits = "1pt";
+our $yunits = "1pt";
+our $units = "pt|pc|in|bp|cm|mm|dd|cc|sp";
+our $defaultsymbol = "\$\\bullet\$";
+our $defaultLFradius = 0;
+use constant PI => atan2(1,1)*4;
+use constant R2D => 180 / PI;
+use constant D2R => PI / 180;
+our $arrowLength = 2;
+our $arrowLengthUnits = "mm";
+our $arrowAngleB = 30;
+our $arrowAngleC = 40; 
+our %DimOfPoint = ();
+our $GlobalDimOfPoints = 0;
+our @Macros = ();
+our $LineThickness = 0.4;
+
+  sub mpp {
+    my $in_line;
+    chomp($in_line = shift);
+    my $LC = shift;
+    my $out_line = $in_line;
+    my $macro_name = "";
+    my @macro_param = ();
+    my $macro_code = "";
+    if ($in_line =~ s/^%def\s*//) {
+      if ($in_line =~ s/^(\w+)\s*//){
+        $macro_name = $1;
+      }
+      else {
+        PrintErrorMessage("No macro name has been found",$LC);
+        return ""
+      }
+      if ($in_line =~ s/^\(\s*//) {
+        # do nothing
+      }
+      else {
+        PrintErrorMessage("No left parenthesis after macro name has been found",$LC);
+        return "";
+      }
+      if ($in_line =~ s/^\)//) {
+        # Macro has no parameters!
+      }
+      else {
+        MACROS: while (1) {
+          if ($in_line =~ s/^(\w+)\s*//) {
+            push (@macro_param, $1);
+          }
+          else {
+            PrintErrorMessage("No macro parameter name has been found",$LC);
+            return "";
+          }
+          if ($in_line =~ s/^,\s*//) {
+            next MACROS;
+          }
+          else {
+            last MACROS;
+          } 
+        }
+        if ($in_line =~ s/^\)//) {
+          # do nothing!
+        }
+        else {
+          PrintErrorMessage("No closing parenthesis after macro parameters",$LC);
+          return "";
+        }
+      }
+      $in_line =~ s/([^%]+)(%.*)/$1/;
+      $macro_code = $in_line;
+      push ( @Macros , { 'macro_name' => $macro_name,
+                         'macro_code' => $macro_code,
+                         'macro_param' => \@macro_param }); 
+      return $out_line;
+    }
+    elsif ($in_line =~ s/^%undef\s*//) {
+      if ($in_line =~ s/^(\w+)//) {
+        my $undef_macro = $1;
+        for(my $i = $#Macros; $i >= 0; $i--) {
+          if ($Macros[$i]->{'macro_name'} eq $undef_macro) {
+           splice(@Macros,$i,1);
+          }   
+        }
+      }
+      return $out_line;
+    }
+    elsif ($in_line =~ s/^\s*%//) {
+      return $out_line;
+    }
+    else {
+      my $comment = $2 if $in_line =~ s/([^%]+)(%.+)/$1/;
+      for(my $i = $#Macros; $i >= 0; $i--) {
+        my $macro_name = $Macros[$i]->{'macro_name'};
+        if ($in_line =~ /&$macro_name/) {       ############################
+          my $num_of_macro_args = @{$Macros[$i]->{'macro_param'}};
+          if ( $num_of_macro_args > 0 ) { 
+           my $pattern = "&$macro_name\\(";
+           foreach my $p ( 1..$num_of_macro_args ) {
+             my $comma = ($p == $num_of_macro_args) ? "\\s*" : "\\s*,\\s*";
+             $pattern .= "\\s*[^\\s\\)]+$comma";
+           }
+           $pattern .= "\\)";
+           while($in_line =~ /\W*&$macro_name\W+/) {
+             if ($in_line =~ /$pattern/) {   
+               my $before = $`;
+                my $after = $';
+                my $match = $&;
+                my $new_code = $Macros[$i]->{'macro_code'};
+                $match =~ s/^&$macro_name\(\s*//;
+               $match =~ s/\)$//; 
+                foreach my $arg ( 0..($num_of_macro_args - 1) ) {
+                 my $old = $Macros[$i]->{'macro_param'}->[$arg];
+                 my $comma = ($arg == ($num_of_macro_args - 1)) ? "" : ",";
+                 $match =~ s/^\s*([^\s,]+)\s*$comma//;
+                 my $new = $1; 
+                 $new_code =~ s/$old/$new/;
+               }
+                $in_line = "$before$new_code$after"; 
+             }
+             else {
+                PrintErrorMessage("Usage of macro &$macro_name does not " .
+                                  "match its definition", $LC); 
+                return "";
+            }
+           }
+         }
+         else {
+          my $replacement = $Macros[$i]->{'macro_code'};
+          $in_line =~ s/&$macro_name/$replacement/g;
+          }
+        }
+      } 
+      return "$in_line$comment";   
+    }
+  }
+
+      sub PrintErrorMessage {
+        my $errormessage = shift;
+        my $error_line   = shift;
+        my ($l,$A);
+        $l = 1+length($command)-length;
+        $A = substr($command,0,$l);
+        $l += 7 +length($error_line);
+
+        for my $fh (STDOUT, LOG) {
+          print $fh "$curr_in_file", "Line $error_line: $A\n";
+          print $fh " " x $l  ,$_,"***Error: $errormessage\n";
+        }
+        if ($comments_on) {  #print to output file file
+          print OUT "%% *** $curr_in_file", "Line $error_line: $A\n";
+          print OUT "%% *** "," " x $l  ,$_,"%% ... Error: $errormessage\n";
+        }
+        $no_errors++;
+      }
+
+      sub PrintWarningMessage {
+        my $warningMessage = shift;
+        my $warning_line   = shift;
+        my ($l,$A);
+        $l = 1+length($command)-length;
+        $A = substr($command,0,$l);
+        $l += 7 +length($warning_line);
+
+        for my $fh (STDOUT, LOG) {
+          print $fh "$curr_in_file", "Line $warning_line: $A\n";
+          print $fh " " x $l  ,$_,"***Warning: $warningMessage\n";
+        }
+        if ($comments_on) {  #print to output file file
+          print OUT "%% *** $curr_in_file", "Line $warning_line: $A\n";
+          print OUT "%% *** "," " x $l  ,$_,"%% ... Warning: $warningMessage\n";
+        }
+      }
+
+      sub PrintFatalError {
+        my $FatalMessage = shift;
+        my $fatal_line   = shift;
+        my ($l,$A);
+        $l = 1+length($command)-length;
+        $A = substr($command,0,$l);
+        $l += 7 +length($fatal_line);
+
+        die "$curr_in_file", "Line $fatal_line: $A\n" .
+            (" " x $l) . $_ . "***Fatal Error: $FatalMessage\n";
+      }
+
+  sub chk_lparen {
+    my $token = $_[0];
+    my $lc    = $_[1];
+    s/\s*//;
+    if (/^[^\(]/) {
+      PrintErrorMessage("Missing ( after $token",$lc);
+    }
+    else {
+      s/^\(\s*//;
+    }
+  }
+
+  sub chk_rparen {
+    my $token = $_[0];
+    my $lc    = $_[1];
+    s/\s*//;
+    if (s/^\)//) {
+      s/\s*//;
+    }
+    else {
+      PrintErrorMessage("Missing ) after $token",$lc);
+    }
+  }
+
+
+  sub chk_lcb {
+    my $token = $_[0];
+    my $lc    = $_[1];
+    s/\s*//;
+    if ($_ !~ /^\{/) {
+      PrintErrorMessage("Missing { after $token",$lc);
+    }
+    else {
+      s/^{\s*//;
+    }
+  }
+
+  sub chk_rcb {
+    my $token = $_[0];
+    my $lc    = $_[1];
+    if ($_ !~ /^\s*\}/) {
+      PrintErrorMessage("Missing } after $token",$lc);
+    }
+    else {
+      s/^\s*}\s*//;
+    }
+  }
+
+  sub chk_lsb {
+     my $token = $_[0];
+     my $lc    = $_[1];
+
+     s/\s*//;
+     if ($_ !~ /^\[/) {
+       PrintErrorMessage("Missing [ after  $token",$lc);
+     }
+     else {
+        s/^\[\s*//;
+     }
+  }
+
+  sub chk_rsb {
+     my $token = $_[0];
+     my $lc    = $_[1];
+
+     s/\s*//;
+     if ($_ !~ /^\]/) {
+       PrintErrorMessage("Missing ] after  $token",$lc);
+     }
+     else {
+       s/^\]\s*//;
+     }
+  }
+
+  sub chk_comma {
+    my $lc  = $_[0];
+
+    s/\s*//;
+    if (/^[^,]/) {
+      PrintErrorMessage("Did not find expected comma",$lc);
+    }
+    else {
+      s/^,\s*//;
+    }
+  }
+
+  sub chk_comment {
+     my $lc = $_[0];
+
+     s/\s*//;
+     if (/^%/) {
+       # do nothing!
+     }
+     elsif (/^[^%]/) {
+       PrintWarningMessage("Trailing text is ignored",$lc);
+     }
+  }
+
+    sub print_headers
+    {
+       my ($sec,$min,$hour,$mday,$mon,$year,$wday,$yday,$isdst) = localtime;
+       $year+=1900;
+       $mon+=1;
+       $now_string = "$year/" . ($mon>9 ? "$mon/" : "0$mon/") .
+                     ($mday>9 ? "$mday    " : "0$mday    ") .
+                     ($hour>9 ? "$hour:" : "0$hour:") .
+                     ($min>9 ? "$min:" : "0$min:") .
+                     ($sec>9 ? "$sec" : "0$sec");
+       print OUT "%* -----------------------------------------------\n";
+       print OUT "%* mathspic (Perl version $version_number)\n";
+       print OUT "%* A filter program for use with PiCTeX\n";
+       print OUT "%* Copyright (c) 2005 A Syropoulos & RWD Nickalls \n";
+       print OUT "%* Command line: $0 $commandLineArgs\n";
+       print OUT "%* Input filename : $source_file\n";
+       print OUT "%* Output filename: $out_file\n";
+       print OUT "%* Date & time: $now_string\n";
+       print OUT "%* -----------------------------------------------\n";
+       #
+       print LOG "----\n";
+       print LOG "$now_string\n";
+       print LOG "mathspic (Perl version $version_number)\n";
+       print LOG "Copyright (c) 2005 A Syropoulos & RWD Nickalls \n";
+       print LOG "Input file  = $source_file\n";
+       print LOG "Output file = $out_file\n";
+       print LOG "Log file    = $log_file\n";
+       print LOG "----\n";
+    }
+
+     sub get_point {
+
+         my ($lc) = $_[0];
+         my ($PointName);
+
+         if (s/^([^\W\d_]\d{0,3})\s*//i) { #point name
+           $PointName = $1;
+           if (!exists($PointTable{lc($PointName)})) {
+             PrintErrorMessage("Undefined point $PointName",$lc);
+             return "_undef_";
+           }
+           else {
+             return lc($PointName);
+           }
+         }
+         else {
+           PrintErrorMessage("Point name expected",$lc);
+           return "_undef_";
+         }
+     }
+
+    sub perpendicular {
+        my ($xP, $yP, $xA, $yA, $xB, $yB) = @_;
+        my ($xF, $yF, $deltax, $deltay, $m1, $m2, $c1, $c2, $factor);
+
+        $deltax = $xA - $xB;
+        return ($xA, $yP) if abs($deltax) < 0.0000001;
+        $deltay = $yA - $yB;
+        return ($xP, $yA) if abs($deltay) < 0.0000001;
+        $m1 = $deltay / $deltax;
+        eval { $m2 = (-1) / $m1;};
+        PrintFatalError("Division by zero",$lc) if $@;
+        $c1 = $yA - $m1 * $xA;
+        $c2 = $yP - $m2 * $xP;
+        eval { $factor = 1 / ($m1 - $m2)};
+        PrintFatalError("Division by zero",$lc) if $@;
+        return (($c2 - $c1) * $factor, ($m1 * $c2 - $m2 * $c1) * $factor);
+    }
+
+    sub Length {
+        my ($xA, $yA, $xB, $yB)=@_;
+        return sqrt(($xB - $xA)**2 + ($yB - $yA)**2);
+    }
+
+    sub triangleArea {
+        my ($xA, $yA, $xB, $yB, $xC, $yC)=@_;
+        my ($lenAB, $lenBC, $lenCA, $s);
+
+        $lenAB = Length($xA,$yA,$xB,$yB);
+        $lenBC = Length($xB,$yB,$xC,$yC);
+        $lenCA = Length($xC,$yC,$xA,$yA);
+        $s = ($lenAB + $lenBC + $lenCA) / 2;
+        return sqrt($s * ($s - $lenAB)*($s - $lenBC)*($s - $lenCA));
+    }
+
+    sub pointOnLine {
+        my ($xA, $yA, $xB, $yB, $dist)=@_;
+        my ($deltax, $deltay, $xPol, $yPol);
+
+        $deltax = $xB - $xA;
+        $deltay = $yB - $yA;
+        $xPol = $xA + ($dist * $deltax / &Length($xA,$yA,$xB,$yB));
+        $yPol = $yA + ($dist * $deltay / &Length($xA,$yA,$xB,$yB));
+        return ($xPol, $yPol);
+    }
+
+
+  
+    sub circumCircleCenter {
+       my ($xA, $yA, $xB, $yB, $xC, $yC, $lc)=@_;
+       my ($deltay12, $deltax12, $xs12, $ys12);
+       my ($deltay23, $deltax23, $xs23, $ys23);
+       my ($xcc, $ycc);
+       my ($m23, $mr23, $c23, $m12, $mr12, $c12);
+       my ($sideA, $sideB, $sideC, $a, $radius);
+
+       if (abs(triangleArea($xA, $yA, $xB, $yB, $xC, $yC)) < 0.0000001)
+       {
+          PrintErrorMessage("Area of triangle is zero!",$lc);
+          return (0,0,0);
+       }
+       $deltay12 = $yB - $yA;
+       $deltax12 = $xB - $xA;
+       $xs12 = $xA + $deltax12 / 2;
+       $ys12 = $yA + $deltay12 / 2;
+       #
+       $deltay23 = $yC - $yB;
+       $deltax23 = $xC - $xB;
+       $xs23 = $xB + $deltax23 / 2;
+       $ys23 = $yB + $deltay23 / 2;
+       #
+       CCXYLINE:{
+       if (abs($deltay12) < 0.0000001)
+       {
+          $xcc = $xs12;
+          if (abs($deltax23) < 0.0000001)
+          {
+             $ycc = $ys23;
+             last CCXYLINE;
+          }
+          else
+          {
+             $m23 = $deltay23 / $deltax23;
+             $mr23 = -1 / $m23;
+             $c23 = $ys23 - $mr23 * $xs23;
+             $ycc = $mr23 * $xs12 + $c23;
+             last CCXYLINE;
+          }
+       }
+       if (abs($deltax12) < 0.0000001)
+       {
+          $ycc = $ys12;
+          if (abs($deltay23) < 0.0000001)
+          {
+             $xcc = $xs23;
+             last CCXYLINE;
+          }
+          else
+          {
+             $m23 = $deltay23 / $deltax23;
+             $mr23 = -1 / $m23;
+             $c23 = $ys23 - $mr23 * $xs23;
+             $xcc = ($ys12 - $c23) / $mr23;
+             last CCXYLINE;
+          }
+       }
+       if (abs($deltay23) < 0.0000001)
+       {
+          $xcc = $xs23;
+          if (abs($deltax12) < 0.0000001)
+          {
+             $ycc = $ys12;
+             last CCXYLINE;
+          }
+          else
+          {
+             $m12 = $deltay12 / $deltax12;
+             $mr12 = -1 / $m12;
+             $c12 = $ys12 - $mr12 * $xs12;
+             $ycc = $mr12 * $xcc + $c12;
+             last CCXYLINE;
+          }
+       }
+       if (abs($deltax23) < 0.0000001)
+       {
+          $ycc = $ys23;
+          if (abs($deltay12) < 0.0000001)
+          {
+             $xcc = $xs12;
+             last CCXYLINE;
+          }
+          else
+          {
+             $m12 = $deltay12 / $deltax12;
+             $mr12 = -1 / $m12;
+             $c12 = $ys12 - $mr12 * $xs12;
+             $xcc = ($ycc - $c12) / $mr12;
+             last CCXYLINE;
+          }
+       }
+       $m12 = $deltay12 / $deltax12;
+       $mr12 = -1 / $m12;
+       $c12 = $ys12 - $mr12 * $xs12;
+       #-----
+       $m23 = $deltay23 / $deltax23;
+       $mr23 = -1 / $m23;
+       $c23 = $ys23 - $mr23 * $xs23;
+       $xcc = ($c23 - $c12) / ($mr12 - $mr23);
+       $ycc = ($c23 * $mr12 - $c12 * $mr23) / ($mr12 - $mr23);
+       }
+       #
+       $sideA = &Length($xA,$yA,$xB,$yB);
+       $sideB = &Length($xB,$yB,$xC,$yC);
+       $sideC = &Length($xC,$yC,$xA,$yA);
+       $a = triangleArea($xA, $yA, $xB, $yB, $xC, $yC);
+       $radius = ($sideA * $sideB * $sideC) / (4 * $a);
+       #
+       return ($xcc, $ycc, $radius);
+    }
+
+    sub ComputeDist {
+       my ($lc) = $_[0];
+       my ($v1, $v2);
+
+       if (s/^((\+|-)?\d+(\.\d+)?([eE](\+|-)?\d+)?)//) #is it a number?
+       {
+           return ($1, 1);
+       }
+       elsif (/^[^\W\d_]\d{0,3}[^\W\d_]\d{0,3}/) #it is a pair of IDs?
+       {
+          s/^([^\W\d_]\d{0,3})//i;
+          $v1 = $1;
+          if (!exists($PointTable{lc($v1)})) {  
+             if (exists($VarTable{lc($v1)})) {
+                return ($VarTable{lc($v1)}, 1);
+             }
+             PrintErrorMessage("Point $v1 has not been defined", $lc);
+             s/^\s*[^\W\d_]\d{0,3}//i;
+             return (0,0);
+          }
+          $v1 = lc($v1);
+          s/^\s*([^\W\d_]\d{0,3})//i;
+          $v2 = $1;
+          if (!exists($PointTable{lc($v2)}))
+          {
+             PrintErrorMessage("Point $v2 has not been defined", $lc);
+             return (0,0);
+          }
+          $v2 = lc($v2);
+          my ($x1,$y1,$pSV1,$pS1) = unpack("d3A*",$PointTable{$v1});
+          my ($x2,$y2,$pSV2,$pS2) = unpack("d3A*",$PointTable{$v2});
+          return (Length($x1,$y1,$x2,$y2), 1);
+       }
+       elsif (s/^([^\W\d_]\d{0,3})//i) # it is a single id
+       {
+         $v1 = $1;
+         if (!exists($VarTable{lc($v1)})) #it isn't a variable
+         {
+           PrintErrorMessage("Variable $v1 has not been defined", $lc);
+           return (0,0);
+         }
+         return ($VarTable{lc($v1)}, 1);
+       }
+       else
+       {
+          PrintErrorMessage("Unexpected token", $lc);
+          return (0,0);
+       }
+     }
+
+     sub intersection4points {
+       my ($x1, $y1, $x2, $y2, $x3, $y3, $x4, $y4) = @_;
+       my ($deltay12, $deltax12, $deltay34, $deltax34);
+       my ($xcc, $ycc, $m34, $c34, $m12, $c12);
+
+       $deltay12 = $y2 - $y1;
+       $deltax12 = $x2 - $x1;
+       #
+       $deltay34 = $y4 - $y3;
+       $deltax34 = $x4 - $x3;
+       I4PXYLINE:{
+          if (abs($deltay12) < 0.0000001)
+          {
+             $ycc = $y1;
+             if (abs($deltax34) < 0.0000001)
+             {
+                $xcc = $x3;
+                last I4PXYLINE;
+             }
+             else
+             {
+                $m34 = $deltay34 / $deltax34;
+                $c34 = $y3 - $m34 * $x3;
+                $xcc = ($ycc - $c34) / $m34;
+                last I4PXYLINE;
+             }
+          }
+          if (abs($deltax12) < 0.0000001)
+          {
+             $xcc = $x1;
+             if (abs($deltay34) < 0.0000001)
+             {
+                $ycc = $y3;
+                last I4PXYLINE;
+             }
+             else
+             {
+                $m34 = $deltay34 / $deltax34;
+                $c34 = $y3 - $m34 * $x3;
+                $ycc = $m34 * $xcc + $c34;
+                last I4PXYLINE;
+             }
+          }
+          if (abs($deltay34) < 0.0000001)
+          {
+             $ycc = $y3;
+             if (abs($deltax12) < 0.0000001)
+             {
+                $xcc = $x1;
+                last I4PXYLINE;
+             }
+             else
+             {
+                $m12 = $deltay12 / $deltax12;
+                $c12 = $y1 - $m12 * $x1;
+                $xcc = ($ycc - $c12) / $m12;
+                last I4PXYLINE;
+             }
+          }
+          if (abs($deltax34) < 0.0000001)
+          {
+             $xcc = $x3;
+             if (abs($deltay12) < 0.0000001)
+             {
+                $ycc = $y1;
+                last I4PXYLINE;
+             }
+             else
+             {
+                $m12 = $deltay12 / $deltax12;
+                $c12 = $y1 - $m12 * $x1;
+                $ycc = $m12 * $xcc + $c12;
+                last I4PXYLINE;
+             }
+          }
+          $m12 = $deltay12 / $deltax12;
+          $c12 = $y1 - $m12 * $x1;
+          $m34 = $deltay34 / $deltax34;
+          $c34 = $y3 - $m34 * $x3;
+          $xcc = ($c34 - $c12) / ($m12 - $m34);
+          $ycc = ($c34 * $m12 - $c12 * $m34) / ($m12 - $m34);
+       }
+       return ($xcc, $ycc);
+     }
+
+    sub IncircleCenter {
+       my ($Ax, $Ay, $Bx, $By, $Cx, $Cy) = @_;
+       my ($sideA, $sideB, $sideC);
+       my ($ba1, $xA1, $yA1, $cb1, $ac1, $xB1, $yB1, $xC1, $yC1, $a, $s, $r);
+
+       #determine the lengths of the sides
+       $sideA = Length($Bx, $By, $Cx, $Cy);
+       $sideB = Length($Cx, $Cy, $Ax, $Ay);
+       $sideC = Length($Ax, $Ay, $Bx, $By);
+       #
+       $ba1 = ($sideC * $sideA) / ($sideB + $sideC);
+       ($xA1, $yA1) = pointOnLine($Bx, $By, $Cx, $Cy, $ba1);
+       $cb1 = ($sideA * $sideB) / ($sideC + $sideA);
+       ($xB1, $yB1) = pointOnLine($Cx, $Cy, $Ax, $Ay, $cb1);
+       $ac1 = ($sideB * $sideC) / ($sideA + $sideB);
+       ($xC1, $yC1) = pointOnLine($Ax, $Ay, $Bx, $By, $ac1);
+       ($xcenter, $ycenter) = &intersection4points($Ax, $Ay, $xA1, $yA1,
+                                                   $Bx, $By, $xB1, $yB1);
+       # get radius
+       $a = &triangleArea($Ax, $Ay, $Bx, $By, $Cx, $Cy);
+       $s = ($sideA + $sideB +$sideC) / 2;
+       $r = $a / $s;
+       return ($xcenter, $ycenter, $r);
+    }
+
+    sub Angle {
+      my ($Ax, $Ay, $Bx, $By, $Cx, $Cy) = @_;
+      my ($RAx, $RAy, $RBx, $RBy, $RCx, $RCy, $deltax, $deltay);
+      my ($lineBA, $lineBC, $lineAC, $k, $kk, $angle);
+      my ($T, $cosT, $sinT) = (0.3, cos(0.3), sin(0.3));
+
+      $RAx = $Ax * $cosT + $Ay * $sinT;
+      $RAy = -$Ax * $sinT + $Ay * $cosT;
+      $RBx = $Bx * $cosT + $By * $sinT;
+      $RBy = -$Bx * $sinT + $By * $cosT;
+      $RCx = $Cx * $cosT + $Cy * $sinT;
+      $RCy = -$Cx * $sinT + $Cy * $cosT;
+      $deltax = $RBx - $RAx;
+      $deltay = $RBy - $RAy;
+      $lineBA = sqrt($deltax*$deltax + $deltay*$deltay);
+      if ($lineBA < 0.0000001)
+      {
+         return -500;
+      }
+      $deltax = $RBx - $RCx;
+      $deltay = $RBy - $RCy;
+      $lineBC = sqrt($deltax*$deltax + $deltay*$deltay);
+      if ($lineBC < 0.0000001)
+      {
+         return -500;
+      }
+      $deltax = $RAx - $RCx;
+      $deltay = $RAy - $RCy;
+      $lineAC = sqrt($deltax*$deltax + $deltay*$deltay);
+      if ($lineAC < 0.0000001)
+      {
+         return -500;
+      }
+      $k = ($lineBA*$lineBA + $lineBC*$lineBC - $lineAC*$lineAC ) /
+           (2 * $lineBA * $lineBC);
+      $k = -1 if $k < -0.99999;
+      $k = 1 if $k > 0.99999;
+      $kk = $k * $k;
+      if (($kk * $kk) == 1)
+      {
+         $angle = PI if $k == -1;
+         $angle = 0 if $k == 1;
+      }
+      else
+      {
+         $angle = (PI / 2) - atan2($k / sqrt(1 - $kk),1);
+      }
+      return $angle * 180 / PI;
+    }
+
+    sub excircle {
+      my ($A, $B, $C, $D, $E) = @_;
+      my ($Ax,$Ay,$Bx,$By,$Dx,$Dy,$Ex,$Ey,$ASVA,$ASA);
+      ($Ax,$Ay,$ASVA,$ASA)=unpack("d3A*",$PointTable{$A});
+      ($Bx,$By,$ASVA,$ASA)=unpack("d3A*",$PointTable{$B});
+      ($Cx,$Cy,$ASVA,$ASA)=unpack("d3A*",$PointTable{$C});
+      ($Dx,$Dy,$ASVA,$ASA)=unpack("d3A*",$PointTable{$D});
+      ($Ex,$Ey,$ASVA,$ASA)=unpack("d3A*",$PointTable{$E});
+      my ($sideA, $sideB, $sideC, $s, $R, $theAdeg, $d);
+      my ($Xmypoint, $Ymypoint, $deltax, $deltay, $mylength, $xc, $yc);
+
+      $sideA = &Length($Bx, $By, $Cx, $Cy);
+      $sideB = &Length($Cx, $Cy, $Ax, $Ay);
+      $sideC = &Length($Ax, $Ay, $Bx, $By);
+      $s = ($sideA + $sideB + $sideC) / 2;
+      $R = triangleArea($Ax, $Ay, $Bx, $By, $Cx, $Cy) /
+           ($s - &Length($Dx, $Dy, $Ex, $Ey));
+      if (($D eq $A && $E eq $B) || ($D eq $B && $E eq $A))
+      {
+        $theAdeg = &Angle($Bx, $By, $Cx, $Cy, $Ax, $Ay);
+        $Xmypoint = $Cx;
+        $Ymypoint = $Cy;
+      }
+      elsif (($D eq $B && $E eq $C) || ($D eq $C && $E eq $B))
+      {
+        $theAdeg = &Angle($Cx, $Cy, $Ax, $Ay, $Bx, $By);
+        $Xmypoint = $Ax;
+        $Ymypoint = $Ay;
+      }
+      elsif (($D eq $C && $E eq $A) || ($D eq $A && $E eq $C))
+      {
+        $theAdeg = &Angle($Ax, $Ay, $Bx, $By, $Cx, $Cy);
+        $Xmypoint = $Bx;
+        $Ymypoint = $By;
+      }
+      else
+      {
+         return (0,0,0);
+      }
+      $d = $R  / sin($theAdeg * PI / 180 / 2);
+      my ($xIn, $yIn, $rin) = &IncircleCenter($Ax, $Ay, $Bx, $By, $Cx, $Cy);
+      $deltax = $xIn - $Xmypoint;
+      $deltay = $yIn - $Ymypoint;
+      $mylength = sqrt($deltax*$deltax + $deltay*$deltay);
+      $xc = $Xmypoint + $d * $deltax / $mylength;
+      $yc = $Ymypoint + $d * $deltay / $mylength;
+      return ($xc, $yc, $R);
+    }
+
+     sub DrawLineOrArrow {      
+       my $draw_Line = shift;
+       my $lc = shift;
+       my $lineLength = -1;
+       my $stacklen = 0;
+       my @PP = ();
+  #     if ($draw_Line != 2) {
+  #       s/\s*//;
+  #       if (s/^\[\s*//) { # optional length specifier 
+  #         $lineLength = expr($lc);
+  #         if ($lineLength <= 0) {
+  #           PrintErrorMessage("length must greater than zero",$lc);
+  #           $lineLength = -1;
+  #         } 
+  #         chk_rsb("optional part",$lc);
+  #       }
+  #     } 
+       chk_lparen("$cmd",$lc);
+       DRAWLINES:while(1) {
+         @PP = () ;
+         while(1) {
+           if (s/^([^\W\d_]\d{0,3})\s*//i) { #point name
+             $P = $1;
+             if (!exists($PointTable{lc($P)})) {
+               PrintErrorMessage("Undefined point $P",$lc);
+             }
+             else {
+               push (@PP,$P);
+             }
+           }
+           else {
+             $stacklen = @PP;
+             if ($draw_Line != 2) {
+               if ($stacklen <= 1) {
+                 PrintErrorMessage("Wrong number of points",$lc);
+               }
+               else {
+                 push(@PP,$lc);
+                 if ($draw_Line == 0) {
+                   drawarrows(@PP); 
+                 }
+                 elsif ($draw_Line == 1) {
+                   drawlines(@PP); 
+                 }
+               }
+             }
+             if (s/^,\s*// and $draw_Line != 2) {
+               next DRAWLINES;
+             }
+             else {
+               last DRAWLINES;
+             }
+           }
+         }
+       }
+       if ($draw_Line == 2) {
+         $stacklen =  @PP;
+         if ($stacklen < 2) {
+           PrintErrorMessage("Wrong number of points",$lc);
+         }
+         elsif ($stacklen % 2 == 0) {
+           PrintErrorMessage("Number of points must be odd",$lc);
+         }
+         else {
+           drawCurve(@PP);        
+         }
+       }
+       chk_rparen("arguments of $cmd",$lc);
+       chk_comment($lc);
+     }
+
+    sub drawarrows {
+      my ($NoArgs);
+      $NoArgs = @_;
+      my ($lc) = $_[$NoArgs-1]; #line number is the last argument
+      my ($NumberOfPoints, $p, $q, $r12, $d12);
+      my ($px,$py,$pSV,$pS, $qx,$qy,$qSV,$qS);
+
+      $NumberOfPoints = $NoArgs - 1;
+      LOOP: for(my $i=0; $i < $NumberOfPoints - 1; $i++)
+      {
+         $p = $_[$i];
+         $q = $_[$i+1];
+         ($px,$py,$pSV,$pS) = unpack("d3A*",$PointTable{lc($p)});
+         ($qx,$qy,$qSV,$qS) = unpack("d3A*",$PointTable{lc($q)});
+         $pSV = $defaultLFradius if $pSV == 0;
+         $qSV = $defaultLFradius if $qSV == 0;
+         $r12 = $pSV + $qSV;
+         $d12 = Length($px,$py,$qx,$qy);
+         if ($d12 <= $r12)
+         {
+            if($d12 == 0)
+            {
+               PrintErrorMessage("points $p and $q are the same", $lc);
+               next LOOP;
+            }
+            PrintWarningMessage("arrow $p$q not drawn: points too close or ".
+                                  "radii too big", $lc);
+            next LOOP;
+         }
+         ($px, $py) = pointOnLine($px, $py, $qx, $qy, $pSV) if $pSV > 0;
+         ($qx, $qy) = pointOnLine($qx, $qy, $px, $py, $qSV) if $qSV > 0;
+         my ($beta, $gamma);
+         $beta  = tan($arrowAngleB * D2R / 2);
+         $gamma = 2 * tan($arrowAngleC * D2R / 2);
+         printf OUT "\\arrow <%.5f%s> [%.5f,%.5f] from %.5f %.5f to %.5f %.5f\n",
+                $arrowLength, $arrowLengthUnits, $beta, $gamma, $px, $py, $qx, $qy;
+      }
+    }
+
+    sub drawlines {
+      my ($NoArgs);
+      $NoArgs = @_;
+      my ($lc) = $_[$NoArgs-1]; #line number is the last argument
+      my ($NumberOfPoints, $p, $q, $r12, $d12);
+      my ($px,$py,$pSV,$pS, $qx,$qy,$qSV,$qS);
+
+      $NumberOfPoints = $NoArgs - 1;
+      LOOP: for(my $i=0; $i < $NumberOfPoints - 1; $i++)
+      {
+         $p = $_[$i];
+         $q = $_[$i+1];
+         ($px,$py,$pSV,$pS) = unpack("d3A*",$PointTable{lc($p)});
+         ($qx,$qy,$qSV,$qS) = unpack("d3A*",$PointTable{lc($q)});
+         $pSV = $defaultLFradius if $pSV == 0;
+         $qSV = $defaultLFradius if $qSV == 0;
+         $r12 = $pSV + $qSV;
+         $d12 = Length($px,$py,$qx,$qy);
+         if ($d12 <= $r12)
+         {
+            if($d12 == 0)
+            {
+               PrintErrorMessage("points $p and $q are the same", $lc);
+               next LOOP;
+            }
+            PrintWarningMessage("line $p$q not drawn: points too close or ".
+                                  "radii too big", $lc);
+            next LOOP;
+         }
+         ($px, $py) = pointOnLine($px, $py, $qx, $qy, $pSV) if $pSV > 0;
+         ($qx, $qy) = pointOnLine($qx, $qy, $px, $py, $qSV) if $qSV > 0;
+         if ($px == $qx || $py == $qy)
+         {
+            printf OUT "\\putrule from %.5f %.5f to %.5f %.5f %%%% %s%s\n",
+                       $px,$py,$qx,$qy,$p,$q;
+         }
+         else
+         {
+            printf OUT "\\plot %.5f %.5f\t%.5f %.5f / %%%% %s%s\n",
+                       $px, $py,$qx,$qy,$p,$q;
+         }
+      }
+    }
+
+    sub drawCurve {
+      my ($NoArgs);
+      $NoArgs = @_;
+      my ($lc) = $_[$NoArgs-1]; #line number is the last argument
+      my ($NumberOfPoints, $p);
+
+      $NumberOfPoints = $NoArgs - 1;
+      print OUT "\\setquadratic\n\\plot\n";
+      for(my $i=0; $i <= $NumberOfPoints; $i++)
+      {
+         $p = $_[$i];
+         my ($px,$py,$pSV,$pS) = unpack("d3A*",$PointTable{lc($p)});
+         printf OUT "\t%0.5f  %0.5f", $px, $py;
+         print OUT (($i == $NumberOfPoints) ? " / %$p\n" : " %$p\n");  
+      }
+      print OUT "\\setlinear\n";
+    }
+
+    sub drawpoints {
+      my ($NumberOfPoints,$p);
+      $NumberOfPoints = @_;
+      my ($px,$py,$pSV,$pS);
+
+      for($i=0; $i < $NumberOfPoints; $i++)
+      {
+         $p = $_[$i];
+         ($px,$py,$pSV,$pS) = unpack("d3A*",$PointTable{lc($p)});
+         if ($pS eq "" and $defaultsymbol =~ /circle|square/) {
+           $pS = $defaultsymbol;
+         }
+         POINTSWITCH: {
+           if ($pS eq "") # no plot symbol specified
+           {
+             printf OUT "\\put {%s} at %.5f %.5f %%%% %s\n",
+                        $defaultsymbol, $px, $py, $p;
+             last POINTSWITCH;
+           }
+           if ($pS eq "circle") # plot symbol is a circle
+           { 
+             my $radius = (defined($DimOfPoint{lc($p)})) ? $DimOfPoint{lc($p)} : 
+                           $GlobalDimOfPoints; 
+             if ($radius > 0) # draw a circle using the current units
+             {
+                if ($radius == 1.5) # use \bigcirc
+                {
+                  printf OUT "\\put {\$\\bigcirc\$} at %.5f %.5f  %%%% %s\n",
+                             $px, $py, $p;
+                }
+                else
+                {
+                  printf OUT "\\circulararc 360 degrees from %.5f %.5f center at %.5f %.5f %%%% %s\n",
+                             $px+$radius, $py, $px, $py, $p;
+                }
+             }
+             else #use \circ symbol
+             {
+               printf OUT "\\put {\$\\circ\$} at %.5f %.5f %%%% %s\n",
+                           $px,$py,$p;
+             }
+             last POINTSWITCH;
+           }
+           if ($pS eq "square")
+           {
+             my $side = (defined($DimOfPoint{lc($p)})) ? $DimOfPoint{lc($p)} : 
+                           $GlobalDimOfPoints;
+             printf OUT "\\put {%s} at %.5f %.5f %%%% %s\n",
+                           drawsquare($side), $px, $py, $p;
+             last POINTSWITCH;
+           }
+           printf OUT "\\put {%s} at %.5f %.5f %%%% %s\n", $pS,$px,$py,$p;
+         }
+      }
+    }
+
+  sub drawAngleArc {
+    my ($P1, $P2, $P3, $radius, $inout, $direction) = @_;
+    my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$P1});
+    my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$P2});
+    my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$P3});
+
+    my $internalAngle = Angle($x1, $y1, $x2, $y2, $x3, $y3);
+    my $externalAngle = 360 - $internalAngle;
+    my ($x, $y) = pointOnLine($x2, $y2, $x1, $y1, $radius);
+    my $code = "";
+    if ($inout eq "internal" and $direction eq "clockwise" ) {
+       $code = sprintf "\\circulararc %.5f degrees from %.5f %.5f center at %.5f %.5f\n",
+               -1 * $internalAngle, $x, $y, $x2, $y2;
+    }
+    elsif ($inout eq "internal" and $direction eq "anticlockwise" ) {
+       $code = sprintf "\\circulararc %.5f degrees from %.5f %.5f center at %.5f %.5f\n",
+               $internalAngle, $x, $y, $x2, $y2;
+    }
+    elsif ($inout eq "external" and $direction eq "clockwise" ) {
+       $code = sprintf "\\circulararc %.5f degrees from %.5f %.5f center at %.5f %.5f\n",
+               -1 * $externalAngle, $x, $y, $x2, $y2;
+    }
+    elsif ($inout eq "external" and $direction eq "anticlockwise" ) {
+       $code = sprintf "\\circulararc %.5f degrees from %.5f %.5f center at %.5f %.5f\n",
+               $externalAngle, $x, $y, $x2, $y2;
+    }
+    return $code;
+  }
+
+  sub drawAngleArrow {
+     my ($P1, $P2, $P3, $radius, $inout, $direction) = @_;
+     my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$P1});
+     my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$P2});
+     my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$P3});
+
+     my $code = drawAngleArc($P1, $P2, $P3, $radius, $inout, $direction);
+
+     my ($xqp, $yqp) = pointOnLine($x2, $y2, $x1, $y1, $radius);
+     my ($deltax, $deltay) = ($x1 - $x2, $y1 - $y2);
+     my $AL;
+
+     if ($xunits =~ /mm/) {
+       $AL = 1;
+     }
+     elsif ($xunits =~ /cm/) {
+       $AL = 0.1;
+     }
+     elsif ($xunits =~ /pt/) {
+       $AL = 2.845;
+     }
+     elsif ($xunits =~ /bp/) {
+       $AL = 2.835;
+     }
+     elsif ($xunits =~ /pc/) {
+       $AL = 0.2371;
+     }
+     elsif ($xunits =~ /in/) {
+       $AL = 0.03937;
+     }
+     elsif ($xunits =~ /dd/) {
+       $AL = 2.659;
+     }
+     elsif ($xunits =~ /cc/) {
+       $AL = 0.2216;
+     }
+     elsif ($xunits =~ /sp/) {
+       $AL = 186467.98;
+     }
+     my $halfAL = $AL / 2;
+     my $d = sqrt($radius * $radius - $halfAL * $halfAL);
+     my $alpha = atan2($d / $halfAL, 1) * R2D;
+     my $beta = 2 * (90 - $alpha);
+     my $thetaqr;
+     if (abs($deltay) < 0.00001) {
+       if ($deltax > 0 ) {$thetaqr = 0 }
+       elsif ($deltax < 0) {$thetaqr = -180}
+     }
+     else {
+       if (abs($deltax) < 0.00001) {
+          $thetaqr = 90;
+       }
+       else {
+         $thetaqr = atan2($deltay / $deltax, 1) * R2D;
+       }
+     }
+     my ($xqr, $yqr) = pointOnLine($x2, $y2, $x3, $y3, $radius);
+     $deltax = $x3 - $x2;
+     $deltay = $y3 - $y2;
+     $alpha = atan2(sqrt($radius * $radius - $halfAL * $halfAL) / $halfAL, 1) /
+              D2R;
+     $beta = 2 * (90 - $alpha);
+     LINE2 : {
+       if (abs($deltax) < 0.00001) {
+         if ($deltay > 0) { $thetaqr = 90 }
+         elsif ($deltay < 0) { $thetaqr = - 90 }
+         last LINE2;
+       }
+       else {
+         $thetaqr = atan2($deltay / $deltax, 1) * R2D;
+       }
+       if (abs($deltay) < 0.00001) {
+         if ($deltax > 0)    { $thetaqr = 0 }
+         elsif ($deltax < 0) { $thetaqr = -180 }
+         last LINE2;
+       }
+       else {
+         $thetaqr = atan2($deltay / $deltax, 1) * R2D;
+       }
+       if ($deltax < 0 and $deltay > 0) { $thetaqr += 180 }
+       elsif ($deltax < 0 and $deltay < 0) { $thetaqr += 180 }
+       elsif ($deltax > 0 and $deltay < 0) { $thetaqr += 360 }
+     }
+     my $xqrleft  = $x2 + $radius * cos(($thetaqr + $beta) * D2R);
+     my $yqrleft  = $y2 + $radius * sin(($thetaqr + $beta) * D2R);
+     my $xqrright = $x2 + $radius * cos(($thetaqr - $beta) * D2R);
+     my $yqrright = $y2 + $radius * sin(($thetaqr - $beta) * D2R);
+     if ($inout eq "internal" and $direction eq "clockwise") {
+       $code .= sprintf "\\arrow <1.5mm> [0.5, 1] from %.5f %.5f to %.5f %.5f\n",
+                $xqrleft, $yqrleft, $xqr, $yqr;
+     }
+     elsif ($inout eq "internal" and $direction eq "anticlockwise") {
+       $code .= sprintf "\\arrow <1.5mm> [0.5, 1] from %.5f %.5f to %.5f %.5f\n",
+                $xqrright, $yqrright, $xqr, $yqr;
+     }
+     elsif ($inout eq "external" and $direction eq "clockwise") {
+       $code .= sprintf "\\arrow <1.5mm> [0.5, 1] from %.5f %.5f to %.5f %.5f\n",
+                $xqrleft, $yqrleft, $xqr, $yqr;
+     }
+     elsif ($inout eq "external" and $direction eq "anticlockwise") {
+       $code .= sprintf "\\arrow <1.5mm> [0.5, 1] from %.5f %.5f to %.5f %.5f\n",
+                $xqrright, $yqrright, $xqr, $yqr;
+     }
+     return $code;
+  }
+
+  sub expr {
+    my $lc = $_[0];
+    my($left,$op,$right);
+
+    $left = term($lc);
+    while ($op = addop()) {
+      $right = term($lc);
+      if ($op eq '+')
+        { $left += $right }
+      else
+        { $left -= $right }
+    }
+    return $left;
+  }
+
+  sub addop {
+    s/^([+-])// && $1;
+  }
+
+  sub term {
+    my $lc = $_[0];
+    my ($left, $op, $right);
+    $left = factor($lc);
+    while ($op = mulop()) {
+      $right = factor($lc);
+      if ($op eq '*')
+        { $left *= $right }
+      elsif ($op =~ /rem/i) {
+        eval {$left %= $right};
+        PrintFatalError("Division by zero", $lc) if $@;
+      }
+      else {
+        eval {$left /= $right};
+        PrintFatalError("Division by zero", $lc) if $@;
+      }
+    }
+    return $left;
+  }
+
+  sub mulop {
+    (s#^([*/])## || s/^(rem)//i) && lc($1);
+  }
+
+  sub factor {
+    my $lc = $_[0];
+    my ($left);
+
+    $left = primitive($lc);
+    if (s/^\*\*//) {
+      $left **= factor($lc);
+    }
+    return $left;
+  }
+
+  sub primitive {
+    my $lc = $_[0];
+    my $val;
+    s/\s*//;
+    if (s/^\(//) {  #is it an expr in parentheses
+      $val = expr($lc);
+      s/^\)// || PrintErrorMessage("Missing right parenthesis", $lc);
+    }
+    elsif (s/^-//) { # is it a negated primitive
+       $val = - primitive();
+    }
+    elsif (s/^\+//) { # is it a negated primitive
+       $val = primitive();
+    }
+    elsif (s/^angledeg//i) {
+       chk_lparen("angledeg",$lc);
+       my $point_1 = get_point($lc);
+       my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point_1});
+       my $point_2 = get_point($lc);
+       my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$point_2});
+       my $point_3 = get_point($lc);
+       my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$point_3});
+       my $d12 = Length($x1, $y1, $x2, $y2);
+       my $d23 = Length($x2, $y2, $x3, $y3);
+       my $d31 = Length($x3, $y3, $x1, $y1);
+       if ( $d12 == 0 ) {
+         PrintErrorMessage("points `$point_1' and `$point_2' are the same", $lc);
+         $val = 0;
+       }
+       elsif ( $d23 == 0 ) {
+         PrintErrorMessage("points `$point_2' and `$point_3' are the same", $lc);
+         $val = 0;
+       }
+       elsif ( $d31 == 0 ) {
+         PrintErrorMessage("points `$point_1' and `$point_3' are the same", $lc);
+         $val = 0;
+       }
+       else {  
+         $val = Angle($x1, $y1, $x2, $y2, $x3, $y3);
+         $val = 0 if $val == -500;
+       }
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^angle//i) {
+       chk_lparen("angle".$lc);
+       my $point_1 = get_point($lc);
+       my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point_1});
+       my $point_2 = get_point($lc);
+       my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$point_2});
+       my $point_3 = get_point($lc);
+       my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$point_3});
+       my $d12 = Length($x1, $y1, $x2, $y2);
+       my $d23 = Length($x2, $y2, $x3, $y3);
+       my $d31 = Length($x3, $y3, $x1, $y1);
+       if ( $d12 == 0 ) {
+         PrintErrorMessage("points `$point_1' and `$point_2' are the same", $lc);
+         $val = 0;
+       }
+       elsif ( $d23 == 0 ) {
+         PrintErrorMessage("points `$point_2' and `$point_3' are the same", $lc);
+         $val = 0;
+       }
+       elsif ( $d31 == 0 ) {
+         PrintErrorMessage("points `$point_1' and `$point_3' are the same", $lc);
+         $val = 0;
+       }
+       else {
+         $val =  Angle($x1, $y1, $x2, $y2, $x3, $y3);
+         if ($val == -500) {
+           $val = 0;
+         }
+         else {
+           $val = D2R * $val;
+         }
+       }
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^area//i) {
+       chk_lparen("angledeg",$lc);
+       my $point_1 = get_point($lc);
+       my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point_1});
+       my $point_2 = get_point($lc);
+       my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$point_2});
+       my $point_3 = get_point($lc);
+       my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$point_3});
+       $val = triangleArea($x1, $y1, $x2, $y2, $x3, $y3);
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^asin//i) {
+       chk_lparen("asin");
+       $val = expr();
+       PrintFatalError("Can't take asin of $val", $lc) if $val < -1 || $val > 1;
+       $val = asin($val);
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^acos//i) {
+       chk_lparen("acos");
+       $val = expr();
+       PrintFatalError("Can't take acos of $val", $lc) if $val < -1 || $val > 1 ;
+       $val = acos($val);
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^atan//i) {
+       chk_lparen("atan");
+       $val = expr();
+       $val = atan($val);
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^cos//i) {
+       chk_lparen("cos");
+       $val = expr();
+       $val = cos($val);
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^directiondeg//i) {
+       chk_lparen("directiondeg",$lc);
+       my $point_1 = get_point($lc);
+       my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point_1});
+       my $point_2 = get_point($lc);
+       my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$point_2});
+       my $x3 = $x1+1;
+       if ( ($y2 - $y1) >= 0) {
+         $val = Angle($x3, $y1, $x1, $y1, $x2, $y2);
+         $val = 0 if $val == -500;
+       }
+       else {
+         $val = 360 - Angle($x3, $y1, $x1, $y1, $x2, $y2);
+         $val = 0 if $val == -500;  
+       }
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^direction//i) {
+       chk_lparen("direction",$lc);
+       my $point_1 = get_point($lc);
+       my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point_1});
+       my $point_2 = get_point($lc);
+       my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$point_2});
+       my $x3 = $x1+1;
+       if ( ($y2 - $y1) >= 0) {
+         $val = Angle($x3, $y1, $x1, $y1, $x2, $y2);
+         $val = 0 if $val == -500;
+         $val = D2R * $val;
+       }
+       else {
+         $val = 360 - Angle($x3, $y1, $x1, $y1, $x2, $y2);
+         $val = 0 if $val == -500;
+         $val = D2R * $val;  
+       }
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^exp//i) {
+       chk_lparen("exp");
+       $val = expr();
+       $val = exp($val);
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^int//i) {
+       chk_lparen("int");
+       $val = expr();
+       $val = int($val);
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^log//i) {
+       chk_lparen("log");
+       $val = expr();
+       PrintFatalError("Can't take log of $val", $lc) if $val <= 0;
+       $val = log($val);
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^sin//i) {
+       chk_lparen("sin");
+       $val = expr();
+       $val = sin($val);
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^sgn//i) {
+       chk_lparen("sgn");
+       $val = expr();
+       if ($val > 0) { 
+         $val = 1;
+       }
+       elsif ($val == 0) {
+         $val = 0;
+       }
+       else {
+         $val = -1;
+       }
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^sqrt//i) {
+       chk_lparen("sqrt");
+       $val = expr();
+       $val = sqrt($val);
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^tan//i) {
+       chk_lparen("tan");
+       $val = expr();
+       $val = sin($val)/cos($val);
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^xcoord//i) {
+       chk_lparen("xcoord");
+       my $point_name = get_point;
+       my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point_name});
+       $val = $x1;
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^ycoord//i) {
+       chk_lparen("ycoord");
+       my $point_name = get_point;
+       my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point_name});
+       $val = $y1;
+       chk_rparen("Missing right parenthesis", $lc);
+    }
+    elsif (s/^_pi_//i) {
+       $val = PI;
+    }
+    elsif (s/^_e_//i) {
+       $val = 2.71828182845905;
+    }
+    elsif (s/^_linethickness_//i) {
+       $val = $LineThickness / $xunits;
+    }
+    else  {
+      my $err_code;
+      ($val,$err_code) = ComputeDist($lc);
+    }
+    s/\s*//;
+    return $val;
+  }
+
+  sub memberOf {
+     my $elem = shift(@_);
+
+     my $found = 0;
+     foreach $item (@_){
+         if ($item eq $elem){
+            $found = 1;
+            last;
+         }
+     }
+     return $found;
+  }
+
+  sub tand {
+    my $d = $_[0];
+    $d = $d * PI / 180;
+    return sin($d)/cos($d);
+  }
+
+  sub get_string {
+    my $string = shift;
+    my $lc = shift;
+
+    $string =~ s/^\s+//;
+    if ($string !~ s/^\"//) {
+      PrintErrorMessage("No starting \" found",$lc);
+      return (1,$string,$string);
+    }
+    my @ch = split //,$string;
+    my @cmd;
+    while (@ch and $ch[0] ne "\"") {
+      if ($ch[0] eq "\\" and (defined $ch[1]) and $ch[1] eq "\"") {
+         shift @ch;
+         push @cmd, $ch[0];
+         shift @ch;
+      }
+      else {
+        push @cmd, $ch[0];
+        shift @ch;
+      }
+    }
+    if (! defined $ch[0]) {
+       PrintErrorMessage("No closing \" found",$lc);
+       return (1,join("",@cmd), join("",@ch))
+    }
+    else {
+      shift @ch;
+      return (0, join("",@cmd), join("",@ch))
+    }
+  }
+
+  sub is_tainted {
+    my $arg = shift;
+    my $nada = substr($arg,0,0);
+    local $@;
+    eval { eval "# $nada"};
+    return length($@) != 0;
+  }
+
+  sub noOfDigits {
+    my $num = $_[0];
+
+    if ($num =~ /^[\+-]?\d+(?!\.)/) {
+      return 0;
+    }
+    elsif ($num =~ /^[\+-]\d+\.(\d+)?/) {
+      return length($1);
+    }
+  }
+
+  sub drawsquare {
+     my $s = $_[0];
+     #$s *= sqrt(2);
+     $s = sprintf "%.5f", $s;
+     my $code = "\\setlength{\\unitlength}{$xunits}%\n";
+     $code .= "\\begin{picture}($s,$s)\\put(0,0)" .
+              "{\\framebox($s,$s){}}\\end{picture}";
+     return $code;
+  }
+
+  sub X2sp {
+     my $LT = shift;
+     my $units = shift;
+
+     if ($units eq "pc") {
+       return $LT * 786432;
+     }
+     elsif ($units eq "pt") {
+       return $LT * 65536;
+     }
+     elsif ($units eq "in") {
+       return $LT * 4736286.72;
+     }
+     elsif ($units eq "bp") {
+       return $LT * 65781.76;
+     }
+     elsif ($units eq "cm") {
+       return $LT * 1864679.811023622;
+     }
+     elsif ($units eq "mm") {
+       return $LT * 186467.981102362;
+     }
+     elsif ($units eq "dd") {
+       return $LT * 70124.086430424;
+     }
+     elsif ($units eq "cc") {
+       return $LT * 841489.037165082;
+     }
+     elsif ($units eq "sp") {
+       return $LT;
+     }
+  }
+
+
+  sub sp2X {
+     my $LT = shift;
+     my $units = shift;
+
+     if ($units eq "pc") {
+       return $LT / 786432;
+     }
+     elsif ($units eq "pt") {
+       return $LT / 65536;
+     }
+     elsif ($units eq "in") {
+       return $LT / 4736286.72;
+     }
+     elsif ($units eq "bp") {
+       return $LT / 65781.76;
+     }
+     elsif ($units eq "cm") {
+       return $LT / 1864679.811023622;
+     }
+     elsif ($units eq "mm") {
+       return $LT / 186467.981102362;
+     }
+     elsif ($units eq "dd") {
+       return $LT / 70124.086430424;
+     }
+     elsif ($units eq "cc") {
+       return $LT / 841489.037165082;
+     }
+     elsif ($units eq "sp") {
+       return $LT;
+     }
+  }
+
+  sub setLineThickness {
+    my $Xunits = shift;
+    my $LT = shift;
+    $Xunits =~ s/^((\+|-)?\d+(\.\d+)?([eE](\+|-)?\d+)?)//;
+    my $xlength = "$1";
+    $Xunits =~ s/^\s*($units)//;
+    my $x_in_units = $1;
+    $LT =~ s/^((\+|-)?\d+(\.\d+)?([eE](\+|-)?\d+)?)//;
+    my $LTlength = "$1";
+    $LT =~ s/^\s*($units)//;
+    my $LT_in_units = $1;
+    $LTlength = X2sp($LTlength,$LT_in_units);
+    $LTlength = sp2X($LTlength,$x_in_units);
+    return $LTlength;
+  }
+
+    sub process_input {
+      my ($INFILE,$currInFile) = @_;
+      my $lc = 0;
+      my $no_output = 0;
+      $curr_in_file = $currInFile;
+      LINE: while(<$INFILE>) {
+        $lc++;
+        chomp($command = $_);
+        s/^\s+//;
+        if (/^beginSkip\s*/i) {
+          $no_output = 1;
+          print OUT "%%$_" if $comments_on;
+          next LINE;
+        }
+        elsif (/^endSkip\s*/i) {
+          if ($no_output == 0) {
+            PrintErrorMessage("endSkip without beginSkip",$lc);
+          }
+          else {
+            $no_output = 0;
+          }
+          print OUT "%%$_" if $comments_on and !$no_output;
+          next LINE;
+       }
+       elsif ($no_output == 1) {
+         next LINE;
+       }
+       else {
+          if (/^[^\\]/) {
+            my $out_line  = mpp($command,$lc) unless /^\\/;  #call macro pre-processor
+            $_ = "$out_line\n";
+          }
+         
+           if (/^\s*%/)
+           {
+              print OUT "$_" if $comments_on;
+           }
+           elsif (s/^\s*(beginloop(?=\W))//i) {
+             s/\s+//;
+             my $times = expr($lc);
+             print OUT "%% BEGINLOOP $times\n" if $comments_on;
+             my @C = ();
+             REPEATCOMMS: while (<$INFILE>) {
+               if (/^\s*endloop/i) {
+                 last REPEATCOMMS;
+               }
+               else {
+                 push @C, $_;
+               }
+             }
+             if (! /^\s*endloop/i) {
+               PrintFatalError("unexpected end of file",$lc);
+             }
+             else {
+               s/^\s*endloop//i;
+               for(my $i=1; $i<=$times; $i++) {
+                 tie *DUMMY, 'DummyFH', \@C;
+                 process_input(DUMMY, $currInFile);
+                 untie *DUMMY;
+               }
+               print OUT "%% ENDLOOP\n" if $comments_on;
+             }
+           }
+           elsif (s/^\s*(ArrowShape(?=\W))//i)
+           {
+               my $cmd = $1;
+               print OUT "%% $cmd$_" if $comments_on;
+              
+                      chk_lparen("$cmd",$lc);
+                      if (s/^default//i) {
+                        $arrowLength = 2;
+                        $arrowLengthUnits = "mm";
+                        $arrowAngleB = 30;
+                        $arrowAngleC = 40;
+                      }
+                      else {
+                        my ($LocalArrowLength, $LocalArrowAngleB ,$LocalArrowAngleC) = (0,0,0);
+                        $LocalArrowLength = expr($lc);
+                        if (s/^\s*($units)//i) {
+                          $arrowLengthUnits = "$1";
+                        }
+                        else {
+                          $xunits =~ /(\d+(\.\d+)?)\s*($units)/;
+                          $LocalArrowLength *= $1;
+                          $arrowLengthUnits = "$3";
+                        }
+                        chk_comma($lc);
+                        $LocalArrowAngleB = expr($lc);
+                        chk_comma($lc);
+                        $LocalArrowAngleC = expr($lc);
+                        $arrowLength = ($LocalArrowLength == 0 ? 2  : $LocalArrowLength);
+                        $arrowLengthUnits = ($LocalArrowLength == 0 ? "mm" : $arrowLengthUnits);
+                        $arrowAngleB = ($LocalArrowAngleB == 0 ? 30 : $LocalArrowAngleB);
+                        $arrowAngleC = ($LocalArrowAngleC == 0 ? 40 : $LocalArrowAngleC);
+                      }
+                      chk_rparen("after $cmd arguments",$lc);
+                      chk_comment("after $cmd command",$lc);
+                      print OUT "%% arrowLength = $arrowLength$arrowLengthUnits, ",
+                                "arrowAngleB = $arrowAngleB ",
+                                "and arrowAngleC = $arrowAngleC\n" if $comments_on;
+
+           }
+           elsif (s/^\s*(const(?=\W))//i)
+           {
+              print OUT "%% $1$_" if $comments_on;
+                  do{
+                    s/\s*//;
+                    PrintErrorMessage("no identifier found after token const",$lc)
+                      if $_ !~ s/^([^\W\d_]\d{0,3})//i;
+                    my $Constname = $1;
+                    my $constname = lc($Constname);
+                    if (exists $ConstTable{$constname}) {
+                      PrintErrorMessage("Redefinition of constant $constname",$lc);
+                    }
+                    s/\s*//; #remove leading white space
+                    PrintErrorMessage("did not find expected = sign",$lc)
+                      if $_ !~ s/^[=]//i;
+                    my $val = expr($lc);
+                    $VarTable{$constname} = $val;
+                    $ConstTable{$constname} = 1;
+                    print OUT "%% $Constname = $val\n" if $comments_on;
+                  }while (s/^,//);
+                  chk_comment($lc);
+                  s/\s*//;
+                  if (/^[^%]/) {
+                    PrintWarningMessage("Trailing text is ignored",$lc);
+                  }
+
+           }
+           elsif (s/^\s*(dasharray(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+                    chk_lparen($cmd,$lc);
+                    my @DashArray = ();
+                    my $dash = "";
+                    my $dashpattern = "";
+                    PATTERN: while (1) {
+                      $dash = sprintf("%.5f", expr($lc));
+                      if (s/^\s*($units)//i) {
+                         push (@DashArray, "$dash$1");
+                      }
+                      else {
+                        PrintErrorMessage("Did not found unit after expression", $lc);
+                      }
+                      s/\s*//;
+                      if (/^[^,]/) {
+                        last PATTERN;
+                      }
+                      else {
+                        s/^,\s*//;
+                      }
+                    }
+                    print OUT "\\setdashpattern <";
+                    while (@DashArray) {
+                      $dashpattern .= shift @DashArray;
+                      $dashpattern .= ",";
+                    }
+                    $dashpattern =~ s/,$//;
+                    print OUT $dashpattern, ">\n";
+                    chk_rparen("arguments of $cmd",$lc);
+                    chk_comment($lc);
+
+           }
+           elsif (s/^\s*(drawAngleArc(?=\W))//i or s/^\s*(drawAngleArrow(?=\W))//i )
+           {
+              my $cmd = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+              
+                    chk_lcb($cmd,$lc);
+                          my ($P1, $P2, $P3);
+                          if (s/^angle(?=\W)//i) {
+                            chk_lparen("token angle of command $cmd",$lc);
+                            $P1 = get_point($lc);
+                            next LINE if $P1 eq "_undef_";
+                            $P2 = get_point($lc);
+                            next LINE if $P2 eq "_undef_";
+                            $P3 = get_point($lc);
+                            next LINE if $P3 eq "_undef_";
+                            my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$P1});
+                            my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$P2});
+                            my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$P3});
+                            my $Angle = Angle($x1, $y1, $x2, $y2, $x3, $y3);
+                            if ($Angle <= 0) {
+                              if ($Angle == 0) {
+                                PrintErrorMessage("Angle is equal to zero",$lc);
+                                next LINE;
+                              }
+                              elsif ($Angle < -400) {
+                                PrintErrorMessage("Something is wrong with the points",$lc);
+                                next LINE;
+                              }
+                            }
+                            chk_rparen("angle part of command $cmd",$lc);
+                          }
+                          else {
+                            PrintErrorMessage("Did not find expected angle part",$lc);
+                            next LINE;
+                          }
+
+                    s/^,\s*// or s/\s*//; #parse optional comma
+                    
+                         my $radius;
+                         if (s/^radius(?=\W)//i) {
+                           chk_lparen("token radius of command $cmd",$lc);
+                           $radius = expr($lc);
+                           chk_rparen("radius part of command $cmd",$lc);
+                         }
+                         else {
+                           PrintErrorMessage("Did not found expected angle part",$lc);
+                           next LINE;
+                         }
+
+                    s/^,\s*// or s/\s*//; #parse optional comma
+                    my $inout = "";
+                    if (s/^(internal(?=\W))//i or s/^(external(?=\W))//i) {
+                      $inout = $1;
+                    }
+                    else {
+                      PrintErrorMessage("Did not find expected 'internal' specifier", $lc);
+                      next LINE;
+                    }
+                    s/^,\s*// or s/\s*//; #parse optional comma
+                    my $direction = "";
+                    if (s/^(clockwise(?=\W))//i or s/^(anticlockwise(?=\W))//i) {
+                      $direction = $1;
+                    }
+                    else {
+                      PrintErrorMessage("Did not find expected 'direction' specifier", $lc);
+                      next LINE;
+                    }
+                    chk_rcb("arguments of $cmd",$lc);
+                    chk_comment($lc);
+                    my $code;
+                    if (lc($cmd) eq "drawanglearc") {
+                      $code = drawAngleArc($P1, $P2, $P3, $radius, $inout, $direction);
+                    }
+                    else {
+                      $code = drawAngleArrow($P1, $P2, $P3, $radius, $inout, $direction);
+                    }
+                    print OUT $code if $code ne "";
+
+           }
+           elsif (s/^\s*(drawArrow(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+              DrawLineOrArrow(0,$lc);
+           }
+           elsif (s/^\s*(drawcircle(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+                    chk_lparen("drawcircle",$lc);
+                    my $Point = get_point($lc);
+                    next LINE if $Point eq "_undef_";
+                    chk_comma($lc);
+                    my $R = expr($lc);
+                    if ($R <= 0) {
+                       PrintErrorMessage("Radius must be greater than zero",$lc);
+                       next LINE;
+                    }
+                    my ($x,$y,$pSV,$pS)=unpack("d3A*",$PointTable{lc($Point)});
+                    printf OUT "\\circulararc 360 degrees from %.5f %.5f center at %.5f %.5f\n",
+                                $x+$R, $y, $x, $y;
+                    chk_rparen("arguments of $cmd",$lc);
+                    chk_comment($lc);
+
+           }
+           elsif (s/^\s*(drawcurve(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+              DrawLineOrArrow(2,$lc);
+           }
+           elsif (s/^\s*(drawcircumcircle(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+                    chk_lparen("drawcircumcircle",$lc);
+                    my $point1 = get_point($lc);
+                    next LINE if $point1 eq "_undef_";
+                    my $point2 = get_point($lc);
+                    next LINE if $point2 eq "_undef_";
+                    my $point3 = get_point($lc);
+                    next LINE if $point3 eq "_undef_";
+                    my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point1});
+                    my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$point2});
+                    my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$point3});
+                    my ($xc, $yc,$r) = circumCircleCenter($x1,$y1,$x2,$y2,$x3,$y3,$lc);
+                    next LINE if $xc == 0 and  $yc == 0 and $r == 0;
+                    print OUT "%% circumcircle center = ($xc,$yc), radius = $r\n" if $comments_on;
+                    printf OUT "\\circulararc 360 degrees from %.5f %.5f center at %.5f %.5f\n",
+                                $xc+$r, $yc, $xc, $yc;
+                    chk_rparen("arguments of $cmd",$lc);
+                    chk_comment($lc);
+
+           }
+           elsif (s/^\s*(drawexcircle(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+                   chk_lparen("drawexcircle",$lc);
+                   my $point1 = get_point($lc);
+                   next LINE if $point1 eq "_undef_";
+                   my $point2 = get_point($lc);
+                   next LINE if $point2 eq "_undef_";
+                   my $point3 = get_point($lc);
+                   next LINE if $point3 eq "_undef_";
+                   my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point1});
+                   my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$point2});
+                   my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$point3});
+                   if (triangleArea($x1, $y1, $x2, $y2, $x3, $y3) < 0.0001) {
+                     PrintErrorMessage("Area of triangle is zero!",$lc);
+                     next LINE;
+                   }
+                   chk_comma($lc);
+                   my $point4 = get_point($lc);
+                   if (!memberOf($point4, $point1, $point2, $point3)) {
+                     PrintErrorMessage("Current point isn't a side point",$lc);
+                     next LINE;
+                   }
+                   next LINE if $point4 eq "_undef_";
+                   my $point5 = get_point($lc);
+                   next LINE if $point5 eq "_undef_";
+                   if (!memberOf($point5, $point1, $point2, $point3)) {
+                     PrintErrorMessage("Current point isn't a side point",$lc);
+                     next LINE;
+                   }
+                   if ($point4 eq $point5) {
+                     PrintErrorMessage("Side points are identical",$lc);
+                     next LINE;
+                   }
+                   chk_rparen("arguments of $cmd",$lc);
+                   my ($xc, $yc, $r) = excircle($point1, $point2, $point3,
+                                                $point4, $point5);
+                   my $R=$r;
+                   if (s/^\s*\[\s*//) {
+                      $R += expr($lc);
+                      if ($R < 0.0001) {
+                        PrintErrorMessage("Radius has become equal to zero!",$lc);
+                        next LINE;
+                      }
+                      chk_rsb($lc);
+                   }
+                   if ($R > (500 / 2.845)) {
+                     PrintErrorMessage("Radius is greater than 175mm!",$lc);
+                     next LINE;
+                   }
+                   print OUT "%% excircle center = ($xc,$yc) radius = $R\n" if $comments_on;
+                   printf OUT "\\circulararc 360 degrees from %.5f %.5f center at %.5f %.5f\n",
+                                $xc+$R, $yc, $xc, $yc;
+                   chk_comment($lc);
+
+           }
+           elsif (s/^\s*(drawincircle(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+                   chk_lparen("drawincircle",$lc);
+                   my $point1 = get_point($lc);
+                   next LINE if $point1 eq "_undef_";
+                   my $point2 = get_point($lc);
+                   next LINE if $point2 eq "_undef_";
+                   my $point3 = get_point($lc);
+                   next LINE if $point3 eq "_undef_";
+                   my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point1});
+                   my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$point2});
+                   my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$point3});
+                   if (triangleArea($x1, $y1, $x2, $y2, $x3, $y3) < 0.0001) {
+                     PrintErrorMessage("Area of triangle is zero!",$lc);
+                     next LINE;
+                   }
+                   my ($xc, $yc, $r) = IncircleCenter($x1,$y1,$x2,$y2,$x3,$y3);
+                   my $R=$r;
+                   if (s/^\s*\[\s*//) {
+                      $R += expr($lc);
+                      if ($R < 0.0001) {
+                        PrintErrorMessage("Radius has become equal to zero!",$lc);
+                        next LINE;
+                      }
+                      chk_rsb($lc);
+                   }
+                   if ($R > (500 / 2.845)) {
+                     PrintErrorMessage("Radius is greater than 175mm!",$lc);
+                     next LINE;
+                   }
+                   print OUT "%% incircle center = ($xc,$yc) radius = $R\n" if $comments_on;
+                   printf OUT "\\circulararc 360 degrees from %.5f %.5f center at %.5f %.5f\n",
+                                $xc+$R, $yc, $xc, $yc;
+                   chk_rparen("arguments of $cmd",$lc);
+                   chk_comment($lc);
+
+           }
+           elsif (s/^\s*(drawline(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+              DrawLineOrArrow(1,$lc);
+           }
+           elsif (s/^\s*(drawthickarrow(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+              print OUT "\\setplotsymbol  ({\\usefont{OT1}{cmr}{m}{n}\\large .})%\n";
+              print OUT "{\\setbox1=\\hbox{\\usefont{OT1}{cmr}{m}{n}\\large .}%\n";
+              print OUT " \\global\\linethickness=0.31\\wd1}%\n";
+              DrawLineOrArrow(0,$lc);
+              print OUT "\\setlength{\\linethickness}{0.4pt}%\n";
+              print OUT "\\setplotsymbol  ({\\usefont{OT1}{cmr}{m}{n}\\tiny .})%\n";
+           }
+           elsif (s/^\s*(drawthickline(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+              print OUT "\\setplotsymbol  ({\\usefont{OT1}{cmr}{m}{n}\\large .})%\n";
+              print OUT "{\\setbox1=\\hbox{\\usefont{OT1}{cmr}{m}{n}\\large .}%\n";
+              print OUT " \\global\\linethickness=0.31\\wd1}%\n";
+              DrawLineOrArrow(1,$lc);
+              print OUT "\\setlength{\\linethickness}{0.4pt}%\n";
+              print OUT "\\setplotsymbol  ({\\usefont{OT1}{cmr}{m}{n}\\tiny .})%\n";
+           }
+           elsif (s/^\s*(drawperpendicular(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+              
+                    chk_lparen($cmd,$lc);
+                    my $A = get_point($lc);
+                    next LINE if $A eq "_undef_";
+                    chk_comma($lc);
+                    my $B = get_point($lc);
+                    next LINE if $A eq "_undef_";
+                    s/\s*//; #ignore white space
+                    my $C = get_point($lc);
+                    next LINE if $A eq "_undef_";
+                    chk_rparen("arguments of $cmd",$lc);
+                    chk_comment($lc);
+                    #
+                    #start actual computation
+                    #
+                    my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$A});
+                    my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$B});
+                    my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$C});
+                    my ($xF, $yF) = perpendicular($x1, $y1, $x2, $y2, $x3, $y3);
+                    printf OUT "\\plot  %.5f %.5f    %.5f %.5f  /\n",
+                           $x1, $y1, $xF, $yF;
+
+           }
+           elsif (s/^\s*(drawpoint(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+                    my ($stacklen);
+                    chk_lparen("$cmd",$lc);
+                    if (/^\)/) {
+                      PrintErrorMessage("There are no point to draw",$lc);
+                      next LINE;
+                    }
+                    my(@PP);
+                    DRAWPOINTS:while(1) {
+                      if (s/^([^\W\d_]\d{0,3})//i) { #point name
+                        $P = $1;
+                        if (!exists($PointTable{lc($P)})) {
+                          PrintErrorMessage("Undefined point $P",$lc);
+                          next DRAWPOINTS;
+                        }
+                        else {
+                          push (@PP,$P);
+                          s/\s*//;
+                        }
+                      }
+                      else {
+                        last DRAWPOINTS;
+                      }
+                    }
+                    drawpoints(@PP);
+                    chk_rparen("arguments of $cmd",$lc);
+                    chk_comment($lc);
+
+           }
+           elsif (s/^\s*(drawRightAngle(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+              
+                   chk_lparen("drawRightAngle",lc);
+                   my $point1 = get_point($lc);
+                   next LINE if $point1 eq "_undef_";
+                   my $point2 = get_point($lc);
+                   next LINE if $point2 eq "_undef_";
+                   my $point3 = get_point($lc);
+                   next LINE if $point3 eq "_undef_";
+                   my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point1});
+                   my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$point2});
+                   my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$point3});
+                   chk_comma($lc);
+                   my $dist = expr($lc);
+                   chk_rparen("arguments of $cmd",$lc);
+                   chk_comment($lc);
+                   #
+                   #actual computation
+                   #
+                   my $lengthAC = Length($x1, $y1, $x3, $y3);
+                   my $angleBAC = Angle($x2, $y2, $x1, $y1, $x3, $y3);
+                   my $line1 = $lengthAC / (1 + tand($angleBAC));
+                   #
+                   # determine coordinates of point Q
+                   #
+                   my ($xQ, $yQ) = pointOnLine($x1, $y1, $x3, $y3, $line1);
+                   #
+                   # determine coordinates of point S
+                   #
+                   my $deltax = $xQ - $x2;
+                   my $deltay = $yQ - $y2;
+                   my $lengthBQ = sqrt($deltax * $deltax + $deltay * $deltay);
+                   my $xS = $x2 + ($dist * sqrt(2) * $deltax / $lengthBQ);
+                   my $yS = $y2 + ($dist * sqrt(2) * $deltay / $lengthBQ);
+                   #
+                   # determine coordinates of points S1 and S2
+                   #
+                   my ($xS1, $yS1) = pointOnLine($x2, $y2, $x3, $y3, $dist);
+                   ($xS2, $yS2) = ($xS, $yS);
+                   #
+                   # emit PiCTeX code
+                   #
+                   printf OUT "\\plot  %.5f %.5f    %.5f %.5f   /\n",
+                          $xS1, $yS1, $xS2, $yS2;
+                   ($xS1, $yS1) = pointOnLine($x2, $y2, $x1, $y1, $dist);
+                   printf OUT "\\plot  %.5f %.5f    %.5f %.5f   /\n",
+                          $xS1, $yS1, $xS2, $yS2;
+
+           }
+           elsif (s/^\s*(drawsquare(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+                 chk_lparen("drawSquare",$lc);
+                 my $p = get_point($lc);
+                 chk_comma($lc);
+                 my $side = expr($lc);
+                 $side = $side - (1.1 * $LineThickness/$xunits); #Suggested by RWDN
+                 my ($x,$y,$pSV,$pS) = unpack("d3A*",$PointTable{$p});
+                 printf OUT "\\put {%s} at %.5f %.5f %%drawsquare\n", drawsquare($side), $x, $y;
+                 chk_rparen("arguments of $cmd",$lc);
+                 chk_comment($lc);
+
+           }
+           elsif (s/^\s*inputfile\*//i)
+           {
+                   chk_lparen("inputfile*",$lc);
+                   my $row_in = "";
+                   if (s/^((\w|-|\.)+)//) {
+                     $row_in = $1;
+                   }
+                   else {
+                     PrintErrorMessage("No input file name found",$lc);
+                     next LINE;
+                   }
+                   if (!(-e $row_in)) {
+                     PrintErrorMessage("File $row_in does not exist",$lc);
+                     next LINE;
+                   }
+                   open(ROW, "$row_in")|| die "Can't open file $row_in\n";
+                   while (defined($in_line=<ROW>)) { print OUT $in_line; }
+                   print OUT "%% ... end of input file <$row_in>\n";
+                   close ROW;
+                   chk_rparen("input file name",$lc);
+                   chk_comment($lc);
+
+
+           }
+           elsif (s/^\s*(inputfile(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+              
+                   chk_lparen("inputfile",$lc);
+                   my $comm_in = "";
+                   if (s/^((\w|-|\.)+)//) {
+                     $comm_in = $1;
+                   }
+                   else {
+                     PrintErrorMessage("No input file name found",$lc);
+                     next LINE;
+                   }
+                   if (!(-e $comm_in)) {
+                     PrintErrorMessage("File $comm_in does not exist",$lc);
+                     next LINE;
+                   }
+                   chk_rparen("input file name",$lc);
+                   my $input_times = 1; #default value
+                   if (s/^\[//) {
+                     $input_times = expr($lc);       
+                     chk_rsb("optional argument",$lc);
+                   }
+                   print OUT "%% ... start of file <$comm_in> loop [$input_times]\n";
+                   for (my $i=0; $i<int($input_times); $i++) {
+                     open(COMM,"$comm_in") or die "Can't open file $comm_in\n";
+                     print OUT "%%% Iteration number: ",$i+1,"\n";
+                     my $old_file_name = $curr_in_file;
+                     process_input(COMM,"File $comm_in, ");
+                     $curr_in_file = $old_file_name;
+                     close COMM;
+                   }
+                   print OUT "%% ... end of file <$comm_in> loop [$input_times]\n";
+                   chk_comment($lc);
+
+           }
+           elsif (s/^\s*(linethickness(?=\W))//i)
+           {
+              my $cmd = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+                   chk_lparen("linethickness", $lc);
+                   if (s/^default//i) {
+                     print OUT "\\linethickness=0.4pt\\Linethickness{0.4pt}%%\n";
+                     print OUT "\\setplotsymbol  ({\\usefont{OT1}{cmr}{m}{n}\\tiny .})%\n";
+                     $LineThickness = setLineThickness($xunits,"0.4pt");
+                   }
+                   else {
+                     my $length = expr($lc);
+                     if (s/^\s*($units)//i) {
+                       my $units = $1;
+                       printf OUT "\\linethickness=%.5f%s\\Linethickness{%.5f%s}%%\n",
+                              $length, $units, $length, $units;
+                       $LineThickness = setLineThickness($xunits,"$length$units");
+                       my $mag;
+                       if ($units eq "pc") {
+                         $mag = $length * 12;
+                       }
+                       elsif ($units eq "in") {
+                         $mag = $length * 72.27;
+                       }
+                       elsif ($units eq "bp") {
+                         $mag = $length * 1.00375;
+                       }
+                       elsif ($units eq "cm") {
+                         $mag = $length * 28.45275;
+                       }
+                       elsif ($units eq "mm") {
+                         $mag = $length * 2.845275;
+                       }
+                       elsif ($units eq "dd") {
+                         $mag = $length * 1.07001;
+                       }
+                       elsif ($units eq "cc") {
+                         $mag = $length * 0.08917;
+                       }
+                       elsif ($units eq "sp") {
+                         $mag = $length * 0.000015259;
+                       }
+                       elsif ($units eq "pt") {
+                         $mag = $length;
+                       }
+                       $mag = 10 * $mag / 1.00278219;
+                       printf OUT "\\font\\CM=cmr10 at %.5fpt%%\n", $mag;
+                       print OUT "\\setplotsymbol  ({\\CM .})%\n";
+                     }
+                     else {
+                       PrintErrorMessage("Did not found expect units part",$lc);
+                     }
+                   }
+                   chk_rparen("linethickness", $lc);
+                   chk_comment($lc);
+
+
+           }
+           elsif (s/^\s*(paper(?=\W))//i)
+           {
+              my ($cmd) = $1;
+              print OUT "%% $cmd$_" if $comments_on;
+                   chk_lcb("paper", $lc);
+                   if (s/^units(?=\W)//i)
+                   {
+                      
+                              chk_lparen("units",$lc);
+                              if(s/^\)//)
+                              {
+                                PrintWarningMessage("Missing value in \"units\"--default is 1pt",
+                                                       $lc);
+                                $xunits = "1pt";          
+                              }
+                              else {
+                                $xunits = expr($lc);
+                                s/\s*//;
+                                if (s/^($units)//i) {
+                                  $xunits .= "$1";
+                                  $LineThickness = setLineThickness($xunits,"0.4pt");
+                                }
+                                elsif(s/^(\w)+//i) {
+                                  PrintErrorMessage("$1 is not a valid mathspic unit",$lc);
+                                  $xunits = "1pt";
+                                }
+                                else {
+                                  PrintErrorMessage("No x-units found",$lc);
+                                  $xunits = "1pt";
+                                }
+                                s/\s*//; #ignore white space
+                                if (s/^,//) {  # there is a comma so expect an y-units
+                                  s/\s*//; #ignore white space
+                                  $yunits = expr($lc);
+                                  s/\s*//; #ignore white space
+                                  if (s/^($units)//i) {
+                                    $yunits .= "$1";
+                                  }
+                                  elsif(s/^(\w)+//i) {
+                                    PrintErrorMessage("$1 is not a valid mathspic unit",$lc);
+                                    $yunits = "1pt";
+                                  }
+                                  else {
+                                    PrintErrorMessage("No y-units found",$lc);
+                                    $yunits = $xunits;
+                                  }
+                                }
+                                else {
+                                  $yunits = $xunits;
+                                }
+                                chk_rparen("units",$lc);
+                              }
+
+                      $nounits = 0;
+                   }
+                   else
+                   {
+                      $nounits = 1;
+                   }
+                   s/^,\s*// or s/\s*//;
+                   if (s/^xrange//i)
+                   {
+                      
+                              chk_lparen("xrange",$lc);
+                              my $ec;
+                              ($xlow,$ec) = ComputeDist($lc);
+                              next LINE if $ec == 0;
+                              chk_comma($lc);
+                              ($xhigh,$ec) = ComputeDist($lc);
+                              next LINE if $ec == 0;
+                              if ($xlow >= $xhigh)
+                              {
+                                 PrintErrorMessage("xlow >= xhigh in xrange",$lc);
+                                 next LINE;
+                              }
+                              chk_rparen("$xhigh",$lc);
+
+                      $noxrange = 0;
+                   }
+                   else
+                   {
+                      $noxrange = 1;
+                   }
+                   s/^,\s*// or s/\s*//;
+                   if (s/^yrange//i)
+                   {
+                      
+                              chk_lparen("yrange",$lc);
+                              my $ec;
+                              ($ylow,$ec) = ComputeDist($lc);
+                              next LINE if $ec == 0;
+                              chk_comma($lc);
+                              ($yhigh,$ec) = ComputeDist($lc);
+                              next LINE if $ec == 0;
+                              if ($ylow >= $yhigh)
+                              {
+                                 PrintErrorMessage("ylow >= yhigh in yrange",$lc);
+                                 next LINE;
+                              }
+                              chk_rparen("$yhigh",$lc);
+
+                      $noyrange = 0;
+                   }
+                   else
+                   {
+                      $noyrange = 1;
+                   }
+                   
+                       if (!$nounits)
+                       {
+                          printf OUT "\\setcoordinatesystem units <%s,%s>\n",
+                                     $xunits,$yunits;
+                       }
+                       if(!$noxrange && !$noyrange)
+                       {
+                          printf OUT "\\setplotarea x from %.5f to %.5f, y from %.5f to %.5f\n",
+                                     $xlow, $xhigh, $ylow, $yhigh;
+
+                       }
+
+                   s/^,\s*// or s/\s*//;
+                   $axis = "";
+                   if (s/^ax[ei]s(?=\W)//i)
+                   {
+                      
+                              chk_lparen("axis",$lc);
+                              while(/^[^\)]/)
+                              {
+                                  if (s/^([lrtbxy]{1}\*?)//i)
+                                  {
+                                     $axis .= $1;
+                                  }
+                                  elsif (s/^([^lrtbxy])//i)
+                                  {
+                                     PrintErrorMessage("Non-valid character \"$1\" in axis()",$lc);
+                                  }
+                                  s/\s*//;
+                              }
+                              chk_rparen("axis(arguments",$lc);
+
+                   }
+                   $axis = uc($axis);
+                   s/^,\s*// or s/\s*//;
+                   if (s/^ticks(?=\W)//i)
+                   {
+                              chk_lparen("ticks",$lc);
+                              my $ec;
+                              ($xticks,$ec) = ComputeDist($lc);
+                              next LINE if $ec == 0;
+                              chk_comma($lc);
+                              ($yticks,$ec) = ComputeDist($lc);
+                              next LINE if $ec == 0;
+                              chk_rparen("ticks(arguments",$lc);
+
+                   }
+                   else
+                   {
+                      $xticks = $yticks = 0;
+                   }
+                   chk_rcb("paper", $lc);
+                   YBRANCH: {
+                      if (index($axis, "Y")>-1)
+                      {
+                         if (index($axis, "Y*")>-1)
+                         {
+                            print OUT "\\axis left shiftedto x=0 / \n";
+                            last YBRANCH;
+                         }
+                         if ($yticks > 0)
+                         {
+                            if (index($axis, "T")>-1 && index($axis, "B")==-1)
+                            {
+                               print OUT "\\axis left shiftedto x=0 ticks numbered from ";
+                               print OUT "$ylow to -$yticks by $yticks\n      from $yticks to ";
+                               print OUT $yhigh-$yticks," by $yticks /\n";
+                            }
+                            elsif (index($axis, "T")==-1 && index($axis, "B")>-1)
+                            {
+                               print OUT "\\axis left shiftedto x=0 ticks numbered from ";
+                               print OUT $ylow+$yticks," to -$yticks by $yticks\n      from ";
+                               print OUT "$yticks to $yhigh by $yticks /\n";
+                            }
+                            elsif (index($axis, "T")>-1 && index($axis, "B")>-1)
+                            {
+                               print OUT "\\axis left shiftedto x=0 ticks numbered from ";
+                               print OUT $ylow+$yticks," to -$yticks by $yticks\n      from ";
+                               print OUT "$yticks to ",$yhigh-$yticks," by $yticks /\n";
+                            }
+                            else
+                            {
+                               print OUT "\\axis left shiftedto x=0 ticks numbered from ";
+                               print OUT "$ylow to -$yticks by $yticks\n      from ";
+                               print OUT "$yticks to $yhigh by $yticks /\n";
+                            }
+                         }
+                         else
+                         {
+                            print OUT "\\axis left shiftedto x=0 /\n";
+                         }
+                      }
+                      }
+                      XBRANCH: { if (index($axis, "X")>-1)
+                      {
+                         if (index($axis, "X*")>-1)
+                         {
+                            print OUT "\\axis bottom shiftedto y=0 /\n";
+                            last XBRANCH;
+                         }
+                         if ($xticks > 0)
+                         {
+                            if (index($axis, "L")>-1 && index($axis, "R")==1)
+                            {
+                               print OUT "\\axis bottom shiftedto y=0 ticks numbered from ";
+                               print OUT $xlow + $xticks," to -$xticks by $xticks\n      from";
+                               print OUT " $xticks to $xhigh by $xticks /\n";
+                            }
+                            elsif (index($axis, "L")==-1 && index($axis, "R")>-1)
+                            {
+                               print OUT "\\axis bottom shiftedto y=0 ticks numbered from ";
+                               print OUT "$xlow to -$xticks by $xticks\n      from ";
+                               print OUT "$xticks to ",$xhigh-$xticks," by $xticks /\n";
+                            }
+                            elsif (index($axis, "L")>-1 && index($axis, "R")>-1)
+                            {
+                               print OUT "\\axis bottom shiftedto y=0 ticks numbered from ";
+                               print OUT $xlow + $xticks," to -$xticks by $xticks\n      from ";
+                               print OUT "$xticks to ",$xhigh - $xticks," by $xticks /\n";
+                            }
+                            else
+                            {
+                               print OUT "\\axis bottom shiftedto y=0 ticks numbered from ";
+                               print OUT "$xlow to -$xticks by  $xticks\n      from ";
+                               print OUT "$xticks to $xhigh by $xticks /\n";
+                            }
+                         }
+                         else
+                         {
+                            print OUT "\\axis bottom shiftedto y=0 /\n";
+                         }
+                      } }
+                      LBRANCH: {if (index($axis, "L")>-1)
+                      {
+                         if (index($axis, "L")>-1)
+                         {
+                            if (index($axis, "L*")>-1)
+                            {
+                               print OUT "\\axis left /\n";
+                               last LBRANCH;
+                            }
+                            if ($yticks > 0)
+                            {
+                               print OUT "\\axis left ticks numbered from ";
+                               print OUT "$ylow to $yhigh by $yticks /\n";
+                            }
+                            else
+                            {
+                               print OUT "\\axis left /\n";
+                            }
+                         }
+                      } }
+                      RBRANCH: { if (index($axis, "R")>-1)
+                      {
+                         if (index($axis, "R*")>-1)
+                         {
+                            print OUT "\\axis right /\n";
+                            last RBRANCH;
+                         }
+                         if ($yticks > 0)
+                         {
+                            print OUT "\\axis right ticks numbered from $ylow to $yhigh by ";
+                            print OUT "$yticks /\n";
+                         }
+                         else
+                         {
+                            print OUT "\\axis right /\n";
+                         }
+                      } }
+                      TBRANCH: { if (index($axis, "T")>-1)
+                      {
+                         if (index($axis, "T*")>-1)
+                         {
+                            print OUT "\\axis top /\n";
+                            last TBRANCH;
+                         }
+                         if ($xticks > 0)
+                         {
+                            print OUT "\\axis top ticks numbered from $xlow to $xhigh by ";
+                            print OUT "$xticks /\n";
+                         }
+                         else
+                         {
+                            print OUT "\\axis top /\n";
+                         }
+                      } }
+                      BBRANCH: { if (index($axis, "B")>-1)
+                      {
+                         if (index($axis, "B*")>-1)
+                         {
+                            print OUT "\\axis bottom /\n";
+                            last BBRANCH;
+                         }
+                         if ($xticks > 0)
+                         {
+                            print OUT "\\axis bottom ticks numbered from $xlow to $xhigh by ";
+                            print OUT "$xticks /\n";
+                         }
+                         else
+                         {
+                            print OUT "\\axis bottom /\n";
+                         }
+                      } }
+
+
+                   chk_comment($lc);
+
+           }
+           elsif (s/^\s*(PointSymbol(?=\W))//i)
+           {
+               my $cmd = $1;
+               print OUT "%% $cmd$_" if $comments_on;
+              
+                      chk_lparen("$cmd",$lc);
+                      if (s/^default//i) #default point symbol
+                      {
+                         $defaultsymbol = "\$\\bullet\$";
+                      }
+                      elsif (s/^(circle|square)//i) {
+                        $defaultsymbol = $1;
+                        chk_lparen($defaultsymbol, $lc);
+                        $GlobalDimOfPoints = expr($lc);
+                        chk_rparen("expression", $lc);          
+                      }
+                      elsif (s/^(((\\,){1}|(\\\)){1}|(\\\s){1}|[^\),\s])+)//) #arbitrary LaTeX point
+                      {
+                         $defaultsymbol = $1;
+                         $defaultsymbol=~ s/\\\)/\)/g;
+                         $defaultsymbol=~ s/\\,/,/g;
+                         $defaultsymbol=~ s/\\ / /g;
+                      }
+                      else
+                      {
+                         PrintErrorMessage("unrecognized point symbol",$lc);
+                      }
+                      if (s/\s*,\s*//) {
+                        $defaultLFradius = expr($lc);
+                      }
+                      chk_rparen("after $cmd arguments",$lc);
+                      chk_comment("after $cmd command",$lc);
+
+           }
+           elsif (s/^\s*point(?=\W)//i)
+           {
+              my ($Point_Line);
+              chomp($Point_Line=$_);
+                   my ($pointStar, $PointName, $origPN);
+                   $pointStar = 0; # default value: he have a point command
+                   $pointStar = 1 if s/^\*//;
+                   chk_lparen("point" . (($pointStar)?"*":""),$lc);
+                   if (s/^([^\W\d_](?![^\W\d_])\d{0,3})//i) {
+                   #
+                   # Note: the regular expression (foo)(?!bar) means that we are
+                   # looking a foo not followed by a bar. Moreover, the regular
+                   # expression [^\W\d_] means that we are looking for letter.
+                   #
+                     $origPN = $1;
+                     $PointName = lc($1);
+                   }
+                   else {
+                     PrintErrorMessage("Invalid point name",$lc);
+                     next LINE;
+                   }
+                   #if ($pointStar and !exists($PointTable{$PointName})) {
+                   #  PrintWarningMessage("Point $origPN has not been defined",$lc);       
+                   #}
+                   if (!$pointStar and exists($PointTable{$PointName})) {
+                      PrintWarningMessage("Point $origPN has been used already",$lc);
+                   }
+                   chk_rparen("point" . (($pointStar)?"*":""). "($origPN",$lc);
+                   chk_lcb("point" . (($pointStar)?"*":""). "($origPN)",$lc);
+                   my ($Px, $Py);
+                          if (s/^perpendicular(?=\W)//i) {
+                                      chk_lparen("perpendicular",$lc);
+                                      my $FirstPoint = &get_point($lc);
+                                      next LINE if $FirstPoint eq "_undef_";
+                                      chk_comma($lc);
+                                      my $SecondPoint = &get_point($lc);
+                                      next LINE if $SecondPoint eq "_undef_";
+                                      my $ThirdPoint = &get_point($lc);
+                                      next LINE if $ThirdPoint eq "_undef_";
+                                      chk_rparen("No closing parenthesis found",$lc);
+                                      my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$FirstPoint});
+                                      my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$SecondPoint});
+                                      my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$ThirdPoint});
+                                      ($Px, $Py) = perpendicular($x1,$y1,$x2,$y2,$x3,$y3);
+
+                          }
+                          elsif (s/^intersection(?=\W)//i) {
+                                       chk_lparen("intersection",$lc);
+                                       my $FirstPoint = get_point($lc);
+                                       next LINE if $FirstPoint eq "_undef_";
+                                       my $SecondPoint = get_point($lc);
+                                       next LINE if $SecondPoint eq "_undef_";
+                                       chk_comma($lc);
+                                       my $ThirdPoint = get_point($lc);
+                                       next LINE if $ThirdPoint eq "_undef_";
+                                       my $ForthPoint = get_point($lc);
+                                       next LINE if $ForthPoint eq "_undef_";
+                                       chk_rparen("No closing parenthesis found",$lc);
+                                       my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$FirstPoint});
+                                       my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$SecondPoint});
+                                       my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$ThirdPoint});
+                                       my ($x4,$y4,$pSV4,$pS4)=unpack("d3A*",$PointTable{$ForthPoint});
+                                       ($Px, $Py) = intersection4points($x1,$y1,$x2,$y2,$x3,$y3,$x4,$y4);
+
+
+                          }
+                          elsif (s/^midpoint(?=\W)//i) {
+                                      chk_lparen("midpoint",$lc);
+                                      my $FirstPoint = &get_point($lc);
+                                      next LINE if $FirstPoint eq "_undef_";
+                                      my $SecondPoint = &get_point($lc);
+                                      next LINE if $SecondPoint eq "_undef_";
+                                      chk_rparen("No closing parenthesis found",$lc);
+                                      my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$FirstPoint});
+                                      my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$SecondPoint});
+                                      $Px = $x1 + ($x2 - $x1)/2;
+                                      $Py = $y1 + ($y2 - $y1)/2;
+
+                          }
+                          elsif (s/^pointonline(?=\W)//i) {
+                                      chk_lparen("pointonline",$lc);
+                                      my $FirstPoint = &get_point($lc);
+                                      next LINE if $FirstPoint eq "_undef_";
+                                      my $SecondPoint = &get_point($lc);
+                                      next LINE if $SecondPoint eq "_undef_";
+                                      chk_comma($lc);
+                                      # now get the distance
+                                      my $distance = expr($lc);
+                                      chk_rparen("No closing parenthesis found",$lc);
+                                      my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$FirstPoint});
+                                      my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$SecondPoint});
+                                      ($Px, $Py) = pointOnLine($x1,$y1,$x2,$y2,$distance);
+
+                          }
+                          elsif (s/^circumcircleCenter(?=\W)//i) {
+                                       chk_lparen("circumCircleCenter",$lc);
+                                       my $FirstPoint = &get_point($lc);
+                                       next LINE if $FirstPoint eq "_undef_";
+                                       my $SecondPoint = &get_point($lc);
+                                       next LINE if $SecondPoint eq "_undef_";
+                                       my $ThirdPoint = &get_point($lc);
+                                       next LINE if $ThirdPoint eq "_undef_";
+                                       chk_rparen("No closing parenthesis found",$lc);
+                                       my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$FirstPoint});
+                                       my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$SecondPoint});
+                                       my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$ThirdPoint});
+                                       ($Px, $Py,$r) = &circumCircleCenter($x1,$y1,$x2,$y2,$x3,$y3,$lc);
+                                       next LINE if $Px == 0 and $Py == 0 and $r == 0;
+
+                          }
+                          elsif (s/^IncircleCenter(?=\W)//i) {
+                                       chk_lparen("IncircleCenter",$lc);
+                                       my $FirstPoint = &get_point($lc);
+                                       next LINE if $FirstPoint eq "_undef_";
+                                       my $SecondPoint = &get_point($lc);
+                                       next LINE if $SecondPoint eq "_undef_";
+                                       my $ThirdPoint = &get_point($lc);
+                                       next LINE if $ThirdPoint eq "_undef_";
+                                       chk_rparen("No closing parenthesis found",$lc);
+                                       my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$FirstPoint});
+                                       my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$SecondPoint});
+                                       my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$ThirdPoint});
+                                       ($Px, $Py, $r) = IncircleCenter($x1,$y1,$x2,$y2,$x3,$y3);
+
+                          }
+                          elsif (s/^ExcircleCenter(?=\W)//i) {
+                                       chk_lparen("ExcircleCenter",$lc);
+                                       my $PointA = get_point($lc);
+                                       next LINE if $PointA eq "_undef_";
+                                       my $PointB = get_point($lc);
+                                       next LINE if $PointB eq "_undef_";
+                                       my $PointC = get_point($lc);
+                                       next LINE if $PointC eq "_undef_";
+                                       chk_comma($lc);
+                                       my $PointD = &get_point($lc);
+                                       next LINE if $PointD eq "_undef_";
+                                       if (!memberOf($PointD, $PointA, $PointB, $PointC)) {
+                                         PrintErrorMessage("Current point isn't a side point",$lc);
+                                         next LINE;
+                                       }
+                                       my $PointE = get_point($lc);
+                                       next LINE if $PointE eq "_undef_";
+                                       if (!memberOf($PointE, $PointA, $PointB, $PointC)) {
+                                         PrintErrorMessage("Current point isn't a side point",$lc);
+                                         next LINE;
+                                       }
+                                       if ($PointD eq $PointE) {
+                                         PrintErrorMessage("Side points are identical",$lc);
+                                         next LINE;
+                                       }
+                                       ($Px, $Py, $r) = excircle($PointA, $PointB, $PointC,
+                                                                  $PointD, $PointE);
+                                       chk_rparen("after coordinates part",$lc);
+
+                          }
+                          elsif (/^[^\W\d_]\d{0,3}\s*[^,\w]/) {
+                            m/^([^\W\d_]\d{0,3})\s*/i;
+                            if (exists($PointTable{lc($1)})) {         
+                              my $Tcoord = get_point($lc);
+                              my ($x,$y,$pSV,$pS)=unpack("d3A*",$PointTable{$Tcoord});
+                              $Px = $x;
+                              $Py = $y;
+                            }
+                            else {
+                              $Px = expr();
+                              chk_comma($lc);
+                              $Py = expr();
+                            }
+                          }
+                          elsif (/[^\W\d_]\d{0,3}\s*,\s*shift|polar|rotate|vector/i) { #a point?
+                            s/^([^\W\d_]\d{0,3})//i;
+                            my $PointA = $1;
+                            if (exists($PointTable{lc($PointA)})) {
+                              s/\s*//;
+                              if (s/^,//) {
+                                s/\s*//;
+                                if (s/^shift(?=\W)//i) {
+                                  
+                                            chk_lparen("shift",$lc);
+                                            my $dist1 = expr($lc);
+                                            chk_comma($lc);
+                                            my $dist2 = expr($lc);
+                                            my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{lc($PointA)});
+                                            $Px = $x1 + $dist1;
+                                            $Py = $y1 + $dist2;
+                                            chk_rparen("shift part",$lc);
+
+                                }
+                                elsif (s/^polar(?=\W)//i) {
+                                            chk_lparen("polar",$lc);
+                                            my ($R1, $Theta1);
+                                            $R1 = expr($lc);
+                                            chk_comma($lc);
+                                            $Theta1 = expr($lc);
+                                            my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{lc($PointA)});
+                                            s/\s*//;
+                                            if (s/^rad(?=\W)//i) {
+                                               # do nothing!
+                                            }
+                                            elsif (s/^deg(?=\W)//i) {
+                                              $Theta1 = $Theta1 * PI / 180;
+                                            }
+                                            else {
+                                              #$Theta1 = $Theta1 * PI / 180; 
+                                            }
+                                            $Px = $x1 + $R1 * cos($Theta1);
+                                            $Py = $y1 + $R1 * sin($Theta1);
+                                            chk_rparen("after polar part",$lc);
+
+                                }
+                                elsif (s/^rotate(?=\W)//i) {
+                                            chk_lparen("rotate",$lc);
+                                            my $Q = lc($PointA);
+                                            my $P = get_point($lc);
+                                            next LINE if $P eq "_undef_";
+                                            chk_comma($lc);
+                                            my $Theta1 = expr($lc);
+                                            my ($xP,$yP,$pSV1,$pS1)=unpack("d3A*",$PointTable{$P});
+                                            my ($xQ,$yQ,$pSV2,$pS2)=unpack("d3A*",$PointTable{$Q});
+                                            s/\s*//;
+                                            if (s/^rad(?=\W)//i)
+                                            {
+                                               # do nothing!
+                                            }
+                                            elsif (s/^deg(?=\W)//i)
+                                            {
+                                                $Theta1 = $Theta1 * PI / 180;
+                                            }
+                                            else
+                                            {
+                                                $Theta1 = $Theta1 * PI / 180;
+                                            }
+                                            # shift origin to P
+                                            $xQ -= $xP;
+                                            $yQ -= $yP;
+                                            # do the rotation
+                                            $Px = $xQ * cos($Theta1) - $yQ * sin($Theta1);
+                                            $Py = $xQ * sin($Theta1) + $yQ * cos($Theta1);
+                                            # return origin back to original origin
+                                            $Px += $xP;
+                                            $Py += $yP;
+                                            chk_rparen("after rotate part",$lc);
+
+                                }
+                                elsif (s/^vector(?=\W)//i) {
+                                            chk_lparen("vector",$lc);
+                                            my ($x0,$y0,$pSV0,$pS0) = unpack("d3A*",$PointTable{lc($PointA)});
+                                            my $P = get_point($lc);
+                                            my $Q = get_point($lc);
+                                            my ($x1,$y1,$pSV1,$pS1) = unpack("d3A*",$PointTable{$P});
+                                            my ($x2,$y2,$pSV2,$pS2) = unpack("d3A*",$PointTable{$Q});
+                                            $Px = $x0 + $x2 - $x1;
+                                            $Py = $y0 + $y2 - $y1;
+                                            chk_rparen("vector part",$lc);
+
+
+                                }
+                                else {
+                                  PrintErrorMessage("unexpected token",$lc);
+                                  next LINE;
+                                }
+                              }
+                              else {
+                                my ($xA,$yA,$pSVA,$pSA)=unpack("d3A*",$PointTable{lc($PointA)});
+                                $Px = $xA;
+                                $Py = $yA;
+                              }
+                            }
+                            else {
+                              PrintErrorMessage("Undefined point $PointA",$lc);
+                              next LINE;
+                            }
+                          }
+                          else {
+                            $Px = expr();
+                            chk_comma($lc);
+                            $Py = expr();
+                          }
+
+                   chk_rcb("coordinates part",$lc);
+                   my $sv = $defaultsymbol;
+                   my $sh = $defaultLFradius;
+                   my $side_or_radius = undef;
+                   if (s/^\[\s*//) { # the user has opted to specify the optional part
+                                if (/^(symbol|radius|side)\s*/i) {
+                                  my @previous_options = ();
+                                  my $number_of_options = 1;
+                                  my $symbol_set = 0;        
+                                  while (s/^(symbol|radius)\s*//i and $number_of_options <= 2) {
+                                    my $option = lc($1);
+                                    if (s/^=\s*//) {
+                                      if (memberOf($option,@previous_options)) {
+                                        PrintErrorMessage("Option \"$option\" has been already defined", $lc);
+                                        my $dummy = expr($lc);
+                                      } 
+                                      elsif ($option eq "radius") {
+                                        $sh = expr($lc);
+                                        $sv = $defaultsymbol if ! $symbol_set;
+                                      }
+                                      elsif ($option eq "symbol") {
+                                        if (s/^circle\s*//i) {
+                                          $sv = "circle";
+                                          chk_lparen("after token circle",$lc);
+                                          $side_or_radius = expr($lc);
+                                          chk_rparen("expression",$lc);
+                                        }
+                                        elsif (s/^square\s*//i) {
+                                          $sv = "square";
+                                          chk_lparen("after token square",$lc);
+                                          $side_or_radius = expr($lc);
+                                          chk_rparen("expression",$lc);
+                                        }
+                                        elsif (s/^(((\\\]){1}|(\\,){1}|(\\\s){1}|[^\],\s])+)//) {
+                                          $sv = $1;
+                                          $sv =~ s/\\\]/\]/g;
+                                          $sv =~ s/\\,/,/g;
+                                          $sv =~ s/\\ / /g;
+                                          s/\s*//;
+                                        }
+                                        $symbol_set = 1;
+                                      }
+                                    }
+                                    else {
+                                      PrintErrorMessage("unexpected token", $lc);
+                                      next LINE;
+                                    }
+                                    $number_of_options++;
+                                    push (@previous_options, $option);
+                                    s/^,\s*//;
+                                  }
+                                }
+                                else {
+                                  PrintErrorMessage("unexpected token", $lc);
+                                  next LINE;
+                                }
+
+                      chk_rsb("optional part",$lc);
+                   }
+                   # to avoid truncation problems introduced by the pack function, we
+                   # round each number up to five decimal digits
+                   $Px = sprintf("%.5f", $Px);
+                   $Py = sprintf("%.5f", $Py);
+                   print OUT "%% point$Point_Line \t$origPN = ($Px, $Py)\n" if $comments_on;
+                   chk_comment($lc);
+                   $PointTable{$PointName} = pack("d3A*",$Px,$Py,$sh,$sv);
+                   if (defined($side_or_radius)) {
+                     $DimOfPoint{$PointName} = $side_or_radius;
+                   } 
+
+           }
+           elsif (/^\s*setPointNumber(?=\W)/i)
+           {
+              PrintWarningMessage("Command setPointNumber is ignored",$lc);
+              next LINE;
+           }
+           elsif (s/^\s*(showAngle(?=\W))//i)
+           {
+                     chk_lparen("showangle",$lc);
+                     my $point_1 = get_point($lc);
+                     my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point_1});
+                     my $point_2 = get_point($lc);
+                     my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$point_2});
+                     my $point_3 = get_point($lc);
+                     my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$point_3});
+                     my $angle = Angle($x1, $y1, $x2, $y2, $x3, $y3);
+                     $angle = 0 if $angle == -500;
+                     printf OUT "%%%% angle(%s%s%s) = %.5f deg ( %.5f rad)\n", $point_1,
+                            $point_2, $point_3, $angle, $angle*D2R;
+                     chk_rparen("Missing right parenthesis", $lc);
+
+           }
+           elsif (s/^\s*(showArea(?=\W))//i)
+           {
+                     chk_lparen("showarea",$lc);
+                     my $point_1 = get_point($lc);
+                     my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point_1});
+                     my $point_2 = get_point($lc);
+                     my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$point_2});
+                     my $point_3 = get_point($lc);
+                     my ($x3,$y3,$pSV3,$pS3)=unpack("d3A*",$PointTable{$point_3});
+                     print OUT "%% area($point_1$point_2$point_3) = ",
+                           triangleArea($x1, $y1, $x2, $y2, $x3, $y3), "\n";
+                     chk_rparen("Missing right parenthesis", $lc);
+
+           }
+           elsif (s/^\s*(showLength(?=\W))//i)
+           {
+                     chk_lparen("showlength",$lc);
+                     my $point_1 = get_point($lc);
+                     my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{$point_1});
+                     my $point_2 = get_point($lc);
+                     my ($x2,$y2,$pSV2,$pS2)=unpack("d3A*",$PointTable{$point_2});
+                     print OUT "%% length($point_1$point_2) = ",
+                           Length($x1, $y1, $x2, $y2), "\n";
+                     chk_rparen("Missing right parenthesis", $lc);
+
+
+           }
+           elsif (/^\s*showPoints(?=\W)/i)
+           {
+              print OUT "%%-------------------------------------------------\n";
+              print OUT "%%            L I S T  O F  P O I N T S            \n";
+              print OUT "%%-------------------------------------------------\n";
+              foreach my $p (keys(%PointTable)) {
+                my ($x, $y, $pSV, $pS) = unpack("d3A*",$PointTable{$p});
+                printf OUT "%%%%\t%s\t= ( %.5f, %.5f ), LF-radius = %.5f, symbol = %s\n",
+                        $p, $x, $y, $pSV, $pS;
+              }
+              print OUT "%%-------------------------------------------------\n";
+              print OUT "%%      E N D  O F  L I S T  O F  P O I N T S      \n";
+              print OUT "%%-------------------------------------------------\n";
+              next LINE;
+           }
+           elsif (/^\s*showVariables(?=\W)/i)
+           {
+              print OUT "%%-------------------------------------------------\n";
+              print OUT "%%       L I S T  O F  V A R I A B L E S           \n";
+              print OUT "%%-------------------------------------------------\n";
+              foreach my $var (keys(%VarTable)) {
+                print OUT "%%\t", $var, "\t=\t", $VarTable{$var}, "\n";
+              }
+              print OUT "%%-------------------------------------------------\n";
+              print OUT "%%   E N D  O F  L I S T  O F  V A R I A B L E S   \n";
+              print OUT "%%-------------------------------------------------\n";
+              next LINE;
+           }
+           elsif (s/^\s*(system(?=\W))//i)
+           {
+              print OUT "%% $1$_" if $comments_on;
+              
+                   chk_lparen("$cmd",$lc);
+                   my ($error, $command, $rest) = get_string($_);
+                   next LINE if $error == 1;
+                   $_ = $rest;
+                   if (! is_tainted($command)) {
+                      system($command);
+                   }
+                   else {
+                      PrintErrorMessage("String \"$command\" has tainted data", $lc);
+                      next LINE;
+                   }
+                   chk_rparen("after $cmd arguments",$lc);
+                   chk_comment("after $cmd command",$lc);
+
+           }
+           elsif (s/^\s*(text(?=\W))//i)
+           {
+              print OUT "%% $1$_" if $comments_on;
+              
+                      chk_lparen("text",$lc);
+                      my ($level,$text)=(1,"");
+                      TEXTLOOP: while (1)
+                      {
+                        $level++ if /^\(/;
+                        $level-- if /^\)/;
+                        s/^(.)//;
+                        last TEXTLOOP if $level==0;
+                        $text .= $1;
+                      }
+                      chk_lcb("text part",$lc);
+                      my ($Px, $Py,$dummy,$pos);
+                      $pos="";
+                      s/\s*//;
+                      
+                        if (/^[^\W\d_]\d{0,3}\s*[^,\w]/) {
+                          my $Tcoord = get_point($lc);
+                          my ($x,$y,$pSV,$pS)=unpack("d3A*",$PointTable{$Tcoord});
+                          $Px = $x;
+                          $Py = $y;
+                        }
+                        elsif (/[^\W\d_]\d{0,3}\s*,\s*shift|polar/i) {
+                          s/^([^\W\d_]\d{0,3})//i;
+                          my $PointA = $1;
+                          if (exists($PointTable{lc($PointA)})) {
+                            s/\s*//;
+                            if (s/^,//) {
+                              s/\s*//;
+                              if (s/^shift(?=\W)//i) {
+                                
+                                          chk_lparen("shift",$lc);
+                                          my $dist1 = expr($lc);
+                                          chk_comma($lc);
+                                          my $dist2 = expr($lc);
+                                          my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{lc($PointA)});
+                                          $Px = $x1 + $dist1;
+                                          $Py = $y1 + $dist2;
+                                          chk_rparen("shift part",$lc);
+
+                              }
+                              elsif (s/^polar(?=\W)//i) {
+                                          chk_lparen("polar",$lc);
+                                          my ($R1, $Theta1);
+                                          $R1 = expr($lc);
+                                          chk_comma($lc);
+                                          $Theta1 = expr($lc);
+                                          my ($x1,$y1,$pSV1,$pS1)=unpack("d3A*",$PointTable{lc($PointA)});
+                                          s/\s*//;
+                                          if (s/^rad(?=\W)//i) {
+                                             # do nothing!
+                                          }
+                                          elsif (s/^deg(?=\W)//i) {
+                                            $Theta1 = $Theta1 * PI / 180;
+                                          }
+                                          else {
+                                            #$Theta1 = $Theta1 * PI / 180; 
+                                          }
+                                          $Px = $x1 + $R1 * cos($Theta1);
+                                          $Py = $y1 + $R1 * sin($Theta1);
+                                          chk_rparen("after polar part",$lc);
+
+                              }
+                            }
+                          }
+                          else {
+                            PrintErrorMessage("undefined point/var",$lc);
+                            next LINE;
+                          }
+                        }
+                        else {
+                          $Px = expr();
+                          chk_comma($lc);
+                          $Py = expr();
+                        }
+
+                      chk_rcb("coordinates part of text command",$lc);
+                      if (s/^\[//)
+                      {
+                        s/\s*//;
+                        
+                          if (s/^(\w{1})\s*//) {
+                            $pos .= $1;
+                            if (memberOf($pos, "l", "r")) {
+                              if (s/^(\w{1})\s*//) {
+                                my $np = $1;
+                                if (memberOf($np, "t", "b", "B")) {
+                                  $pos .= $np;
+                                }
+                                else {
+                                  if (memberOf($np, "l", "r")) {
+                                     PrintErrorMessage("$np can't follow 'l' or 'r'", $lc);
+                                  }
+                                  else {
+                                    PrintErrorMessage("$np is not a valid positioning option", $lc);
+                                  }
+                                  next LINE;
+                                }
+                              }
+                            }
+                            elsif (memberOf($pos, "t", "b", "B")) {
+                              if (s/^(\w{1})\s*//) {
+                                my $np = $1;
+                                if (memberOf($np, "l", "r")) {
+                                  $pos .= $np;
+                                }
+                                else {
+                                  if (memberOf($np, "t", "b", "B")) {
+                                     PrintErrorMessage("$np can't follow 't', 'b', or 'B'", $lc);
+                                  }
+                                  else {
+                                    PrintErrorMessage("$np is not a valid positioning option", $lc);
+                                  }
+                                  next LINE;
+                                }
+                              }
+                            }
+                            else {
+                              PrintErrorMessage("$pos is not a valid positioning option", $lc);
+                              next LINE;
+                            }
+                          }
+                          else {
+                            PrintErrorMessage("illegal token in optional part of text command",$lc);
+                            next LINE;
+                          }
+
+                        s/\s*//;
+                        chk_rsb("optional part of text command",$lc);
+                      }
+                      chk_comment($lc);
+                      if ($pos eq "")
+                      {
+                         printf OUT "\\put {%s} at %f %f\n", $text, $Px, $Py;
+                      }
+                      else
+                      {
+                         printf OUT "\\put {%s} [%s] at %f %f\n", $text, $pos, $Px, $Py;
+                      }
+
+           }
+           elsif (s/^\s*(var(?=\W))//i)
+           {
+              print OUT "%% $1$_" if $comments_on;
+                  do{
+                    s/\s*//;
+                    PrintErrorMessage("no identifier found after token var",$lc)
+                      if $_ !~ s/^([^\W\d_]\d{0,3})//i;
+                    my $Varname = $1;
+                    my $varname = lc($Varname);
+                    if (exists $ConstTable{$varname}) {
+                      PrintErrorMessage("Redefinition of constant $varname",$lc);
+                    }
+                    s/\s*//; #remove leading white space
+                    PrintErrorMessage("did not find expected = sign",$lc)
+                      if $_ !~ s/^[=]//i;
+                    my $val = expr($lc);
+                    $VarTable{$varname} = $val;
+                    print OUT "%% $Varname = $val\n" if $comments_on;
+                  }while (s/^,//);
+                  chk_comment($lc);
+                  s/\s*//;
+                  if (/^[^%]/) {
+                    PrintWarningMessage("Trailing text is ignored",$lc);
+                  }
+           }
+           elsif (/^\s*\\(.+)/)
+           {
+              my $line = $1;
+              if ($line =~ /^\s+(.+)/)
+              {
+                 print OUT " $line\n";
+              }
+              else
+              {
+                print OUT "\\$line\n";
+              }
+              next LINE;
+           }
+           elsif (0==length) #empty line
+           {
+              next LINE;
+           }
+           else {
+             PrintErrorMessage("command not recognized",$lc);
+             next LINE;
+           }
+
+       }
+     }
+   }
+
+
+our $alarm="";
+our $comments_on=1;
+our $out_file="default";
+our $argc=@ARGV;
+if ($argc == 0 || $argc > 5 ){ # no command line arguments or more than 4
+                               # arguments
+  die "\nmathspic version $version_number\n" .
+      "Usage: mathspic [-h] [-b] [-c] [-o <out file>] <in file>\n\n";
+}
+else {
+     our $file = "";
+     SWITCHES:
+     while($_ = $ARGV[0]) {
+       shift;     
+       if (/^-h$/) {
+         die "\nThis is mathspic version $version_number\n" .
+             "Type \"man mathspic\" for detailed help\n".
+             "Usage:\tmathspic  [-h] [-b] [-c] [-o <out file>] <in file>\n" . 
+             "\twhere,\n" . 
+             "\t[-b]\tenables bell sound if error exists\n" .
+             "\t[-c]\tdisables comments in ouput file\n" .
+             "\t[-h]\tgives this help listing\n" .
+             "\t[-o]\tcreates specified output file\n\n"; 
+       }
+       elsif (/^-b$/) {
+         $alarm = chr(7);
+       }
+       elsif (/^-c$/) {
+         $comments_on = 0;
+       }
+       elsif (/^-o$/) {
+         die "No output file specified!\n" if !@ARGV;
+         $out_file = $ARGV[0];
+         shift;
+       }
+       elsif (/^-\w+/) {
+         die "$_: Illegal command line switch!\n";
+       }
+       else {
+         $file = $_;
+       }
+     }
+     die "No input file specified!\n" if $file eq "";
+
+  print "This is mathspic version $version_number\n";
+}
+   our ($source_file, $log_file);
+   if (! -e $file) {
+      die "$file: no such file!\n" if (! (-e "$file.m"));
+      $source_file = "$file.m";
+   }
+   else {
+      $source_file = $file;
+      $file = $1 if $file =~ /(\w[\w-\.]+)\.\w+/;
+   }
+   $out_file= "$file.mt" if $out_file eq "default";
+   $log_file= "$file.mlg";
+
+
+  open(IN,"$source_file")||die "Can't open source file: $source_file\n";
+  open(OUT,">$out_file")||die "Can't open output file: $out_file\n";
+  open(LOG,">$log_file")||die "Can't open log file: $log_file\n";
+  print_headers;
+  process_input(IN,"");
+
+print $alarm if $no_errors > 0;
+__END__
+