Abstract. This paper presents mesh-free procedures for solving linear di erential equations (ODEs and elliptic PDEs) based on Multiquadric (MQ) Radial Basis Function Networks (RBFNs). Based on our study of approximation of function and its derivatives using RBFNs that was reported in an earlier paper (Mai-Duy and Tran-Cong, 1999), new RBFN approximation procedures are developed in this paper for solving DEs, which can also be classi ed into two t ypes: a direct (DRBFN) and an indirect (IRBFN) RBFN procedure. In the present procedures the width of the RBFs is the only adjustable parameter according to a (i) = d (i) , where d (i) is the distance from the ith centre to the nearest centre. The IRBFN method is more accuarte than the DRBFN one and experience so far shows that can be chosen in the range 7 10 for the former. Di erent c o m binations of RBF centres and collocation points (uniformly and randomly distributed) are tested on both regularly and irregularly shaped domains. The results for a 1D Poisson's equation show that the DRBFN and the IRBFN procedures achieve a norm of error of at least O(1:0e ; 4) and O(1:0e ; 8), respectively, with a centre density of 50. Similarly, the results for a 2D Poisson's equation show t h a t t h e DRBFN and the IRBFN procedures achieve a norm of error of at least O(1:0e;3) and O(1:0e ; 6), respectively, with a centre density o f 1 2 12.