p=[-1.63601e-02+0.00000e+00i;
    -2.93921e-02+0.00000e+00i;
    -3.18408e-02+0.00000e+00i;
    -3.97838e-02+0.00000e+00i;
    -3.97866e-02+0.00000e+00i;
    -6.15560e-02+0.00000e+00i;
    -7.04742e-02+0.00000e+00i;
    -1.64334e-03-5.68434e+06i;
    -1.64334e-03+5.68434e+06i;
    -4.06044e+06-8.02865e+07i;
    -4.06044e+06+8.02865e+07i;
    -2.55686e+02+2.21772e+09i;
    -2.55686e+02-2.21772e+09i;
    -6.21109e+03-2.31261e+09i;
    -6.21109e+03+2.31261e+09i;
    -2.30225e+02-2.64997e+09i;
    -2.30225e+02+2.64997e+09i;
    -2.83489e+03-3.04281e+09i;
    -2.83489e+03+3.04281e+09i;
    -2.10592e+03-3.06838e+09i;
    -2.10592e+03+3.06838e+09i;
    -9.66426e+02-3.29534e+09i;
    -9.66426e+02+3.29534e+09i;
    -2.59758e+02-3.45230e+09i;
    -2.59758e+02+3.45230e+09i;
    -1.35539e+02-3.79774e+09i;
    -1.35539e+02+3.79774e+09i;
    -2.93730e+01-4.61102e+09i;
    -2.93730e+01+4.61102e+09i;
    -4.42050e+03-6.30144e+09i;
    -4.42050e+03+6.30144e+09i;
    -3.80191e+02-6.54472e+09i;
    -3.80191e+02+6.54472e+09i;
    -7.80846e+02-6.60252e+09i;
    -7.80846e+02+6.60252e+09i;
    -5.74363e+01-6.92779e+09i;
    -5.74363e+01+6.92779e+09i;
     -3.97559e+02-7.04464e+09i;
      -3.97559e+02+7.04464e+09i;
      -1.32259e+03-1.16857e+10i;
      -1.32259e+03+1.16857e+10i]

  
  z=[-2.52933e-04+0.00000e+00i;
      -3.12784e-02-1.78701e-03i;
      -3.12784e-02+1.78701e-03i;
      -3.35785e-02+0.00000e+00i;
      -3.94756e-02+0.00000e+00i;
      -6.16305e-02+0.00000e+00i;
      -1.56380e-01+0.00000e+00i;
      -4.61387e-01+0.00000e+00i;
      4.67146e-01+0.00000e+00i;
      -4.08025e+06-8.04815e+07i;
      -4.08025e+06+8.04815e+07i;
      -1.28022e-02-2.21092e+09i;
      -1.28022e-02+2.21092e+09i;
      -2.77470e-02-2.55409e+09i;
      -2.77470e-02+2.55409e+09i;
      -2.62383e-02-3.06692e+09i;
      -2.62383e-02+3.06692e+09i;
      -1.16794e+08-3.07766e+09i;
      -1.16794e+08+3.07766e+09i;
      1.16793e+08-3.07766e+09i;
      1.16793e+08+3.07766e+09i;
      -2.69539e-02-3.43574e+09i;
      -2.69539e-02+3.43574e+09i;
      -1.59476e-02-4.24387e+09i;
      -1.59476e-02+4.24387e+09i;
      1.42797e-06-5.03292e+09i;
      1.42797e-06+5.03292e+09i;
      -3.29387e-03-6.51065e+09i;
      -3.29387e-03+6.51065e+09i;
      -3.47367e+02-6.62013e+09i;
      -3.47367e+02+6.62013e+09i;
      3.75036e+02-6.88285e+09i;
      3.75036e+02+6.88285e+09i;
      -7.23296e-03-7.00346e+09i;
      -7.23296e-03+7.00346e+09i;]
k=8.54677e+59;
sys=zpk(z,p,k)
bode(sys)