def kershawtri_determinant_test(): #*****************************************************************************80 # ## KERSHAWTRI_DETERMINANT_TEST tests KERSHAWTRI_DETERMINANT. # # Licensing: # # This code is distributed under the GNU LGPL license. # # Modified: # # 17 February 2015 # # Author: # # John Burkardt # from kershawtri import kershawtri from r8vec_uniform_ab import r8vec_uniform_ab from r8mat_print import r8mat_print print '' print 'KERSHAWTRI_DETERMINANT_TEST' print ' KERSHAWTRI_DETERMINANT computes the KERSHAWTRI determinant.' n = 5 x_n = ((n + 1) // 2) r8_lo = -5.0 r8_hi = +5.0 seed = 123456789 x, seed = r8vec_uniform_ab(x_n, r8_lo, r8_hi, seed) a = kershawtri(n, x) m = n r8mat_print(m, n, a, ' KERSHAWTRI matrix:') value = kershawtri_determinant(n, x) print '' print ' Value = %g' % (value) print '' print 'KERSHAWTRI_DETERMINANT_TEST' print ' Normal end of execution.' return
def kershawtri_determinant_test ( ): #*****************************************************************************80 # ## KERSHAWTRI_DETERMINANT_TEST tests KERSHAWTRI_DETERMINANT. # # Licensing: # # This code is distributed under the GNU LGPL license. # # Modified: # # 17 February 2015 # # Author: # # John Burkardt # from kershawtri import kershawtri from r8vec_uniform_ab import r8vec_uniform_ab from r8mat_print import r8mat_print print '' print 'KERSHAWTRI_DETERMINANT_TEST' print ' KERSHAWTRI_DETERMINANT computes the KERSHAWTRI determinant.' n = 5 x_n = ( ( n + 1 ) // 2 ) r8_lo = -5.0 r8_hi = +5.0 seed = 123456789 x, seed = r8vec_uniform_ab ( x_n, r8_lo, r8_hi, seed ) a = kershawtri ( n, x ) m = n r8mat_print ( m, n, a, ' KERSHAWTRI matrix:' ) value = kershawtri_determinant ( n, x ) print '' print ' Value = %g' % ( value ) print '' print 'KERSHAWTRI_DETERMINANT_TEST' print ' Normal end of execution.' return
def kershawtri_test(): #*****************************************************************************80 # ## KERSHAWTRI_TEST tests KERSHAWTRI. # # Licensing: # # This code is distributed under the GNU LGPL license. # # Modified: # # 17 February 2015 # # Author: # # John Burkardt # from r8vec_uniform_ab import r8vec_uniform_ab from r8mat_print import r8mat_print print '' print 'KERSHAWTRI_TEST' print ' KERSHAWTRI computes the KERSHAWTRI matrix.' n = 5 x_n = ((n + 1) // 2) x_lo = -5.0 x_hi = +5.0 seed = 123456789 x, seed = r8vec_uniform_ab(x_n, x_lo, x_hi, seed) a = kershawtri(n, x) m = n r8mat_print(m, n, a, ' KERSHAWTRI matrix:') print '' print 'KERSHAWTRI_TEST' print ' Normal end of execution.' return
def kershawtri_test ( ): #*****************************************************************************80 # ## KERSHAWTRI_TEST tests KERSHAWTRI. # # Licensing: # # This code is distributed under the GNU LGPL license. # # Modified: # # 17 February 2015 # # Author: # # John Burkardt # from r8vec_uniform_ab import r8vec_uniform_ab from r8mat_print import r8mat_print print '' print 'KERSHAWTRI_TEST' print ' KERSHAWTRI computes the KERSHAWTRI matrix.' n = 5 x_n = ( ( n + 1 ) // 2 ) x_lo = -5.0 x_hi = +5.0 seed = 123456789 x, seed = r8vec_uniform_ab ( x_n, x_lo, x_hi, seed ) a = kershawtri ( n, x ) m = n r8mat_print ( m, n, a, ' KERSHAWTRI matrix:' ) print '' print 'KERSHAWTRI_TEST' print ' Normal end of execution.' return