ttp://ime.sakura.ne.jp/cgibin/source/subupload.cgi/ime0005.bmp
c=0.0000001
R=1000
E=5
初期値x(0)=0
で上の画像の方形波3周期分をあたえたときにでる
Vrを求めるのが目的で
x=∫idt
x'=(-1/c*r)*x+(1/r)*e(t)
Vr=(-1/c)*x+e(t)
(x'はxを微分したもの)
この式をオイラー法で解くんですが
周期Tは0.00005と0.001の2つから選べるようにするんです
結果は時刻tとe(t)とVrが出るようにするんですが
どんなプログラムを組めばよろしいでしょうか?
どんなプログラムを書けばよろしいですか