[en] A new MHD-equilibrium/current-drive analysis code was developed to analyse the high beta tokamak equilibria consistent with the beam driven current profiles. In this new code, the critical beta equilibrium, which is stable against the ballooning mode, the kink mode and the Mercier mode, is determined first using MHD equilibrium and stability analysis codes (EQLAUS/ERATO). Then, the current drive parameters and the plasma parameters, required to sustain this critical beta equilibrium, are determined by iterative calculations. The beam driven current profiles are evaluated by the Fokker-Planck calculations on individual flux surfaces, where the toroidal effects on the beam ion and plasma electron trajectories are considered. The pressure calculation takes into account the beam ion and fast alpha components. A peculiarity of our new method is that the obtained solution is not only consistent with the MHD equilibrium but also consistent with the critical beta limit conditions, in the current profile and the pressure profile. Using this new method, β ∼ 21 % bean and β ∼ 6 % D-type critical beta equilibria were scanned for various parameters; the major radius, magnetic field, temperature, injection energy, etc. It was found that the achievable Q value for the bean type was always about 30 % larger than for the D-type cases, where Q = fusion power/beam power. With strong beanness, Q ∼ 6 for DEMO type tokamaks (∼500 MWth) and Q ∼ 20 for power reactor size (4.5 GWth) are achievable. On the other hand, the Q value would not exceed sixteen for the D-type machines. (author)