clear;

%プラントの定義
m=0.1;		%質量m
tuta=0.5392;	%減衰比ζ
w=3.8623;	%固有振動数ω

k=m*w^2;	%バネ定数k
c=m*2*w*tuta;	%減衰係数c

%%%%%%%%プラントの伝達関数表現
pnum=[0 0 1/m];
pden=[1 c/m k/m];
pole_p=roots(pden)

%%%%%%%%%グラフのプロット
plot(real(pole_p),imag(pole_p),'k*')
grid
xlabel('Real number'),ylabel('Imaginary number')
title('placement of pole')

