Discussion :

[legyptien]

Bonjour,

je cherche à calculer l integrale d'une fonction qui est tres compliquee. J'ai eu des erreurs à cause de dimensions de matrice alors j'ai enormément simplifié ma fonction pour comprendre comment ca marche.

PROGRAMME PRINCIPAL:

Code :

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

1
2
clear all;
im = quad(@(x)zt(x),1,3)

PROGRAMME FONCTION:

Code :

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

1
2
3
4
%fonction
function f=zt(x)
size(x)
f=1./(x.^3);

-----------------------------------------------

Alors pourquoi j'obtiens:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
ans =
 
     1     7
 
 
ans =
 
     1     2
 
 
ans =
 
     1     2.....

La dimension de x reste 1 2 jusqu au bout. Pourquoi lors du premier appel à la fonction j'obtiens une dimensions de 1 7.

Merci

[Fiabulu]

Bonjour,

je pense qu'il s'agit d'une particularité de la méthode utilisé, si tu entre dans le code de la fonction quad, tu peux voir le commentaire suivant:

Code :

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

1
2
3
4
5
% Initialize with three unequal subintervals.
h = 0.13579*(b-a);
x = [a a+h a+2*h (a+b)/2 b-2*h b-h b];
y = feval(f, x, varargin{:});
fcnt = 7;

Après, pour comprendre pourquoi, il faut se plonger dans le détail de la méthode de Simpson