calcul numérique à travers le dblquad

**sergy01** · 23/05/2012, 12h21

Bonjour, j'ai un programme matlab dont j'ai fait le calcul formel d'une équation a deux variable t et u ,j'aimerais avoir le calcul par le dblquad c'est à dire le calcul numérique afin de comparer les deux résultats,

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
clear all
clc
 
%function F=binorm(t,u)
uu=sym(20);
%su=1.25;
su=sym(5/4);
ut=sym(8);
st=sym(2);
%a=0.9;
a=sym(9/10);
ro=0;
n=1;
syms t u
 
A1=1/(((a^(n-1))^2)*2*pi*st*su*sqrt(1-ro^2));
A2=((t-ut*a^(n-1))^2)/(((a^(n-1))^2)*st^2);
A3=2*ro*(((t-ut*a^(n-1))*(u-uu*a^(n-1)))/(st*su*((a^(n-1))^2)));
A4=((u-uu*a^(n-1))^2)/(((a^(n-1))^2)*su^2);
A5=-(1/(2*(1-(ro^2))))*(A2-A3+A4);
 
Fn=A1*exp(A5);
figure(1)
ezmesh(Fn,[1 15 10 25])
 
F=int(int(Fn,u,10,u),t,1,t);
figure(2)
ezmesh(F,[1 15 10 25])