bug 6744 fix
[scilab.git] / scilab / modules / cacsd / tests / nonreg_tests / bug_6744.dia.ref
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2010 - INRIA - Serge Steer
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 // <-- Non-regression test for bug 6744 -->
8 // <-- JVM NOT MANDATORY --> 
9 // <-- Bugzilla URL -->
10 // http://bugzilla.scilab.org/show_bug.cgi?id=6744
11 //
12 // <-- Short Description -->
13 //  p_margin() returns an erroneous result
14 z=poly(0,'z');
15 num=(0.0373327-0.0405775*z-0.1413155*z^2+0.1783403*z^3);
16 den=(-0.0373327+0.0405775*z+0.1413155*z^2-0.1783403*z^3+0.04867*z^4-0.04007*z^5-0.17482*z^6+0.2*z^7);
17 H=syslin(1.0,num/den);
18 [m,fr]=p_margin(H);
19 r=repfreq(H,fr);
20 //check if fr give a response with modulus=1
21 if abs(abs(r)-1)>1d-10 then bugmes();quit;end
22 if abs(m-58.03)>0.01 then bugmes();quit;end
23 if abs(m-pmodulo(atand(imag(r),real(r))-180,360))>1d-10 then bugmes();quit;end