- 固定ランキング！
18 ：学生さんは名前がない[]：2011/05/02(月) 00:23:39.72 ID:FCP0/jky0 - Ａ ドラドリ
Ｂ 読書、法３、講師
Ｃ べんじー、桃
|
- トリコってアニメ化されてるけどそんなに面白いか？
4 ：ドラドリ♀д♀[]：2011/05/02(月) 00:27:14.92 ID:FCP0/jky0 - 今のジャンプの中では3番目に面白い 一番はスケットダンス
|
- トリコってアニメ化されてるけどそんなに面白いか？
8 ：ドラドリ♀д♀[]：2011/05/02(月) 00:29:17.53 ID:FCP0/jky0 - スケットダンスは最初バカにしてたけど、ロマンの話みてハマった。
|
- トリコってアニメ化されてるけどそんなに面白いか？
13 ：ドラドリ♀д♀[]：2011/05/02(月) 00:34:21.44 ID:FCP0/jky0 - 何の能力もない小松が活躍してるおかげで
単純なバトルものとは違って面白くなってる
でもハンターハンター再開したら、たぶん読まなくなると思う
|
- 東風平←お前ら読めるか？
5 ：ドラドリ♀д♀[]：2011/05/02(月) 22:23:53.89 ID:FCP0/jky0 - >>1 数値計算プログラミングスレ立てて
|
- 数値計算プログラミングスレ
2 ：ドラドリ♀д♀[]：2011/05/02(月) 22:46:25.81 ID:FCP0/jky0 - ルンゲクッタで方程式
x^-2(d/dx)x^2(dy/dx)=-y
っていうのを解きたいんだが、うまくいかない
誰か助けてください
|
- 数値計算プログラミングスレ
4 ：ドラドリ♀д♀[]：2011/05/02(月) 22:48:44.00 ID:FCP0/jky0 - >>3
スレ立てサンクスです
答えはsin(x)/xってのは分かってる
でも一致しない
|
- 数値計算プログラミングスレ
5 ：ドラドリ♀д♀[]：2011/05/02(月) 22:50:53.57 ID:FCP0/jky0 - program main
implicit none
real(8) x,y1,y2,dx,k1,k2,k3,k4,l1,l2,l3,l4,h
real(8) f,g
integer i
open(100, file='data.dat', status='replace')
h = 0.001d0
dx = 0.01d0
! initial condition
x = h
y1 = 1.0d0 - h**2 / 6.0d0 + h**4 / 120.0d0
y2 = - h**3 / 3.0d0 + h**5 / 30.0d0
write(100,*)"initial" ,x , y1, y2
|
- 数値計算プログラミングスレ
6 ：ドラドリ♀д♀[]：2011/05/02(月) 22:51:24.63 ID:FCP0/jky0 - ! Runge-Kutta method
write(100,*)"x, y1, y2"
do i=1,100
k1 = dx * f(y2,x)
k2 = dx * f(y2 + k1 / 2.0d0, x + dx / 2.0d0)
k3 = dx * f(y2 + k2 / 2.0d0, x + dx / 2.0d0)
k4 = dx * f(y2 + k3 , x + dx)
if(i ==1 )then
write(100,*)"a=",y1
endif
y1 = y1 + (k1 + 2.0d0 * k2 + 2.0d0 * k3 ) / 6.0d0
if(i ==1 )then
write(100,*)"b=",y1
endif
l1 = dx * g(y1,x)
l2 = dx * g(y1 + l1 / 2.0d0, x + dx / 2.0d0)
l3 = dx * g(y1 + l2 / 2.0d0, x + dx / 2.0d0)
l4 = dx * g(y1 + l3 , x + dx)
y2 = y2 + (l1 + 2.0d0 * l2 + 2.0d0 * l3 + l4) / 6.0d0
x = x + dx
write(100,10),x,y1,y2,sin(x) / x
10 format(f6.2,2x,f6.2,2x,f6.2,3x,f6.2)
if (y1 <= 0) then
exit
endif
enddo
end
|
- 数値計算プログラミングスレ
7 ：ドラドリ♀д♀[]：2011/05/02(月) 22:52:02.13 ID:FCP0/jky0 - ! definition of function
real(8) function f(y2,x)
implicit none
real(8) y2,x
f = x**(-2) * y2
return
end function f
real(8) function g(y1,x)
implicit none
real(8) y1,x
g = - x**2 * y1
return
end function g
|
- 数値計算プログラミングスレ
8 ：ドラドリ♀д♀[]：2011/05/02(月) 22:53:07.99 ID:FCP0/jky0 - これがソースコードです
どうやら、1step目のk4の値がめちゃくちゃでかくなってしまってるのが
原因のような気がするのだが
何が間違っているのか分からないです
|
- トリコってアニメ化されてるけどそんなに面白いか？
59 ：ドラドリ♀д♀[]：2011/05/02(月) 22:57:05.25 ID:FCP0/jky0 - >>50 みつどもえ謎の長期休載
|
- 数値計算プログラミングスレ
11 ：ドラドリ♀д♀[]：2011/05/02(月) 22:57:21.08 ID:FCP0/jky0 - うおーーー！！何が原因なんだー！！
|
- 大学生にもなって今さら聞けないこと
50 ：学生さんは名前がない[]：2011/05/02(月) 23:00:36.31 ID:FCP0/jky0 - 髪の毛染めに行ったら
ワックスつけたままの状態で薬品塗られたんだけど
化学反応とかしないよね？
|
- 数値計算プログラミングスレ
12 ：ドラドリ♀д♀[]：2011/05/02(月) 23:01:00.52 ID:FCP0/jky0 - オラオラ 大学生の力みせてみろや！！
|
- 数値計算プログラミングスレ
15 ：ドラドリ♀д♀[]：2011/05/02(月) 23:04:19.42 ID:FCP0/jky0 - >>14 この場合どこで質問すればいいんだ？
|
- 数値計算プログラミングスレ
20 ：ドラドリ♀д♀[]：2011/05/02(月) 23:07:36.48 ID:FCP0/jky0 - ム板ってなんだよ
|
- 数値計算プログラミングスレ
22 ：ドラドリ♀д♀[]：2011/05/02(月) 23:09:09.59 ID:FCP0/jky0 - >>21 物理系はほとんどフォートランでしょ Ｃの派閥もあることはあるけど
|
- 数値計算プログラミングスレ
25 ：ドラドリ♀д♀[]：2011/05/02(月) 23:11:52.78 ID:FCP0/jky0 - >>23 なんか数学の質問するのは場違いな気がする
|
- 数値計算プログラミングスレ
29 ：ドラドリ♀д♀[]：2011/05/02(月) 23:19:22.04 ID:FCP0/jky0 - >>26
オイラー法だと大丈夫だった
たぶんk3かk4が間違ってるんだと思う
|
- 数値計算プログラミングスレ
31 ：ドラドリ♀д♀[]：2011/05/02(月) 23:19:53.65 ID:FCP0/jky0 - >>27 誰に向かってクチ聞いてんだコラ 殺すぞ
|
- トリコってアニメ化されてるけどそんなに面白いか？
63 ：ドラドリ♀д♀[]：2011/05/02(月) 23:23:35.66 ID:FCP0/jky0 - >>60 とりあえず半年だ・・・半年だけ待ってやろう
|
- スカトロAVに出たことあるけど何か質問ある？
4 ：ドラドリ♀д♀[]：2011/05/02(月) 23:24:13.64 ID:FCP0/jky0 - 本物のウンチつかってるんですか？
|
- 数値計算プログラミングスレ
35 ：ドラドリ♀д♀[]：2011/05/02(月) 23:28:22.32 ID:FCP0/jky0 - ごめんなさい一か所間違い
! Runge-Kutta method
write(100,*)"x, y1, y2"
do i=1,100
k1 = dx * f(y2,x)
k2 = dx * f(y2 + k1 / 2.0d0, x + dx / 2.0d0)
k3 = dx * f(y2 + k2 / 2.0d0, x + dx / 2.0d0)
k4 = dx * f(y2 + k3 , x + dx)
if(i ==1 )then
write(100,*)"a=",y1
endif
y1 = y1 + (k1 + 2.0d0 * k2 + 2.0d0 * k3 + k4) / 6.0d0 !←k4抜けてた
if(i ==1 )then
write(100,*)"b=",y1
endif
|
- 数値計算プログラミングスレ
36 ：ドラドリ♀д♀[]：2011/05/02(月) 23:37:44.64 ID:FCP0/jky0 - プログラム板で聞いてくるﾉｼ
|
- 数値計算プログラミングスレ
38 ：ドラドリ♀д♀[]：2011/05/02(月) 23:43:29.72 ID:FCP0/jky0 - プログラム板めっちゃ過疎ってるやん
やっぱりここで答えを待つわ
|
- 数値計算プログラミングスレ
39 ：ドラドリ♀д♀[]：2011/05/02(月) 23:45:25.90 ID:FCP0/jky0 - >>37
これエムデン方程式(Lane emden eq)っていうんだけど
ググったらみんなルンゲクッタで解いてたよ
|
- トリコってアニメ化されてるけどそんなに面白いか？
65 ：ドラドリ♀д♀[]：2011/05/02(月) 23:48:02.79 ID:FCP0/jky0 - チャンピオンって素人？みたいなヤツを平気で採用してくるよな
特にショートギャグ枠
バイトなうとかましののとかうれいの日々とか
よくいえば新人発掘に力を入れてるってことになるんだろうけど
個人的にはツギハギ生徒会を復活させてほしい！！
|
- トリコってアニメ化されてるけどそんなに面白いか？
66 ：ドラドリ♀д♀[]：2011/05/02(月) 23:48:38.32 ID:FCP0/jky0 - まあジャンプに話を戻すと
バクマンとスケットダンスが面白いです
|
- 宮崎羽衣ちゃんっていまなにしてんの？
4 ：ドラドリ♀д♀[]：2011/05/02(月) 23:58:57.53 ID:FCP0/jky0 - あやひーのウンチ
|
- 数値計算プログラミングスレ
40 ：ドラドリ♀д♀[]：2011/05/02(月) 23:59:33.61 ID:FCP0/jky0 - なんかいいことねーかなー
|