ne(a,b) = 0 for (a=b)
= 1 for (a not equal to b)
H(x) is the heavyside function
I have an 1D domain in x and is divided into no. of cells. and so on are the cell boundaries with i=1,2…….
is the representative point of the cell (here its the mid point of the cell)
Taking for example =5, how can I write the code in Matlab so that I get separate expressions for for i=1,2,3,4,5?