1. Top » 
  2. スポンサー広告 » 
  3. Java ウサギとキツネ
  4. プログラミング全般 » 
  5. Java » 
  6. Java ウサギとキツネ

スポンサーサイト

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

Java ウサギとキツネ

計算機シミュレーションの時間。
ウサギとキツネについての連立一次微分方程式を適当に解く。


import java.io.*;

public class Report1{
public static void main(String[] args) throws IOException{

// rabbit and fox
double r0 = 3;
double f0 = 1;
double rt,rtt,ft,ftt;
double br = 1;
double df = 1;
double a = 0.2;
double b = 0.02;

// time
double t;
double t0 = 0.0;
double tt = 25.0;
double dt = 0.01;

// init
rt = ( 1.0 + ( br - a * f0 ) * dt ) * r0;
ft = ( 1.0 + ( b * r0 - df ) * dt ) * f0;

File file = new File("test.csv");
PrintWriter pw = new PrintWriter(new FileWriter(file));

//0
pw.println("0,"+rt+","+ft);

for(t = t0; t < tt; t += dt){
double tr = rt;
rtt = ( 1.0 + ( br - a * ft ) * dt ) * rt;
ftt = ( 1.0 + ( b * rt - df ) * dt ) * ft;

pw.println(t+","+rtt+","+ftt);

rt = rtt;
ft = ftt;
}
pw.close();
}
}

関連記事
スポンサーサイト

Comment

コメントフォーム
このエントリへコメントを書く
(任意)
(任意)
(任意)
(必須) HTMLタグは使用できません
(任意) ID生成と編集に使用します
(任意) 非公開コメントにする

Page Top

Trackback

Trackback URI
http://parabola0.blog86.fc2.com/tb.php/48-000b792e この記事にトラックバックする(FC2ブログユーザーのみ)

Page Top

訪問者

お引っこし。 http://parabola.sblog.jp/

プロフィール

parabola0

Author:parabola0
Twitter用ですが…。
プロフィール

最新記事
最新コメント
最新トラックバック
カテゴリ
月別アーカイブ
検索フォーム
リンク

このブログをリンクに追加する

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