2011-08-01から1ヶ月間の記事一覧

円と曲線の交点 #2

文関数で定義すると短く済む。 どんどん現代書法的に非推奨になる。いい感じ。 PROGRAM INTERS IMPLICIT NONE DOUBLE PRECISION Y,YD,Y1,Y2,X,XN,E,H,X0 PARAMETER(E=1.0D-12,H=1.0D-8,X0=0.9) Y1(X)=X*TAN(X) Y2(X)=SQRT(1-X**2) Y(X)=Y1(X)-Y2(X) YD(X)=(Y…