clear all
close all
clc
format short eng

z=roots([1 -4])
p=roots([1 5 6])