The mode you're interested in is a quasi-TEM microstrip mode. There's not going to be any [realistically simple] analytical expression for the fields, sorry. You could spend the next month deriving some more exacting expressions, but what's the point?
The thesis uses approximations (such as assuming only the z-component of the fields), which yields relatively straightforward and somewhat accurate expressions. In real life engineering, these are often all we can really hope for -- I suggest that your results are sufficiently accurate for almost all purposes.