!  Testing 3-value arithmetic
 !  ---------------------------------------------
 !  A simple example of non-symmetric error range

 program test
   use sea3_arithmetic                       !  Error analytic module
   type(sea_abs3)        ::  x, y            !  Declare our vars

   call sea_check_precise()                  !  Verify PRECISE reals exists

   sea_default_error = 0.1_precise                  !  Modify default

   call sea_read_old(fin, x, "Enter x: ")           !  Input x
   call sea_write_nl(fout)                          !  Write a newline

   call sea_write_entry(fout, "x =    ", x, 5)
   call sea_write_nl(fout)                          !  Write a newline

   y = x ** 3

   call sea_write_entry(fout, "x**3 = ", y, 5)      !  Write result
   call sea_write_nl(fout)                          !  Write a newline

 end program test