how to calculate the matrices in finite element analysis

you have to check for when |(x-x_j)/h|<=1, obviouslt when x is in [x_j-h,x_j+h] now x also should be between [x_i-h,x_i+h]
Now check for the next cases:
1. both intervals are inside the interval [0,1]
2. one interval its lower endpoint is smaller than 0 and it's higher endpoint is in [0,1], the second interval is totally inside [0,1] (the same symmetric case by changing the intervals).
3. both intervals have their lower endpoints below 0 and their highest endpoints inside [0,1]. (and the opposite case where the highest endpoints are above 1 and the lower endpoints inside [0,1]).
4. like case 2 but changing 0 with 1 and lower endpoint with higher endpoint.

In each case you should determine the integral limits, by deciding on lower limit to be max{0,x_j-h, x_i-h} and upper limit to be min{1,x_j-h,x_i-h}.

Sorry I can't help you more than that, I am kinda busy myself.