<s> """<STR_LIT>""" <EOL> from __future__ import division <EOL> from PyDSTool import * <EOL> from matplotlib import pyplot as plt <EOL> from numpy . linalg import eig <EOL> def make_A ( tau ) : <EOL> """<STR_LIT>""" <EOL> return array ( [ [ - <NUM_LIT:1.> / <NUM_LIT:20> , <NUM_LIT:0> , <NUM_LIT:0> , - <NUM_LIT:1.> / <NUM_LIT:20> ] , <EOL> [ <NUM_LIT:1.> / tau , - <NUM_LIT:1.> / tau , <NUM_LIT:0> , <NUM_LIT:0> ] , <EOL> [ <NUM_LIT:0> , <NUM_LIT> / <NUM_LIT:50> , - <NUM_LIT:1.> / <NUM_LIT:50> , <NUM_LIT:0> ] , <EOL> [ <NUM_LIT:0> , <NUM_LIT:0> , <NUM_LIT:1.> / tau , - <NUM_LIT:1.> / tau ] ] ) <EOL> cols = [ '<STR_LIT:r>' , '<STR_LIT:g>' , '<STR_LIT:k>' , '<STR_LIT:b>' ] <EOL> plt . plot ( [ - <NUM_LIT> , <NUM_LIT> ] , [ <NUM_LIT:0> , <NUM_LIT:0> ] , '<STR_LIT:k>' ) <EOL> plt . plot ( [ <NUM_LIT:0> , <NUM_LIT:0> ] , [ - <NUM_LIT> , <NUM_LIT> ] , '<STR_LIT:k>' ) <EOL> plt . xlabel ( '<STR_LIT>' ) <EOL> plt . ylabel ( '<STR_LIT>' ) <EOL> plt . show ( ) <EOL> Hopf_found = False <EOL> taus = linspace ( <NUM_LIT> , <NUM_LIT:15> , <NUM_LIT:100> ) <EOL> t = plt . text ( - <NUM_LIT> , <NUM_LIT> , '<STR_LIT>' % <NUM_LIT> , fontsize = <NUM_LIT:20> ) <EOL> plt . title ( '<STR_LIT>' ) <EOL> for i , tau in enumerate ( taus ) : <EOL> evals = eig ( make_A ( tau ) ) [ <NUM_LIT:0> ] <EOL> for j , e in enumerate ( evals ) : <EOL> re = real ( e ) <EOL> im = imag ( e ) <EOL> if re > <NUM_LIT:0> and not Hopf_found : <EOL> print "<STR_LIT>" % ( tau , i ) <EOL> Hopf_found = True <EOL> plt . plot ( re , im , cols [ j ] + '<STR_LIT:.>' ) <EOL> plt . xlim ( [ - <NUM_LIT> , <NUM_LIT> ] ) <EOL> plt . ylim ( [ - <NUM_LIT> , <NUM_LIT> ] ) <EOL> t . set_text ( '<STR_LIT>' % tau ) <EOL> plt . draw ( ) </s>
