Type debug to Turn On Debugging: debug Specify Debugging Destination, screen or file: screen List Functions to Debug - or all - end ends list Names following all will not be debugged - end ends :clock : end Enter NX: 150 Enter NY: 150 Enter NZ: 1 Enter dx: 0.6 Enter dt: 0.25 Enter final time: 600.0 Enter intracellular conductivity, (si_x,si_y,si_z): 0.1 0.1 0.0 Enter extracellular conductivity, (se_x,se_y,se_z): 0.1 0.1 0.0 Two ionic models are available FHN (f) 3V-SIM (3) Enter a model (f(default),3): 3 Enter the initial value for v: 0.9 Enter the initial value for w: 0.4 Enter tau_vp: 3.33 Enter tau_v1m: 19.6 Enter tau_v2m: 1000.0 Enter tau_wp: 667.0 Enter tau_wm: 11.0 Enter tau_d: 0.25 Enter tau_o: 8.3 Enter tau_r: 50.0 Enter tau_si: 45 Enter k: 10.0 Enter Vcsi: 0.85 Enter Vc: 0.13 Enter Vv: 0.055 Three problems are available Monodomain (m) 1d Bidomain (1) 2d or 3d Bidomain (b) Enter a problem (m,1,b(dflt)): b Enter the resistance: 1.0 Enter the surface volume ratio: 1.0 Four solution methods are available. Crank_Nicolson (c) Semi_Implicit (i) Split_Dense (s) Split_ADI (t) Enter a method (c,i,s or t): t Choose an initial condition. Plane Wave starting from left (p) Radial wave starting from a corner (c) Radial wave starting from the center (s) Enter the choice (p,c, or s): p Choose a method for stimulation. Cutting half (h) Stimulting the center (c) Enter a method for stimulation: h Enter the time for the stimulus: 50 Enter the strength of the stimulus: 0.8 Enter the radius of the stimulus: 2