!  Testing 3-value arithmetic
 !  --------------------------

 program test
   use sea3_arithmetic                       !  Error analytic module
   type(sea_abs3) :: a, b, c, d, x           !  Declare our vars

   call sea_check_precise()                  !  Verify PRECISE reals exists

   sea_default_error = 0.1_precise                  !  Modify default

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

   call sea_write_entry(fout, "a = ", a, 5)
   call sea_write_entry(fout, "b = ", b, 5)
   call sea_write_entry(fout, "c = ", c, 5)
   call sea_write_entry(fout, "d = ", d, 5)
   call sea_write_nl(fout)                          !  Write a newline

   x = a * d - b * c

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

 end program test