O presente trabalho teve como meta desenvolver uma rotina computacional para obter a soluÃÃo numÃrica de equaÃÃes diferenciais parciais atravÃs do MÃtodo sem malhas. Esse mÃtodo vem sendo aplicado na engenharia, visto que nos casos prÃticos muitas vezes encontram-se descontinuidades ou condiÃÃes de contorno que impossibilitam a soluÃÃo analÃtica, e muitas vezes dificultam bastante as soluÃÃes numÃricas de equaÃÃes diferenciais parciais. No mÃtodo utilizado, calcula-se as equaÃÃes a partir do mÃtodo numÃrico de Kansa, onde toma-se uma FunÃÃo de Base Radial (RBF) que vai calcular a matriz a partir da qual chegaremos à soluÃÃo analÃtica. A RBF utilizada foi a multiquadrÃtica. Essa funÃÃo irà depender de um parÃmetro de forma âcâ. Esse parÃmetro nÃo tem um valor definido, e o objetivo da rotina à encontrar um valor para esse parÃmetro, de forma que a soluÃÃo analÃtica se aproxime o mÃximo possÃvel da soluÃÃo numÃrica. Para se calcular esse valor do âcâ, a rotina irà calcular os valores do resÃduo para o domÃnio, e os valores do resÃduo para o contorno, para cada valor de âcâ em um intervalo determinado. ApÃs calculados esses valores, o programa compara os mesmos e fornece o ponto em que eles se aproximam mais. Nesse ponto à encontrado o parÃmetro de forma âcâ otimizado, e consequentemente a soluÃÃo numÃrica da equaÃÃo proposta. / This study was aimed to develop a computational routine for the numerical solution of partial differential equations using the method without mesh. This method has been applied in engineering, since in practical cases often are discontinuities or boundary conditions that preclude the analytic solution, and often quite difficult numerical solutions of partial differential equations. In the method used, it is estimated the equations from the numerical method of Kansas, it takes it a Radial Basis Function (RBF) that will calculate the matrix from which will come to the analytical solution. The RBF used was multiquadrÃtica. This function will depend on a parameter in a "c". This parameter does not have a defined value, and the routine goal is to find a value for this parameter, so that the analytical solution is closely match the numerical solution. To calculate this value of "c", the routine will calculate the residue values ​​for the domain, and the residue values ​​for the outline, for each value of "c" in a certain range. After these calculated values, the program compares them and supplies the point where they are closer to. This point is found as the parameter "c" optimized, and therefore the numerical solution of the equation proposed.