|
|
 |
| found 14 articles for 0.000 sec |
  |
Group: comp.lang.fortran · Group Profile · Search for Xy1 in comp.lang.fortran
Author: Ron Ford
Date: May 6, 2008 15:34
...BESSEL_Y0(x) real(kind(x)), parameter :: xy1 = BESSEL_Y1(x) real(kind(x)), ...', xjn write(*,*) 'xy0 = ', xy0 write(*,*) 'xy1 = ', xy1 write(*,*) 'xyn = ', xyn end ...= 0.13737752736232718572 xy0 = 0.18902194392082650675 xy1 = 0.41018841788751188289 xyn = 0.41018841788751188289 ...(x) real(kind(x)), parameter :: xy1 = BESSEL_Y1(x) real(kind(x)),...
|
| Show full article (7.83Kb) · Show article thread |
|
|
Group: comp.lang.fortran · Group Profile · Search for Xy1 in comp.lang.fortran
Author: James Van Buskirk
Date: May 6, 2008 12:50
... xy0 = BESSEL_Y0(x) real(kind(x)), parameter :: xy1 = BESSEL_Y1(x) real(kind(x)), parameter :: ...'xjn = ', xjn write(*,*) 'xy0 = ', xy0 write(*,*) 'xy1 = ', xy1 write(*,*) 'xyn = ', xyn end program test8 C...13737752736232718572 xjn = 0.13737752736232718572 xy0 = 0.18902194392082650675 xy1 = 0.41018841788751188289 xyn = 0.41018841788751188289 So mystery solved...
|
| Show full article (5.93Kb) · Show article thread |
Group: comp.lang.forth · Group Profile · Search for Xy1 in comp.lang.forth
Author: J Thomas
Date: Apr 20, 2007 20:12
Maybe there's a better way to use xy structures. \ y = y0 + (y1 - y0)*(x - x0)/(x1 - x0) \ : interp ( x xy0 xy1 -- y ) R R'@ 2@ ROT - NEGATE \ y0 x-x0 2@ 2@ d- \ y0 x-x0 y1-y0 x1-x0 */ + ; \ works only when y1 = y0 Otherwise y's borrow from x's. : interp ( x y0 x0 y1 x1 -- y ) 2OVER D- \ x y0 x0 y1-y0 x1-x0 R ROT - NEGATE \ y0 x-x0 R: y1...
|
| Show full article (0.47Kb) |
|
 |
|
1 ·
2 ·
next |
|
|