This work presents the development of a computational code in the MATLAB language for geometrically nonlinear analysis of plane trusses. This development was based on a formulation for geometrically nonlinear analysis of plane trusses through the direct stiffness method and the implementation of the Newton-Raphson method with load increments for the solution of the nonlinear equations system that represents the load-displacement relations derived by the aforementioned formulation. The code was then numerically tested through examples found in the literature and comparison of results. Good agrément was noticed between the results of nodal displacements and axial forces in the bars obtained through the code developed by the author and the results available in the literature. In examples 1 and 2, the trusses analyzed underwent the instability effect known as "snap-through" and it was possible to observe the structural behaviour and the influence of such effect on the number of iterations needed for analysis convergence. In example 4, the validity of the linear analysis for structures that undergo small displacements could be noted by comparing the results with a geometrically nonlinear analysis.Este trabalho apresenta o desenvolvimento de um código computacional na linguagem MATLAB para a análise geometricamente não linear de treliças planas. Tal desenvolvimento foi baseado em uma formulação para análise geometricamente não linear de treliças planas através do método da rigidez direta e na implementação do método de Newton-Raphson com incrementos de carga para solução do sistema de equações não lineares que representam as relações carga-deslocamento derivadas da formulação. O código foi, então, numericamente testado através de exemplos encontrados na literatura técnica e comparação de resultados. Observou-se boa concordância entre os resultados de deslocamentos nodais e forças axiais nas barras obtidos através do código desenvolvido pelo autor e os resultados disponíveis na bibliografia consultada. Nos exemplos 1 e 2, as treliças analisadas sofreram o efeito de instabilidade conhecido como “snap-through” e foi possível observar o comportamento das estruturas e a influência de tal efeito sobre o número de iterações necessárias para convergência da análise. No exemplo 4, pôde-se perceber a validade da análise linear para estruturas que sofrem pequenos deslocamentos através da comparação dos resultados com uma análise geometricamente não linear.São Cristóvão, SE