立方根を求める関数と、積分の関数を組み合わせて、x1/3を積分する関数を作りたいのですが、積分範囲は 0から1 0から2 0から3 0から4
4から5 4から6 4から7 4から8 のものを求めたくて、あと一歩のところをどうかいていいか迷っています。
なるべくきれいな形にしたいです。範囲指定をどう表し、また、一番スマートでしょうか。どなたかご指導お願いいたします。。。
(><;)ちなみに下は途中まで出来たプログラムです!!
sub cbrt {
local($a)=@_;
local($x0,$x1);
$x0=1;
$x1=(2*$x0**3+$a)/(3*$x0**2);
while( abs($x1-$x0)> 1.e-15) {
$x0=$x1;
$x1=(2*$x0**3+$a)/(3*$x0**2);
}
return $x0;
}
sub simp { # シンプソン法による積分
local($a, $b, $n)=@_;
local($s)=0;
local($dx)=($b-$a)/$n;
local($x);
for($x=$a+$dx; $x<$b; $x += $dx ) {
$s += (2*cbrt($x-$dx/2)+cbrt($x))/3;
}
$s += cbrt($a)/6 -cbrt($b)/6;
return $s*$dx;
}
print simp(0,1,10000/2), "\n";
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- Visual Basic(VBA) VBAプログラミング 2 2022/11/27 12:07
- Visual Basic(VBA) VBAプログラミング 2 2022/11/27 12:13
- C言語・C++・C# C言語の課題が出たのですが自力でやっても分かりませんでした。 要素数がnであるint型の配列v2の並 3 2022/11/19 17:41
- その他(プログラミング・Web制作) python flask から fastapiへの移行時のエラー対処 1 2023/02/05 12:26
- 物理学 ポテンシャルが有限で不連続の時、右側の波動関数をφ1(x)、左側をφ2(x)とする。境界条件の「波動 2 2023/06/04 13:53
- 数学 重積分の積分領域について D={(x,y)∈R^2 | 0≦y≦x≦∞} で表される領域で、∫[0→ 3 2023/05/05 23:33
- 数学 ハイネボレルの被覆定理、内田伏一著 「集合と位相」定理22.1 1 2022/07/07 10:49
- 数学 三角関数の範囲について ∫1/√(a²-x²)dxをx=a・sin(t)と置いて置換積分する時tの範 3 2022/05/05 04:13
- 数学 解析学の問題がわかりません 1 2023/01/12 22:59
- 数学 3次対称群S3はシロー部分群で因数分解できない この問題の証明が分かりません。できる範囲で教えていた 1 2022/12/13 13:10
関連するカテゴリからQ&Aを探す
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
CLispのloop内の挙動について
-
Pythonのプログラミングの質問...
-
Application.ScreenUpdating = ...
-
パイソンのクラスについて
-
Pythonでターミナルに文字を出...
-
どちのほうがすきですか?
-
パイソンのクラスのブログラム
-
パイソンのクラスについて
-
パイソンプログラミング
-
【メモリ不足で落ちる(python)】
-
パソコン → パャRン
-
vscode 文字化け
-
[python] 文字列を変数名として...
-
メモリをアドレスを直接指定し...
-
JSONで文字列が長い時
-
実行時エラー 3020の対策
-
<SELECT>タグの折り返し
-
c++の構造体について
-
int(input("○○"))の使い方
-
テキスト入力後、エンターを押...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
Perlで複数の値を返す良い方法...
-
HASH(0xほげほげ)
-
プロトタイプ宣言のワーニング...
-
Pythonのプログラミングの質問...
-
Net::SSH::Perlについて
-
perlのスクリプトどこが間違っ...
-
print <<"_HTML_"の中でサブル...
-
内部処理形式が日付型 (Date) ...
-
Perlで use strict して our変数
-
変数に値を代入するのに、まずu...
-
Perlで足し算をするには
-
連想配列とforeach
-
Perl:計算誤差について
-
icmp socket error
-
Perlで自分自身の関数名を知る...
-
PERLの「!=」と「=~」 or...
-
[perl] next と return 同じ意...
-
for文の中でのforeachの問題
-
PerlでCabochaを呼び出すには(...
-
Net::Telnetでタイムアウト時間...
おすすめ情報