Simulations of nematic-isotropic transition of liquid crystals in two dimensions are performed using an O(2) vector model characterized by non linear nearest neighbour spin interaction governed by the fourth Legendre polynomial ...