El objetivo de este trabajo es implementar en el programa FEMLAB ecuaciones matemáticas que encontramos en biología, en especial las que se conocen como ecuaciones de Turing o de reacción difusión, que son ecuaciones diferenciales parciales no lineales que describen un mecanismo de formación de patrones semejantes a los que podemos observar en la naturaleza. Primeramente se presentan los aspectos básicos de la teoría de las ecuaciones de Turing, incluyendo el análisis lineal de estabilidad asintótica de las ecuaciones para obtener los diferentes parámetros que dan lugar a la formación de un patrón. En segundo lugar implementarnos y resolvemos algunas de estas ecuaciones con el programa FEMLAB, el cual esta diseñado para modelar y simular fenómenos físicos fácilmente; nos permite introducir ecuaciones diferenciales parciales de manera independiente y personalizada o el uso de modelos ya determinados. Una gran ventaja que observaremos es que con el uso de este programa podemos utilizar diferentes geometrías en el dominio, combinarlas en los sistemas y utilizar distintas dimensiones del espacio en el mismo problema. Finalmente, exploramos soluciones de algunas ecuaciones de Turing en dominios complejos.