From bug-octave-request at che dot utexas dot edu  Wed Oct 20 09:38:16 1993
Subject: problem still there...
From: "A. Scottedward Hodel" <A dot S dot Hodel at eng dot auburn dot edu>
To: bug-octave at che dot utexas dot edu
Cc: Bruce dot Tenison at eng dot auburn dot edu
Date: Wed, 20 Oct 93 09:38:04 CDT

----------
X-Sun-Data-Type: text
X-Sun-Data-Description: text
X-Sun-Data-Name: text
X-Sun-Content-Lines: 45

John,

The m-file below still fails.  I include the patch that you originally sent
me for your reference...

Test case: SparcStation:
---------------------------------------------------
scotte at strangelove dot 19% ~/bin/octave
Octave, version 0.76.  Copyright (C) 1992, 1993, John W. Eaton.
This is free software with ABSOLUTELY NO WARRANTY.
For details, type `warranty'.

octave:1> freq
Frequency response example
gpoles =

   -1.00000 - 1.00000i   -1.00000 + 1.00000i   -0.75000 + 0.00000i

npoles = 3
gzeros =

   0.50000 - 0.50000i   0.50000 + 0.50000i

nzeros = 2
Frequency response range:
Enter minimum frequency of plot (rad/s):1
fr_min = 1
Enter maximum frequency of plot (rad/s):100
fr_max = 100
Enter number of points in frequency plot:100
warning: implicit conversion of complex value to real value
error: invalid  index = 0
error: evaluating assignment expression near line 18, column 8
fr_npts = 100
octave:2> quit

The line in question is:
13: for i=1:length(om)
14: 15:  gs(i) = 1;
16:  for j=1:nzeros
17:    gs(i) = gs(i)*(om(i) - gzeros(j));
18: endfor


Scotte
----------
X-Sun-Data-Type: mail-file
X-Sun-Data-Name: done.lex.patch
X-Sun-Encoding-Info: uuencode
X-Sun-Content-Lines: 47

begin 600 done.lex.patch
M1G)O;2! 1F5D97 at N37-F8RY dot 87-A+D=O=CIJ=V5 8VAE+G5T97AA<RYE9'4@
M5V5D($]C=" Q,R Q-#HQ,3HR." Q.3DS"E)E='5R;BU0871H.B \0$9E9&5X
M+DUS9F,N3F%S82Y';W8Z:G=E0&-H92YU=&5X87,N961U/ at I296-E:79E9#H@
M9G)O;2!&961E>"Y-<V9C+DYA<V$N1V]V(&)Y(&5N9RYA=6)U<FXN961U(" at T
M+C$O4TU)+30N,"!E9&ES;VX at ,2XR,2D*"6ED($%!,C S.#D[(%=E9"P@,3,@
M3V-T(#DS(#$T.C$Q.C(W($-$5 I296-E:79E9#H at 9G)O;2!#3TY44D],4RY-
M4T9#+DY!4T$N1T]6(&)Y($9E9&5X+DUS9F,N3F%S82Y';W8 at =VET:"!)3E1%
M4DY%5" ["B  at (" @(" @("!7960L(#$S($]C=" Y,R Q-#HQ,CHP,B!#1%0*
M4F5C96EV960Z(&9R;VT at <V-H;V-H+F-H92YU=&5X87,N961U(&)Y($-/3E12
M3TQ3+DU31D,N3D%302Y'3U8 at *#0N,2]334DM-"XQ*0H):60@04$P-S@R,CL@
M5V5D+" Q,R!/8W0 at  dot 3,@,30Z,C<Z,# @0T14"E)E8V5I=F5D.B!F<F]M(&QO
M8V%L:&]S=" H:G=E0&QO8V%L:&]S="D at 8GD@<V-H;V-H+F-H92YU=&5X87,N
M961U(" at X+C8O dot "XV*2!W:71H(%--5% @:60@3T%!,C@U-#<@9F]R(#QH;V1E
M;&%S0$-/3E123TQ3+FUS9F,N;F%S82YG;W8^.R!7960L(#$S($]C=" Q.3DS
M(#$T.C$Q.C$S("TP-3 P"DUE<W-A9V4M260Z(#PQ.3DS,3 Q,S$Y,3$N3T%!
M,C at U-#= <V-H;V-H+F-H92YU=&5X87,N961U/@I4;SH@:&]D96QA<T!#3TY4
M4D],4RYM<V9C+FYA<V$N9V]V("A!+B!38V]T="!(;V1E;"D*4W5B:F5C=#H at 
M4F4Z(&TM9FEL92!Q=65S=&EO:6X at "DEN+5)E<&QY+51O dot B!9;W5R(&UE<W-A
M9V4 at ;V8@5V5D+" Q,R!/8W0@.3,@,#DZ,# Z,#4@0T14"D1A=&4Z(%=E9"P@
M,3, at 3V-T(#DS(#$T dot C$Q dot C$R($5$5 I&<F]M.B!J=V5 8VAE+G5T97AA<RYE
M9'4*0V]N=&5N="U,96YG=& at Z(#$P,C,*6"U,:6YE<SH@-#$*4W1A='5S dot B!2
M3PH* dot B!4:&4 at ;2UF:6QE(&)E;&]W(')U;G,@9FEN92!U;F1E<B!O8W1A=F4M
M,"XW-"P at 8G5T(&=I=F5S(&5R<F]R<R!U;F1E<@HZ(&]C=&%V92TP+C<V+@H*
M5&AI<R!I<R!A;B!/8W1A=F4 at 8G5G+B @5&AE('!R;V)L96T@:7,@=&AE(&QO
M;W *"CH at 9F]R(&D],3IL96YG=&@H;VTI"CH@"CH@(&=S*&DI(#T@,3L*(" @
M+BXN"CH at 96YD9F]R"@I/8W1A=F4@:7,@<')E;6%T=7)E;'D@8FEN9&EN9R!T
M:&4 at ;&]O<"!C;W5N=&5R(&!I)R!T;R!T:&4@9VQO8F%L"G-Y;6)O;"!@:2<N
M" at I(97)E)W,@82!P;W-S:6)L92!F:7@N("!)(&AA=F4@;F]T('1E<W1E9"!I
M="!A(&QO="P at 8G5T(&ET('-E96US('1O"G=O<FLL(&%N9"!)(&1O;B=T('1H
M:6YK(&ET('=I;&P at 8V%U<V4@86YY('!R;V)L96US("AH92!S87ES+"!F:6YG
M97)S"F-R;W-S960I+ at H*:G=E"@HJ*BH@;&5X+FQ^"3$Y dot 3,O,3 O,#4@,#0Z
M,S0Z,C<*+2TM(&QE>"YL"3$Y.3,O,3 O,3, at ,3@Z-#<Z-#0**BHJ*BHJ*BHJ
M*BHJ*BHJ"BHJ*B Q,#(P+#$P,C8 at *BHJ*@H@(" @("!R971U<FX@;'-R.PH@
M( H at (" @<WEM8F]L7W)E8V]R9" J9W-R(#T@9VQO8F%L7W-Y;5]T86(M/FQO
M;VMU<" H;F%M92P at ,"P@,"D["B$@("!I9B H8W5R<E]S>6U?=&%B(#T]('1O
M<%]L979E;%]S>6U?=&%B("8F(&=S<B A/2 H<WEM8F]L7W)E8V]R9" J*2!.
M54Q,*0H at (" @("!R971U<FX@9W-R.PH@( H@(" @<F5T=7)N(&-U<G)?<WEM
M7W1A8BT^;&]O:W5P("AN86UE+" Q+" P*3L*+2TM(#$P,C L,3 R-R M+2TM
M"B  at (" @(')E='5R;B!L<W(["B @"B @("!S>6UB;VQ?<F5C;W)D("IG<W(@
M/2!G;&]B86Q?<WEM7W1A8BT^;&]O:W5P("AN86UE+" P+" P*3L*(2  at (&EF
M("AC=7)R7W-Y;5]T86( at /3T@=&]P7VQE=F5L7W-Y;5]T86(@)B8@9W-R("$]
M("AS>6UB;VQ?<F5C;W)D("HI($Y53$P*(2  at (" @(" F)B A("AL;V]P:6YG
M('Q\(&EF9FEN9RDI"B  at (" @(')E='5R;B!G<W(["B @"B @("!R971U<FX@
E8W5R<E]S>6U?=&%B+3YL;V]K=7  at *&YA;64L(#$L(# I.PH*"B!R
 
end
----------
X-Sun-Data-Type: default
X-Sun-Data-Name: freq.m
X-Sun-Content-Lines: 37

disp('Frequency response example')
gpoles = [ -1-i, -1+i, -0.75]
npoles = length(gpoles)
gzeros = [ (1-i)/2, (1+i)/2]
nzeros = length(gzeros)

disp('Frequency response range:')
fr_min = input('Enter minimum frequency of plot (rad/s):')
fr_max = input('Enter maximum frequency of plot (rad/s):')
fr_npts = input('Enter number of points in frequency plot:')
om = i*logspace(log10(fr_min),log10(fr_max),fr_npts);

for i=1:length(om)

 gs(i) = 1;
 for j=1:nzeros
   gs(i) = gs(i)*(om(i) - gzeros(j));
 endfor

 for j=1:npoles
   gs(i) = gs(i)/(om(i) - gpoles(j));
 endfor

endfor

% output the bode plots
loglog(imag(om),abs(gs))
xlabel('Omega')
ylabel('magnitude - open loop')
title('Magnitude bode plot')
disp('press a key to continue')
pause

semilogx(imag(om),angle(gs))
xlabel('Omega')
ylabel('phase - open loop')
title('phase bode plot')