Hello everyone!
This is probably the first time I'm posting here. I want to compute the spatiotemporal temperature distribution of a liquid droplet which is based on a Fourier series solution to the transient heat conduction equation in spherical coordinates. The solution requires us to compute the eigenfunctions and Fourier coefficients along with a few other terms before we can get to the final expression. These results have already been published, the solution algorithm, results from the paper and the references are given in the algorithm.doc file as attached. The code I have solves all of that but it is not able to replicate the original results even though all the inputs are the same. I haven't been able to figure out where I'm erring. Please let me know if you need more details. I would be grateful if anyone could pitch in to help. I have attached the files you may need during the discourse.
Regards Zuhaib
Attachments: