AUTHORS: Ariane Ernst
CONTACT: Ariane Ernst https://github.com/ahjanee
DATE: June 2023

This code was used for the computations in our article 'Model reduction for Ca2+-induced vesicle fusion dynamics'.

PREREQUISITES:
This code was written using Python 3.6.13. You will need the sympy library.

CODE STRUCTURE AND USAGE:
We provide two executables and one parameter file (ks_and_paras). The parameter file contains the parameters used 
in Kobbersmed et al. 2020 as described in the article. The file rate_expressions.py provides the symbolic expressions
for the reduced model's fusion and unpriming rate. You can insert the paramters by uncommmentin the marked section. 
The file matrix_inversion_and_xqss.py contains the symbolic computation of Eq.9 in the article, including an implementation
of the algorithm for tridiangular matrix inversion by Jia and Li (2015).
