代码:
b = [1, 0, 0, 0, 16+1/16, 0, 0, 0, 1]; broots = roots(b) B1 = real(poly([broots(1), broots(2), broots(5), broots(6)])) B2 = real(poly([broots(3), broots(4), broots(7), broots(8)]))
运行结果:
框图为