First draft of code; run from Bioheat_script.m