修正歐拉法計(jì)算
返回
<![CDATA[
###############################
# 修正歐拉法計(jì)算
# (dx/dt)2=-x
# 分解為 dv/dt = -x
# dx/dt = v
###############################
use strict;
my(%t,$t1,$x,$v,$ref);
$x = 1.0;
$v = 0.0;
$t{dt} = 0.01;
$t{tmax} = 100;
open(OUT,">euler1.csv");
print OUT 'NO,x,v,',"\n";
for ($t1=0; $t1<$t{tmax}; $t1+=$t{dt}) {
$t{k00} = $v;
$t{k01} = -$x;
$t{k10} = $v+$t{dt}*$t{k01};
$t{k11} = -($x+$t{dt}*$t{k00});
$x = $x + $t{dt}*($t{k00}+$t{k10})/2.0;
$v = $v + $t{dt}*($t{k01}+$t{k11})/2.0;
print OUT $t1,',';
print OUT $x,',';
print OUT $v,',';
print OUT "\n";
}
close(OUT);
__END__;
]]>
返回
张家川|
栖霞市|
宝应县|
集安市|
东乡族自治县|
淅川县|
彭泽县|
长岭县|
广灵县|
思南县|
温州市|
永仁县|
大方县|
江都市|
招远市|
江川县|
珠海市|
遂溪县|
三门县|
广丰县|
前郭尔|
自治县|
平和县|
青龙|
黄骅市|
西吉县|
柘城县|
新民市|
万宁市|
樟树市|
洞头县|
富顺县|
江永县|
新余市|
琼海市|
麦盖提县|
丰城市|
马龙县|
泾川县|
克山县|
光山县|