From: Iron on
Hello people

I am trying to create a programm that computes the heat and temperature distribution in a cylinder that gives out heat via conduction, with pdepe. The tricky thing is, that the cylinder consists of 4 layers of different heat capacities, densities and thermal conductivities. So I used a sum of 4 heaviside functions to define the parameter "c" of the problem. The rest of the parameters "f" and "s" are already defined.

But, when I see the results it seems like Matlab didnt take the heaviside functions into account, as if there was a constant heat capacity, density and thermal conductivity throughout the cylinder...
What should I do? Anydody has an idea?

Thanks alot!
Iron89