|
投稿者：しばっち
投稿日：2015年10月28日(水)22時05分25秒
|
|
|
|
!'Lambert's Method
INPUT R
LET D=2
LET X=1
DO
LET X=H(X)
PRINT X
LOOP UNTIL ABS(F(X))<1E-12
PRINT X;X^D
FUNCTION F(X)
LET F=X^D-R
END FUNCTION
FUNCTION H(X)
LET H=((D-1)*X^D+(D+1)*R)/((D+1)*X^D+(D-1)*R)*X
END FUNCTION
END
|
|
|