with Arg;
     with A83_014_Vector; use A83_014_Vector;
     with A83_015_Matrix; use A83_015_Matrix;
     with FloaAK; use FloaAK;
     with MathAK; use MathAK;

     procedure Compare is
     ----------------------------------------------------------------------
     ----------------------------------------------------------------------
       u     : Vector_Type ( 1..2 )    := ( 1.0, 0.99 ); -- neurons
       du_dt : Vector_Type ( u'range ); -- derivative of u
       V     : Vector_Type ( u'range ); -- output of u
       T : Matrix_Type ( u'range, u'range ) -- weight matrix
	 := ( (  0.0, -1.0 ),
	      ( -1.0,  0.0 ) );
       Tau  : constant float :=  1.0; -- RC time constant
       Beta : float := Value ( Arg.Data ( 2 ) ); -- gain
       dt   : float := Value ( Arg.Data ( 3 ) ); -- time delta
     begin
       for Index in 1..integer ( 10.0 / dt ) loop
	 V := TanH ( Beta * u );
	 du_dt := -u / Tau + T * V;
	 u := u + dt * du_dt;
	 Put_Line ( u );
       end loop;
     end Compare;