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.
Merci d'avance.
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
Partager