Hey guys,
I was given a parabolic differential equation with the boundary condition of dT/dR = 0 when R = 0 but I am not quite sure how I am suppose to incorporate that into my Thomas algorithm/Crank-Nicolson program.
I currently have it in the program as u1 = du/dr but it tells me the value must be scalar. When I just get rid of it, there are complex answers in my solution.
Please help?