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)