6コメント2KB

MATLAB詳しい人来て！！

■ このスレッドは過去ログ倉庫に格納されています

0001デフォルトの名無しさん

2009/07/05(日) 03:57:19

dy/dx=f(x,y)=3x+3yでy(0)=1とします。区間[0,2]を20等分した点におけるy=y(x)の値をホイン法で求めるときホイン法のプログラムは

function [x,y]=Heun(a,b,n,y_0)
h=(b-a)/n;
x=a:n:b;
y=0*x;
y(1)=y_0;
for k=1:n
s_1=h*fval([x(k),y(k)]);
s_2=h*fval([x(k)+h,y(k)+s_1]);
y(k+1)=y(k)+(s_1+s_2)/2;
end
plot(x,y,'-g');

になるんですがMATLABではエラーが出ます
s_1=h*fval([x(k),y(k)]);この部分です

[x,y]=Heun(0,2,20,1)とうつと
??? Attempted to access x(2); index out of bounds because numel(x)=1.

エラー ==> Heun at 7
s_1=h*fval([x(k),y(k)]);と表示されます。

関数は
function z=fval(w)
x=w(1);
y=w(2);
z=3*x+3*y;としています。
何が違うのかわからないです。お願いします。

0002デフォルトの名無しさん

2009/07/05(日) 03:58:04

またアイちゃんか

0003デフォルトの名無しさん

2009/07/05(日) 04:16:54

だれかーーー

0004デフォルトの名無しさん

2009/07/05(日) 04:20:04

このスレッドは天才チンパンジー「アイちゃん」が
言語訓練のために立てたものです。

アイと研究員とのやり取りに利用するスレッドなので、
関係者以外は書きこまないで下さい。

　　　　　　　　　　　　　　　　　　京都大学霊長類研究所

0005デフォルトの名無しさん

2009/07/05(日) 07:58:32

>>1
>ホイン法で求めるときホイン法のプログラムは

こいつに聞いたらどうだ？　ホイン法のプロだぞ。

957 名前: デフォルトの名無しさん Mail: 投稿日: 2009/06/01(月) 01:44:36
2次のルンゲクッタ法のホイン法と中点法でプログラムを作っています。

958 名前: デフォルトの名無しさん Mail: 投稿日: 2009/06/01(月) 02:07:00
2次のルンゲクッタ法のホイン法と中点法でプログラムを作っています。

959 名前: デフォルトの名無しさん Mail: 投稿日: 2009/06/01(月) 02:09:28
2次のルンゲクッタ法のホイン法と中点法でプログラムを作っています。

960 名前: デフォルトの名無しさん Mail: sage 投稿日: 2009/06/02(火) 22:18:09
2次のルンゲクッタ法のホイン法と中点法でプログラムを作っています。

0006デフォルトの名無しさん

2009/07/05(日) 09:00:12

ボイン砲なら得意です

■ このスレッドは過去ログ倉庫に格納されています