Update of the localization templates
[scilab.git] / scilab / modules / differential_equations / macros / inttrap.sci
1 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
2 // Copyright (C) INRIA
3 // 
4 // This file must be used under the terms of the CeCILL.
5 // This source file is licensed as described in the file COPYING, which
6 // you should have received as part of this distribution.  The terms
7 // are also available at
8 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
9
10 function v = inttrap(x,y)
11 //inttrap  Trapezoidal numerical integration.
12 //v = inttrap(x,y) computes the integral of y with respect to x using
13 //trapezoidal integration.  x and y must be vectors of the same dimension
14 //v = inttrap(y) computes the trapezoidal integral of y assuming unit
15 //spacing between the data points. 
16
17 [lhs,rhs]=argn(0)
18 if rhs<2 then
19   y=x;
20   v=sum(y(1:$-1) + y(2:$))/2;
21 else
22   if size(x,'*')<>size(y,'*') then
23     error(msprintf(gettext("%s: Wrong size for input arguments: Same size expected.\n"),"inttrap"));
24   end
25   x=x(:);y=y(:);
26   v=(x(2:$)-x(1:$-1))'*(y(1:$-1) + y(2:$))/2;
27 end
28 endfunction