28 lines
529 B
Mathematica
28 lines
529 B
Mathematica
|
function b=forcing(nu, dx1, dx2, N1, N2)
|
|||
|
|
%
|
|||
|
|
% Cette fonction construit le vecteur de for<EFBFBD>age de l'EDP
|
|||
|
|
%
|
|||
|
|
% Inputs
|
|||
|
|
% ------
|
|||
|
|
%
|
|||
|
|
% nu : nu=[nu1;nu2], coefficients de diffusivit<EFBFBD> dans les directions x1 et x2.
|
|||
|
|
%
|
|||
|
|
% dx1 : pas d'espace dans la direction x1.
|
|||
|
|
%
|
|||
|
|
% dx2 : pas d'espace dans la direction x2.
|
|||
|
|
%
|
|||
|
|
% N1 : nombre de points de grille dans la direction x1.
|
|||
|
|
%
|
|||
|
|
% N2 : nombre de points de grilles dans la direction x2.
|
|||
|
|
%
|
|||
|
|
% Outputs:
|
|||
|
|
% -------
|
|||
|
|
%
|
|||
|
|
% b : vecteur de for<EFBFBD>age (dimension N1N2)
|
|||
|
|
%
|
|||
|
|
%
|
|||
|
|
|
|||
|
|
b = -ones(N1*N2, 1);
|
|||
|
|
|
|||
|
|
end
|