Presentamos un marco computacional para modelar una única vesícula inextensible impulsada por la fuerza de Helfrich en un fluido extracelular de Carreau incompresible y no newtoniano. La membrana de la vesícula se captura con una estrategia de conjunto de nivel. La restricción local de inextensibilidad se relaja introduciendo una penalización que permite ahorros computacionales y facilita la implementación. Una aproximación de elementos finitos de Galerkin de alto orden permite cálculos precisos de la fuerza de la membrana con derivadas de alto orden. La discretización temporal se basa en la composición doble del esquema de Euler hacia atrás de un paso, mientras que el tamaño del paso de tiempo se controla de forma flexible utilizando una estimación del error de integración temporal. Se presentan ejemplos numéricos con especial atención a la validación y evaluación de la relevancia del modelo en términos de importancia fisiológica. Se obtienen tasas de convergencia óptimas de la discretización temporal.