Discussion :

[bsangoku]

Bonjour,

[bsangoku]

Je crois que le problème se situe dans les:

Code :

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

  .^

[phryte]

Bonjour.
Oui, il faut mettre ^ correctement.
Qu'est ce que i ?
i est bien le complexe (0,1) ?

[bsangoku]

Comment ça correctement?

Non: en faisant un copier coller j'ai oublier de mettre cette instruction

Code :

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

[x,i]=sort(x);

[phryte]

exemple :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
%function I = gauss(f,n) 
clear
f=inline('x');
n=10;
beta = 0.5./sqrt(1-(2*(1:n)).^(-2)); 
T = diag(beta,1) + diag(beta,-1); 
[V,D] = eig(T); 
x = diag(D); [x] 
[x,i]=sort(x);
w = 2*V(1,i).^2;
I = w*feval(f,x)

Il faut éviter d'utiliser i (réservé aux complexes).

[bsangoku]

Oui mais pourquoi on fait un inline?
Parce que l'utilisateur a le droit de choisir n'importe quel autre fonction du genre '2*x+1'...

[phryte]

Oui mais pourquoi on fait un inline?

J'ai fait le calcul directement sans faire de function.
inline définit une fonction.
par exemple :
f=inline('2*x+1')
Si tu ne le mets pas, tu définis une chaîne de caractères.

[bsangoku]

Pourquoi inline (x*x) ne marche pas?

[phryte]

Code :

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

inline('x.*x')