27a93e1554e3f629ca67ff7147f2774e03ad7bb0
[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 // <-- Non-regression test for bug 954 -->
11 //
12 // <-- Bugzilla URL -->
13 // http://bugzilla.scilab.org/show_bug.cgi?id=954
14 //
15 // <-- Short Description -->
16 // Scilab is dead after running this script. It is the graphical
17 // commands at the end. Some test must be missing for
18 // huge number of points.
19
20 // Test #1
21 function Xd=f(t,X)
22     p=X(1);pi=X(2);xh=X(3);z=X(4);
23     pd=-(p^2+2*p+1/lam);
24     pid=1+2*(1/(lam*p)+1)*pi;
25     xhd=-(1+p)*xh;
26     zd=-1/(lam*p)*(z-xh)-(1+(1/pi))*z;
27     sigd=z^2;
28     Xd=[pd;pid;xhd;zd;sigd];
29 endfunction
30
31 lam=.8;
32
33 p0=1;
34 pi0=inv(p0);
35 xh0=2;
36 z0=xh0;
37 sig0=0;
38 TT=0:.0000001:1;
39 XX=ode([p0;pi0;xh0;z0;sig0],0,TT,f);
40
41 TT=TT(1:size(XX,2));
42
43 K=XX($,:)+XX(2,:).*(XX(3,:)-XX(4,:)).^2;
44 scf(0);
45 plot2d(TT',K');
46 delete(gcf());
47 scf(1);
48 plot2d(TT',XX(2,:)');
49 delete(gcf());
50 scf(2);
51 plot2d(TT',XX(1,:)');
52 delete(gcf());
53
54 // Test #2
55 clear
56 scf();
57 N=5d6;
58 TT=linspace(0,1,N);
59 plot2d(TT,sin(TT));
60 delete(gcf());