The paper describes meshfree boundary point solution of problems based on Trefftz interpolation functions. The solution is a collocation method, which can use any type of Trefftz functions smooth enough for approximation of domain fields. Polynomial, Kelvin and Boussinesq functions are used in our formulations. A multi-domain boundary formulation with continuity of tractions and displacements over discrete points of the inter-domain surface is shown. The method is applied to linear elastostatics. 2D and 3D examples illustrate the formulation. Discussion about advantages and drawbacks is included. As the conservation equations (force and moment equilibrium) are not satisfied in integral sense, some benchmark tests are proposed for meshless formulations to obtain a view about the accuracy.