[graphics] fix tests
[scilab.git] / scilab / modules / graphics / tests / nonreg_tests / bug_954.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2008-2008 - INRIA - Vincent COUVERT
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 //
8 // <-- ENGLISH IMPOSED -->
9 // <-- TEST WITH GRAPHIC -->
10 // <-- NO CHECK REF -->
11 //
12 // <-- Non-regression test for bug 954 -->
13 //
14 // <-- Bugzilla URL -->
15 // http://bugzilla.scilab.org/show_bug.cgi?id=954
16 //
17 // <-- Short Description -->
18 // Scilab is dead after running this script. It is the graphical
19 // commands at the end. Some test must be missing for
20 // huge number of points.
21
22 // Test #1
23 function Xd=f(t,X)
24     p=X(1);pi=X(2);xh=X(3);z=X(4);
25     pd=-(p^2+2*p+1/lam);
26     pid=1+2*(1/(lam*p)+1)*pi;
27     xhd=-(1+p)*xh;
28     zd=-1/(lam*p)*(z-xh)-(1+(1/pi))*z;
29     sigd=z^2;
30     Xd=[pd;pid;xhd;zd;sigd];
31 endfunction
32
33 lam=.8;
34
35 p0=1;
36 pi0=inv(p0);
37 xh0=2;
38 z0=xh0;
39 sig0=0;
40 TT=0:.0000001:0.43;
41 XX=ode([p0;pi0;xh0;z0;sig0],0,TT,f);
42
43 TT=TT(1:size(XX,2));
44
45 K=XX($,:)+XX(2,:).*(XX(3,:)-XX(4,:)).^2;
46 scf(0);
47 plot2d(TT',K');
48 delete(gcf());
49 scf(1);
50 plot2d(TT',XX(2,:)');
51 delete(gcf());
52 scf(2);
53 plot2d(TT',XX(1,:)');
54 delete(gcf());
55
56 // Test #2
57 clear
58 scf();
59 N=5d6;
60 TT=linspace(0,1,N);
61 plot2d(TT,sin(TT));
62 delete(gcf());