新しく発言する EXIT インデックスへ

数値微分式

　　数値微分式 しばっち 2008/03/19 21:46:52 
　　├どのように使うのですか 島村1243 2008/03/20 10:01:11 
　　│└返事 しばっち 2008/03/20 13:22:28 
　　│　└有難うございます。 島村1243 2008/03/21 22:18:20 
　　│　　└!Richardsonの補外法(extrapolation) 山中和義 2008/03/22 13:40:09 
　　│　　　└こんなに短いコードで 島村1243 2008/03/22 17:14:23 
　　├可変精度数値微分 しばっち 2008/03/20 21:07:22  (修正1回)
Re: 数値微分式 返事を書く ノートメニュー
しばっち <dihjvcfsyu> 2008/03/20 21:07:22 ** この記事は1回修正されてます
可変精度数値微分


INPUT PROMPT "X=":X
LET H=1/1024
LET N=5 !'次数
LET M=(N+1)/2 !'中央差分
PRINT "微分値";DIFF(X,H,M,N)
PRINT "関数値";F(X)
END

EXTERNAL FUNCTION DIFF(XX,H,M,N)
DIM X(N),Y(N),A(N)
FOR I=1-M TO N-M
LET X(I+M)=XX+I*H
LET Y(I+M)=F(X(I+M))
NEXT I
FOR I=1 TO N
LET L=1
LET KK=1
FOR J=1 TO N
IF J<>I THEN
LET L=L*(X(I)-X(J))
LET A(KK)=X(J)
LET KK=KK+1
END IF
NEXT J
LET S1=0
FOR J=1 TO N-1
LET S=1
FOR K=1 TO N-1
IF K<>J THEN LET S=S*(XX-A(K))
NEXT K
LET S1=S1+S
NEXT J
LET SS=SS+S1*Y(I)/L
NEXT I
LET DIFF=SS
END FUNCTION

EXTERNAL FUNCTION F(X)
LET F=X^2
END FUNCTION
 
　　│├F=X^2を試したら 島村1243 2008/03/22 10:23:28  (修正1回)
　　││└ご指摘ありがとうございます。こちらのミス... しばっち 2008/03/22 14:59:11 
　　││　└完璧でした。 島村1243 2008/03/22 16:36:45 
　　│└可変精度高階数値微分 しばっち 2008/03/23 09:54:12 
　　└高階数値微分式 しばっち 2008/03/22 17:21:13

インデックスへ EXIT
新規発言を反映させるにはブラウザの更新ボタンを押してください。