Find particular solution given initial condition (#579)

Find particular solution given initial condition (#579)

Substitute x0,y0x_0, y_0 after integration to determine constant.