新しく発言する  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 
Re: 有難うございます。  返事を書く  ノートメニュー
山中和義 <drdlxujciw> 2008/03/22 13:40:09
!Richardsonの補外法(extrapolation)

DEF f(x)=x*EXP(x) !関数の定義

FUNCTION N(j,h) !O(h^(2*j))の近似
IF j=1 THEN
LET N=(f(x0+h)-f(x0-h))/(2*h)
ELSE
LET N=N(j-1,h/2) + (N(j-1,h/2)-N(j-1,h))/(4^(j-1)-1) !j=2,3,…
END IF
END FUNCTION


LET x0=2
LET h=0.2

PRINT N(1,h) !N1(h)近似

PRINT N(2,h) !N2(h)近似

PRINT N(3,h) !N3(h)近似

PRINT N(4,h) !N4(h)近似

PRINT N(5,h) !N5(h)近似


PRINT
PRINT x0*EXP(x0)+EXP(x0) !検算


END
  │   └こんなに短いコードで 島村1243 2008/03/22 17:14:23 
  可変精度数値微分 しばっち 2008/03/20 21:07:22  (修正1回)
  │├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
新規発言を反映させるにはブラウザの更新ボタンを押してください。