2011-08-03から1日間の記事一覧
はさみうち法でやってみる。 PROGRAM INTERS IMPLICIT NONE DOUBLE PRECISION ERR,HIGH0,LOW0,Y1,Y2,X,Y,HIGH,LOW,MID PARAMETER(ERR=1.0D-10,HIGH0=1.0,LOW0=0.0) Y1(X)=X*TAN(X) Y2(X)=SQRT(1-X**2) HIGH=HIGH0 LOW=LOW0 DO MID=(HIGH+LOW)/2 Y=Y2(MID)-Y1…
その日の書き物のsnippets置き場
C|CSS||comp|Graphviz||phys||étoile|off-topic||一覧
C|C++|CSS|FORTRAN|Java|Lua|XML||comp|cairo|GMP/MPFR|gnuplot|Graphviz|GTK+|MTCTM||
math|phys||étoile|memo|off-topic||一覧
はさみうち法でやってみる。 PROGRAM INTERS IMPLICIT NONE DOUBLE PRECISION ERR,HIGH0,LOW0,Y1,Y2,X,Y,HIGH,LOW,MID PARAMETER(ERR=1.0D-10,HIGH0=1.0,LOW0=0.0) Y1(X)=X*TAN(X) Y2(X)=SQRT(1-X**2) HIGH=HIGH0 LOW=LOW0 DO MID=(HIGH+LOW)/2 Y=Y2(MID)-Y1…