Discussion :

[JulienPari]

Bonjour,
Je travaille actuellement sur un projet informatique sur octave, cependant je fais face à un certain problème.
Je sais que cela viens de la taille d'une variable mais je ne vois absolument pas.
Voici mon programme:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
 
 
function [t,u]=euler(f,yint,tinit,tfin,Nh)
  u(1)=yint;
  h=(tfin-tinit)/Nh;
  t(1)=tinit;
  for n=1:Nh
    u(n+1)=u(n)+h*f(t(n),u(n));
    t(n+1)=t(n)+h;
  end
end
 
function fn=f(t,y)
  k1=1;
  k2=2;
  fn=zeros(size(y));
  fn(1)=-k1*y(1);
  fn(2)=k1*y(1)-k2*y(2);
  fn(3)=k2*y(3);
end
 
yint=0;
tinit=0;
tfin=5;
Nh=1000;
[T,U]=euler(@f,yint,tinit,tfin,Nh);
plot(T,U)
hold on;

Je vous remercie.

[Jerome Briot]

Montre nous le message d'erreur complet

[JulienPari]

Voici l'erreur:

error: y(2): out of bound 1
error: called from
f at line 22 column 8
euler at line 12 column 11
testdestes at line 30 column 5

Merci