# # Written by: # -- # John L. Weatherwax 2005-08-14 # # email: wax@alum.mit.edu # # Please send comments and especially bug reports to the # above email address. # #----- k_1 = 1; k_2 = 0.1; k_3u = 10; #k_3u = 1; mu_q = 1; # process noise mean z = c( mu_q + k_3u, -k_1, 0, 0, -k_2 ) pr = polyroot( z ) pr abs(pr)