I think you can derive the backprop equation for softmax part. For maxpool layer, set the activations to 0 for non-max values. For example, if d_out = 10, so in maxpool pass the value 10 as gradient for the maximum value and 0 for others.
For backprop through conv layers, you need to do deconv operation.