5ちゃんねる ★スマホ版★ ■掲示板に戻る■ 全部 1- 最新50  

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

MATLAB詳しい人来て!!

1 :デフォルトの名無しさん: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;としています。
何が違うのかわからないです。お願いします。

2 :デフォルトの名無しさん:2009/07/05(日) 03:58:04
またアイちゃんか

3 :デフォルトの名無しさん:2009/07/05(日) 04:16:54
だれかーーー

4 :デフォルトの名無しさん:2009/07/05(日) 04:20:04
このスレッドは天才チンパンジー「アイちゃん」が
言語訓練のために立てたものです。

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

                  京都大学霊長類研究所

5 :デフォルトの名無しさん: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次のルンゲクッタ法のホイン法と中点法でプログラムを作っています。


6 :デフォルトの名無しさん:2009/07/05(日) 09:00:12
ボイン砲なら得意です

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

★スマホ版★ 掲示板に戻る 全部 前100 次100 最新50

read.cgi ver 05.04.00 2017/10/04 Walang Kapalit ★
FOX ★ DSO(Dynamic Shared Object)