Edit: I have reworded the question.
This may be a basic question but I am having trouble figuring out the correct answer. I want to find a local coordinate chart that fits a d-dimensional submanifold in $\mathbb{R}^N$. I am given two points $p_1, p_2 \in \mathbb{R}^N$ and corresponding orthonormal bases $(\phi_1, \phi_2, \ldots, \phi_d)$, $(\tau_1, \tau_2, \ldots, \tau_d) \subseteq \mathbb{R}^N$ for their tangent spaces. I would like to find an algorithmic method for finding functions $(f_1(x_1,x_2,\ldots,x_d), f_2(x_1,x_2,\ldots,x_d),\ldots f_N(x_1,x_2,\ldots,x_d))$ that satisfy these conditions.