スポンサーサイト

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

ほめてあげたい

きょう、帰ってきて、コタツをつけたら
なんと復活していますた。

よくやった! コタツ!! 
褒めて遣わすでゲソ!!

何にしても、この冬を気合だけで乗り切る事にはならなくてよかった。



あと、今日はベホマとカラオケに行ったんだけども。
その後、何だかめんどくさい式をプログラムして、符号がどうなるか調べた。

ほぼノーミスで1発コンパイル出来た自分を褒めてあげたい。


時間がある人は、この式をnが負の無限大から正の無限大まで
(数値的には無理なので、適当な所で切って)
計算してみて、この値の和が負になるか確かめてみてもいいんだからねっ!!

もしも正になったりしたら、忠告してくださいな。


==============================================================================

+(exp(-ener(UB,UFB,nF,n)/temp))*(
& +1.0/((ener(UB,UFB,nF,n+1)-ener(UB,UFB,nF,n))**2.0*temp*
& (ener(UB,UFB,nF,n+2)-ener(UB,UFB,nF,n)))
c
& -0.5/((ener(UB,UFB,nF,n+1)-ener(UB,UFB,nF,n))*temp**2.0*
& (ener(UB,UFB,nF,n+2)-ener(UB,UFB,nF,n+1)))
c
& +1.0/((ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n+1))**3.0*
& (ener(UB,UFB,nF,n+2)-ener(UB,UFB,nF,n)))
c
& +1.0/(ener(UB,UFB,nF,n+1)-ener(UB,UFB,nF,n))*
& (ener(UB,UFB,nF,n+2)-ener(UB,UFB,nF,n+1))*
& ( 1.0/(ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n+2))
& -1.0/(ener(UB,UFB,nF,n+2)-ener(UB,UFB,nF,n+1)) )* !!×
& ( 1.0/(ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n+1))
& -1.0/(ener(UB,UFB,nF,n+1)-ener(UB,UFB,nF,n+2)) ) ) !!
cc
cc
& -exp(-ener(UB,UFB,nF,n+1)/temp)/
& ((ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n+1))**3.0*
& (ener(UB,UFB,nF,n+2)-ener(UB,UFB,nF,n)))
c
c
& +1.0/(ener(UB,UFB,nF,n+1)-ener(UB,UFB,nF,n))*
& (ener(UB,UFB,nF,n+2)-ener(UB,UFB,nF,n+1))*
& ( 1.0/(ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n+2))
& -1.0/(ener(UB,UFB,nF,n+2)-ener(UB,UFB,nF,n+1)) )* !!×
c
& ( exp(-ener(UB,UFB,nF,n+2)/temp)/
& (ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n+2))
& -exp(-ener(UB,UFB,nF,n+1)/temp)/
& (ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n+1)))
c
c
& +(exp(-ener(UB,UFB,nF,n)/temp))*(
& +1.0/((ener(UB,UFB,nF,n-1)-ener(UB,UFB,nF,n))**2.0*temp*
& (ener(UB,UFB,nF,n-2)-ener(UB,UFB,nF,n)))
c
& -0.5/((ener(UB,UFB,nF,n-1)-ener(UB,UFB,nF,n))*temp**2.0*
& (ener(UB,UFB,nF,n-2)-ener(UB,UFB,nF,n-1)))
c
& +1.0/((ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n-1))**3.0*
& (ener(UB,UFB,nF,n-2)-ener(UB,UFB,nF,n)))
c
& +1.0/(ener(UB,UFB,nF,n-1)-ener(UB,UFB,nF,n))*
& (ener(UB,UFB,nF,n-2)-ener(UB,UFB,nF,n-1))*
& ( 1.0/(ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n-2))
& -1.0/(ener(UB,UFB,nF,n-2)-ener(UB,UFB,nF,n-1)) )* !!×
& ( 1.0/(ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n-1))
& -1.0/(ener(UB,UFB,nF,n-1)-ener(UB,UFB,nF,n-2)) ) ) !!
cc
cc
& -exp(-ener(UB,UFB,nF,n-1)/temp)/
& ((ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n-1))**3.0*
& (ener(UB,UFB,nF,n-2)-ener(UB,UFB,nF,n)))
c
c
& +1.0/(ener(UB,UFB,nF,n-1)-ener(UB,UFB,nF,n))*
& (ener(UB,UFB,nF,n-2)-ener(UB,UFB,nF,n-1))*
& ( 1.0/(ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n-2))
& -1.0/(ener(UB,UFB,nF,n-2)-ener(UB,UFB,nF,n-1)) )* !!×
c
& ( exp(-ener(UB,UFB,nF,n-2)/temp)/
& (ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n-2))
& -exp(-ener(UB,UFB,nF,n-1)/temp)/
& (ener(UB,UFB,nF,n)-ener(UB,UFB,nF,n-1)))


----------------------------
ただし、function: ener(x,y,z,n)=x*0.5*n**2.0 + 6.0*y*z*n


=============================================================================

う~む。それにしても長い……ww

一応、上手く行ったけど、これってそもそも手計算が合ってるんだろうか?w 
 
スポンサーサイト

comment

管理者にだけ表示を許可する

たまには本気返信

cが得意な人にソース見てもらえ
多分10分の1ぐらいになる

Re: たまには本気返信

> cが得意な人にソース見てもらえ
> 多分10分の1ぐらいになる

そうなの? てかcじゃなくfortranだけども。
愚直に式をなぞり過ぎたと今では反省している。
プロフィール

いろはうた

Author:いろはうた
いきおいだけで作ったブログです。
なんの気なしに見ちゃえばいいよっ!!
あと、勢いでコメントしてくれていいよっ!!!
リンクもドンドン貼っちゃえばいいよっ!!!!

twitter
twitterもやってます。 よければ、フォローしてやってください。
最新コメント
最新記事
マクロスF~名言集~


Spcial Thanks マクロスF
カレンダー
07 | 2017/08 | 09
- - 1 2 3 4 5
6 7 8 9 10 11 12
13 14 15 16 17 18 19
20 21 22 23 24 25 26
27 28 29 30 31 - -
FC2カウンター
カテゴリ
リンク
最新トラックバック
検索フォーム
RSSリンクの表示
ブロとも申請フォーム

この人とブロともになる

QRコード
QRコード
月別アーカイブ
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。