Discussion :

[follichon]

Bonsoir,

Je cherche à intégrer numériquement sous MATLAB une fonction. J'ai écrit ce code pour cela :

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
clear all;
close all;
 
%n=10;
%h=(b-a)/n;
r=-40:1:40;
r=[]
 
%Lambda=2.2;
%wo=4.25;
%zo=1;
I=sqrt(-1);
 
%r=[];
%r=-40:1:40;
S=[]
 
for k=1:length(r)
 
    %f=inline('(x.^y)*') 
 
 
S=dblquad(flux,0,1,0,1);
 
 
end

avec le script de définition de la fonction :

Code :

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

1
2
3
4
5
6
function f = flux(x,y,r)
 
x=[]
y=[]
 
(x.^y)*r;

Ce que je veux faire c'est calculer l'intégrale double de cette fonction pour différentes valeurs de r pour ensuite représenter ces valeurs en fonction de mon paramètre r.

Mon premier problème ici que Matlab m'indique que "x" n'est pas défini et j'ai un souci avec ma fonction flux où Matlab n'arrive pas à intégrer avec dblquad que j'utilise pour ma fonction définie avec mon paramètre r.

Si quelqu'un voit comment modifier le programme pour ce calcul...

Merci

[Invité]

Bonsoir,

Voir Que représente le symbole @ ?

Remarque : I=sqrt(-1); peut s'écrire simplement i ou j