function bai12
n=input('Nhap n:');
m=input('Nhap m:');
w=linspace(0,100,1000);
B=0;
for i=1:(m+1)
b(i)= input(['Nhap b',num2str(i-1),':']);
B=b(i)*(j*w).^(i-1)+B;
end
A=0;
for i=1:(n+1)
a(i)= input(['Nhap a',num2str(i-1),':']);
A=a(i)*(j*w).^(i-1)+A;
end
G=A./B
R=real(G);
Q=imag(G);
C=(R.^2+Q.^2).^(1/2)
phi=atan(Q./R)
subplot(2,1,1)
plot(w,20*log(C))
subplot(2,1,2)
plot(w,phi)