@chartl I had put this on hold and came back to it today.
I implemented the above and it solves. Is there a way to specify the constraint that ϕ=0 when x=0? I tried excluding ϕ where x=0 (by multiplying it by 0) but it still gives a non-zero value to ϕ and relatively high (0.06).
Thanks in advance!