Discussion :

[Mataleb]

Bonjour tout le monde; Merci de m'aider à exécuter cet exemple (itération de Picard) sur Octave:

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
29
30
31
clc;clear all;close all;
pkg load symbolic
syms x y
% y'=2*x*y
fg= 2*x*y; % Eqt we interest to solve
fex=exp(x^2); % Exact solution
f=fg;
a=0;
b=1;
n=2;
h=(b-a)/n;
 
y0=1;  % Initial cdt
f= y0+subs(int(f,x,0,x),y,1);
 
for j=a:h:b
 % i=i+1;
  fn=subs(fg,y,f)
f=y0+int(fn,x);
end
i=0;zp(1)=0; zex(1)=0;
 
for j=a:h:b
    i=i+1;
    zp(i+1)=zp(i)+subs(f,x,j);
    zex(i+1)=zex(i) +subs(fex,x,j);
xx(i)=j;
end
plot(xx,zp(1:n+1),xx,zex(1:n+1))
legend('Picard method', 'Exact solution')
grid on

erreur:

error: operator =: no conversion for assignment of 'class' to indexed 'scalar'

[Jerome Briot]

Je ne suis pas un spécialiste du calcul symbolique mais d'après la documentation (lire au niveau du premier "Warning") :

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
29
30
31
32
33
34
clc;clear all;close all;
pkg load symbolic
syms x y
% y'=2*x*y
fg= 2*x*y; % Eqt we interest to solve
fex=exp(x^2); % Exact solution
f=fg;
a=0;
b=1;
n=2;
h=(b-a)/n;
 
y0=1;  % Initial cdt
f= y0+subs(int(f,x,0,x),y,1);
 
for j=a:h:b
 % i=i+1;
  fn=subs(fg,y,f)
f=y0+int(fn,x);
end
i=0;zp(1)=0; zex(1)=0;
 
fh = function_handle(f)
fexh = function_handle(fex)
 
for j=a:h:b
    i=i+1;
    zp(i+1)=zp(i)+fh(j);
    zex(i+1)=zex(i) + fexh(j);
xx(i)=j;
end
plot(xx,zp(1:n+1),xx,zex(1:n+1))
legend('Picard method', 'Exact solution')
grid on