curve intersection code snippets in matlab