Bonjour,

J'essayes d'appliquer la méthode des éléments finis mais j'obtiens un résultat incohérent.
Ce qui me gene un peu, c'est qu'on passe psi en paramètre sans préciser ses arguments.

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
 
funcprot(0)
 
function y=f(x)
    y=1;
endfunction
 
function I=integrale(a,b,r)
    H1=5/9;
    H2=8/9;
    H3=5/9;
    q=sqrt(3/5);
    e1=a+(1-q)*(b-a)/2;
    e2=(a+b)/2;
    e3=a+(1+q)*(b-a)/2;
    I=((b-a)/2)*(H1*r(e1)+H2*r(e2)+H3*r(e3));
endfunction
 
function y=phi(x)
y=max(0,1-abs(x));
endfunction
 
function y=psi(x,N,i)
h=1./(N+1);
y=phi((x-i*h)./h);
endfunction 
 
function b=CalculSMD(N, f)
    for i=1:N
        b(i)=integrale(0,1,psi)
        +integrale(i/(N+1),(i+1)/(N+1),psi);
    end
endfunction
Merci d'avance.