x=1:1033; y= (x.^4 + 0.6489*x.^3 + 0.5417*x.^2 + 0.5605*x + 0.9077)./(x.^3 + 0.6489 *x.^2 + 0.5389*x + 0.677); plot(x,y) grid on