Modeling of a Type I compensator in matlab

I am trying to get the bode plot of a Type I compensator with a non ideal opamp. The schematic would look like an integrator where a 3.3nF cap is on the feedback of the opamp, a 7k resistor is in series with the inverting input and the non inverting input is grounded.

The opamp transfer function is:

               3.158e13
        ---------------------------
  s^2 + 6.289e06 s + 3.948e10

and I know the feedback transfer function is 

       1
  ----------
  2.31e-05 s

since we have 1/( R*sC)

My question is how to combine the two transfer functions in order to get the system response in a bode plot? I am simulating the circuit in LT spice and so far I haven't been able to match the frequency response of what I have tried in MATLAB and what I get in spice.

Thanks for the help

Hi

I'm sure you know how to work with transfer functions, http://www3.ul.ie/tilde_accs/mlc/www/support/Loughborough website/chap20/20_7.pdf.
Check this example, http://lpsa.swarthmore.edu/Bode/BodeExamples.html