f1='(V3-V1)/R+(V3-V2)/(s*L)=0'
f2='s*C*V2+(V2-V3)/(s*L)=0'
rV=solve(f1,f2,'V2,V3');
rV.V2
rV.V3
P=rV.V2/'V1'
nP=subs(P,{'R','C','L'},{1e3,22.5e-9,11.25e-3})
[num, den]=numden(nP)
format short e
A=sym2poly(den)
B=sym2poly(num)/A(1)
A=A/A(1)
