Este trabalho trata do desenvolvimento e teste da rotina QUASILIN, utilizada na obtencao de solucoes numericas de problemas de condicao de contorno em dois ou mais pontos. Essa rotina foi elaborada tendo como base o metodo da quasilinearizacao, que e um metodo numerico do tipo indireto, onde o problema original nao linear e transformado em uma serie de problemas lineares. Um algoritmo auxiliar que, a partir dos valores iniciais das variaveis de estado, faz a escolha otima dos valores iniciais das variaveis adjuntas, tambem e apresentado. A funcionalidade da rotina e demonstrada atraves de solucoes de problemas tipicos da area aeroespacial, tais como: calculo de trajetoria de reentrada atmosferica, de minimo aquecimento, para um veiculo espacial do tipo "space shuttle"; minimizacao de energia de operacao de um braco manipulador e a maximizacao de altitute para um foguete de sondagem de estagio unico.