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
34
35
36
37
38
39
40
41
42
43
44
45
46
|
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) INRIA
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution. The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt
function []=errbar(x,y,em,ep)
// Rajoute des barres d'erreur sur un graphique 2D
// x et y decrivent les courbes (voir plot2d)
// em et ep sont deux matrices la barre d'erreur au point
// <x(i,j),y(i,j)> va de <x(i,j),y(i,j)-em(i,j)> a <x(i,j),y(i,j)+em(i,j)>
// x,y,em et ep sont donc des matrices (p,q), q courbes contenant chacunes
// p points.
// Exemple : taper errbar()
// x=0:0.1:2*%pi;
// y=<sin(x);cos(x)>';x=<x;x>';plot2d(x,y);
// errbar(x,y,0.05*ones(x),0.03*ones(x));
[lhs,rhs]=argn(0)
if rhs == 0 then // demo
x = 0:0.1:2*%pi;
y = [sin(x);cos(x)]';
x = [x;x]';
plot2d(x,y);
errbar(x,y,0.05*ones(x),0.03*ones(x));
return;
else
if rhs<4 then
error(msprintf(gettext("%s: Wrong number of input argument(s): At least %d expected.\n"), "errbar", 4));
end
end
show_window();
[n1,n2] = size(x);
y1 = matrix(y-em,1,n1*n2);
x1 = matrix(x,1,n1*n2);
y2 = matrix(y+ep,1,n1*n2);
xsegs([x1;x1],[y1;y2]);
my_axe = gca();
my_axe.clip_state = "clipgrf";
endfunction
|
