MATLAB PDEs solver/animation

Nuevamente, haciendo uso de Claude AI, refinando un prompt y dando algunas indicaciones básicas de MATLAB estoy desarrollando una calculadora en MATLAB para visualizar las soluciones de las principales PDEs: ondas, calor y Laplace. En esta versión del código se permite la selección de varios items, como, por ejemplo, el tipo de condiciones de contorno Dirichlet, Neuman, periódicas, entre otras. Aún estas dos versiones que presento necesitan ser revisadas y, en estos días, publicaré la versión final.

La Versión 3 será el producto de mi programación (depurada por Claude Sonnet 4.6) inspirada en los dos productos anteriores, obviamente, revisando posibles errores de las versiones 1 y 2. Próximamente disponible.

Mejoras de la versión 3

1. second_deriv_vec vectorizado, elimina el bucle for interior, hasta ×10 más rápido para grids grandes

2. Solver CG real con matrices sparse — build_laplace_matrix ensambla el Laplaciano 5-puntos como sparse y lo resuelve con Gradiente Conjugado auténtico, en lugar del fallback SOR de V2

3. validate_parameters con assert informativo antes de llamar al solver

4. print_solver_info imprime dx, CFL estimado y advertencias de rigidez

5. Output functions con persistent correctamente inicializadas (ctr=0 al recibir flag='init')