We present non-LTE analyses of four M33 early B-supergiant stars and five Galactic counterparts. This is the first time that B supergiants beyond the Magellanic Clouds are analyzed by means of detailed Non-LTE techniques. Stellar temperatures, gravities, microturbulences and Si abundances are derived for all objects using the Si ionization equilibrium and the Balmer line wings. O abundances are then also derived. Important approximations made during the calculations are described, and their influence on the results is analyzed. It is found that these approximations have no significant effect on the results at any microturbulence. We find that a strict differential analysis is needed to detect abundance differences. Thus we compare results line by line in M33 and Galactic stars of stellar parameters as similar as possible. Three of the four M33 stars turned out to be O deficient as compared to their Galactic counterparts, and only one, close to the center of M33 (M33 1054) is found to be moderately O enriched. From these differential analyses we find that our data are compatible with a radial O gradient in M33 as that derived from HII region data: we obtain -0.19+/-0.13 or -0.20+/-0.07 dex/kpc, depending on whether B133 is included or not. Our data are also consistent with other possibilities such as a steep increase of the O abundance in the inner region (at projected distances less than 9 arcmin from the center of M33), followed by a flat O abundance profile towards the outer parts of M33. Si shows the same pattern, and it is shown that Si and O correlate well, as expected for alpha-elements, supporting then the high value of the O abundance gradient in M33 as compared to the Milky Way and other nearby spiral galaxies.