The First Integral Method for (3+1) dimensional modified KdV-Zakharov-Kuznetsov and Hirota equations Mustafa Inc, Bulent Kılıç and Yavuz Uğurlu [email protected] Firat University, Elazığ We use first integral method to find traveling wave solutions of (3+1) dimensional modified KdV-Zakharov-Kuznetsov and Hirota equations in this paper. This method is based on commutative algebra. This is effective method for finding solutions of nonlinear partial differential equations. Keywords: First integral method, (3+1) dimensional modified KdV-Zakharov-Kuznetsov equation, Hirota equation, traveling wave solutions. References  Parkes, E. J., Duffy, B. R., An automated tanh-function method for finding solitary wave solutions to non-linear evolution equations, Computer Physics Communications, 98 (1996) 288-300.  Wang, M., Li, X. and Zhang, J., The G -expansion method and travelling wave G solutions of nonlinear evolution equations in mathematical physics, Physics Letters A, 372 (2008) 417-423.  Fu, Z., Liu, S., Liu, S. and Zhao, Q., New Jacobi elliptic function expansion and new periodic solutions of nonlinear wave equations, Physics Letters A, 290 (2001) 72-76.  Drazin, P. G. and Jhonson, R. S., Solitons: An Introduction. Cambridge University Press, 1989, Cambridge.  X. B. Hu, W. X. Ma, Application of Hirota's bilinear formalism to the Toeplitz latticesome special soliton-like solutions. Physics Letters A, 293 (2002) 161-165.