U1=[p^2+a1*p  a2    b5      0
    a6*s      b2    b3      0
    -s        s+a4  0       p
    ie*iepsi  0     Tf*s+1  0];


U2=[s^2+a1*s  a2    b5      b6*s
    a6*s      b2    b3      s^2+b1*s
    -s        s+a4  0       -b7*s-b4
    ie*iepsi  0     Tf*s+1  -(myue*s+ie)*(Tf*s+1)];


F=det(U1)/det(U2);