```/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        A I D E   A   L A   C O M P R E H E N S I O N   D U   P R O B L E M E   'v \$xiii/quad_image\$FON 20120618103427'  :         */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '\$xtc/GenDeform.01\$c' :                                                                                          */
/*                                                                                                                                   */
/*                    Jean-Francois Colonna (LACTAMME, 20120618120052).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

#include  "INCLUDES.01.I"

extern    double    cos();
extern    double    sin();

extern    double    drand48();
#define   RANDOM(amplitude)                                                                                                             \
(0.4*((1.0*(drand48()-0.5))*((amplitude)/((double)(N-N0)))))

#define   N0        1
#define   N         30

#define   RHO0      459.097442
#define   RHON      416.097159

#define   THETA0    3.778463
#define   THETAN    3.857411

#define   PHASE     (-PI/2)

main()
{
int       n;
double    XG1=390.545621,YG1=288.045621;

for       (n=N0 ; n<=N ; n++)
{
double    rho=(((RHON*(n-N0))+(RHO0*(N-n)))/(N-N0));
double    theta=(((THETAN*(n-N0))+(THETA0*(N-n)))/(N-N0))+PHASE;
double    x,y;

if        ((n>N0) && (n<N))
{
rho = rho+RANDOM(ABSO(RHON-RHO0));
theta = theta+RANDOM(ABSO(THETAN-THETA0));
}
else
{
}

x = XG1+(rho*cos(theta));
y = YG1+(rho*sin(theta));

printf("{%f,%f} = {%f,%f}\n",rho,theta,x,y);
/* On notera que :                                                                           */
/*                                                                                           */
/*                  PHASE=0        ==>  {x} monotone et {y} chaotique,                       */
/*                  PHASE=-PI/2    ==>  {x} chaotique et {y} monotone,                       */
/*                                                                                           */
/* (avec les bornes de 'rho' et 'theta' ci-dessus). Ainsi, la "chaoticite" change d'axe      */
/* lors d'un dephasage de pi/2).                                                             */
}
}
```