def dsdt( )return def (xt)return parameters can be set as beta - linspace( ln+ dx [ [ zeros( + u_ zeros( + u_ [ ( u_ [ : let us use : we can now call ode_fe and then make an animation on the screen to see how xtdevelops in timet time clock(from ode_system_fe import ode_fe ut ode_fe(rhsu_ dtt= * * time clock(print 'cpu time fs( make movie import os os system('rm tmp_png'import matplotlib pyplot as plt plt ion( [ ,:lines plt plot(xyplt axis([ [ ] [- ] ( )+ ]plt xlabel(' 'plt ylabel(' ( , )'counter plot each of the first framesthen increase speed by change_speed for in range( shape[ ])print [iplot true if <change_speed else = lines[ set_ydata( [ ,:]if change_speedplt legend([' = xt[ ]]elseplt legend([' = ft[ ]]plt draw(if plotplt savefig('tmp_% pngcountercounter + #time sleep(