How to extract pol & zero data from amplifier's data sheet?
I'm not home right now, but here's a basic, ideal way to make it in LTspice, as a subcircuit:

.subckt ideal_opa548 IN- IN+ OUT params: Aol=106 Gbw=850k Rout=50
G1 0 out IN+ IN- {10**(Aol/20)/Rout}
C1 out 0 {10**(Aol/20)/(2*pi*Gbw*Rout)} Rpar={Rout}
.ends ideal_opa548


I didn't check the datasheet for Rout, and only ogled the gain-bandwidth, so add your own. Here's how it would look like in the schematic editor:



The response is with an 8\$\Omega\$ (load) and 8G\$\Omega\$ (no load). In .AC it looks fine, but in .TRAN you'll encounter problems with output swing (not limited). For this you can start the long way to complication by adding an appropiate "totem-pole" with SW, which will draw current, then an appropiate driver, input resistance/impedance, voltage limiting, ...