%% 
%%  Ein Beispiel der DANTE-Edition
%% 
%% 
%%  Copyright (C) 2010 Herbert Voss
%% 
%%  It may be distributed and/or modified under the conditions
%%  of the LaTeX Project Public License, either version 1.3
%%  of this license or (at your option) any later version.
%% 
%%  See http://www.latex-project.org/lppl.txt for details.
%% 
%% 
%% ==== 
% Show page(s) 1
%% 
\documentclass[]{exaarticle}
\pagestyle{empty}
\setlength\textwidth{352.81416pt}
\setlength\parindent{0pt}
\StartShownPreambleCommands
\usepackage{pst-plot,multido}
\usepackage[nomessages]{fp}
\newwrite\DataFileA
\newwrite\DataFileB
%
\newcommand{\DataFileNameA}{MyCurveA.dat}
\newcommand{\DataFileNameB}{MyCurveB.dat}
%
\def\Xpsz{6}
\def\Ypsz{6}
\newcommand{\PerCalc}[2]{%
  \def\ix{#1}\def\iy{#2}%
  \FPeval{\iY}{\Ypsz * \iy / (\iy + 1.0)}  % Calculates Y-coordinate
  \FPeval{\iX}{(\Ypsz - \iY)* \ix / \Xpsz} % Calculates X-coordinate
}% End perspective calculate

\StopShownPreambleCommands
\begin{document}
\begin{pspicture}(-6,0)(6,6)
  \multido{\I=-6+1}{13}{\psline(0,6)(\I,0)}
  % The horizontal perspective grid
  \multido{\IR=0+1}{13}{%
    \PerCalc{-6.0}{\IR}
    \FPneg{\iXX}{\iX}                         % Exploits symmetry
    \psline[linecolor=red](\iXX,\iY)(\iX,\iY)}% End of the horizontal grid
%
  \immediate\openout\DataFileA=\DataFileNameA
  \immediate\openout\DataFileB=\DataFileNameB
%
  \multido{\IRx=0+1}{90}{%
    \FPeval{\Ry}{ (\IRx * \IRx) / 30.0 }   % Defines the quadratic
    \PerCalc{\IRx}{\Ry}
    \immediate\write\DataFileA{(\iX,\iY)}
    \FPneg{\iXX}{\iX} % Exploits symmetry - so we'll have to loose this...
    \immediate\write\DataFileB{(\iXX,\iY)}}
%
  \immediate\closeout\DataFileA
  \immediate\closeout\DataFileB
%
  \readdata{\MyCurve}{\DataFileNameA}
  \dataplot[plotstyle=curve,%
	showpoints=true,%
	linecolor=green,linewidth=1.5pt]{\MyCurve}

  \readdata{\MyCurve}{\DataFileNameB}
  \dataplot[plotstyle=curve,%
	showpoints=true,%
	linecolor=blue,linewidth=1.5pt]{\MyCurve}
\end{pspicture}
\end{document}
