You should use the intersection point as parameters, then perform global fit in Origin.
Assume that the intersection point is (x0, y0), and polynomial functions are parabolic, then fitting functions can be defined as follows:
f(x)=k*(x-x0)+y0;
g(x)=(x-x0)*(a1+b1*x)+y0;
h(x)=(x-x0)*(a2+b2*x)+y0;
Parameters are x0, y0, k, a1, b1, a2 and b2.
For polynomial functions of high order, you can define it in a similar way.
Sam
OriginLab Technical Services