In view of the problem in which most target bodies are continuously changing and the distribution is complex in the actual geological structures,2-D magnetotelluric forward modeling is simulated combined with the advantage of meshfree method where parameter loading is convenient,calculation accuracy is high and adaptive analysis is convenient.The shape function of the continuous medium model is constructed by the radial basis function,the equivalent linear equations are derived and the corresponding program is compiled.The optimal values of the shape parameters are given by numerical experiments.The correctness of meshfree method is verified by a horizontal continuous variation layered model.The root mean square error of the calculation results is not more than 0.36%,and its accuracy is superior to the finite element method.The electromagnetic response of the horizontal layer with the continuous variation of electrical conductivity and the uniform block model is discussed,and a graben terrain model and a reservoir water flooding model whose conductivity is continuously variable are calculated,respectively.The conclusions are given: the electromagnetic responses between the continuous variation model and uniform block model are significantly different,the use of the continuous variation model is helpful to improve the accuracy of inversion in the actual interpretation,and the observation method of the TM mode has better resolution for the inclination distribution of anomalous body.The meshfree method does not need the complicated mesh generation,and it is more suitable to calculate the response of continuous variation and complex distribution models.Therefore,meshfree method is expected to be a new and robust method for the numerical simulation of complex medium and complex structures.