# I am trying to interpolate a function defined over a three-dimensional real space: f :

I am trying to interpolate a function defined over a three-dimensional real space:
$f:{R}^{3}\to R\phantom{\rule{0ex}{0ex}}\left(x,y,z\right)\to f\left(x,y,z\right)$
Let assume I have ${N}_{1}{N}_{2}{N}_{3}$ points in the space which form my grid for this interpolation, and the multivariate series
$F\left(x,y,z\right)=\sum _{a=0}^{{N}_{1}-1}\sum _{b=0}^{{N}_{2}-1}\sum _{c=0}^{{N}_{3}-1}{C}_{abc}{x}^{a}{y}^{b}{z}^{c}$
is the chosen interpolator. In order to find the coefficients I should form the following sets of equations:

and then solve for the coefficients. However, for this I first need to write the above set of linear equations in the standard form
$A\stackrel{\to }{x}=\stackrel{\to }{b}$
wherein, A is the matrix of coefficients, $\stackrel{\to }{x}$ is the vector of unknowns, and $b=\left\{{f}_{i}{\right\}}_{1}^{{N}_{1}{N}_{2}{N}_{3}}$ is the known vector. For this to be done I would require to expand the multivariate power series in the form of a single variable series, that is,
$\sum _{a=0}^{{N}_{1}-1}\sum _{b=0}^{{N}_{2}-1}\sum _{c=0}^{{N}_{3}-1}{C}_{abc}{x}_{i}^{a}{y}_{i}^{b}{z}_{i}^{c}=\sum _{m=0}^{\left({N}_{1}-1\right)\left({N}_{2}-1\right)\left({N}_{3}-1\right)}{C}_{m}{\eta }_{i}^{m}$
wherein ${\eta }_{i}={\eta }_{i}\left({x}_{i},{y}_{i},{z}_{i}\right)$ and probably $m=abc$. OF course the expansion needs to be nontrivial and useful.
Is it possible at all? Any suggestion to find the coefficients more practically?
Regards, owari
UPDATE.
Maybe it appears that the most natural way for solving this problem is comprised of the following steps:
1. first solve for the coefficients of

2. then solve for the coefficients of

3. and finally solve for the coefficients of

which gives the coefficients we were looking for. However, this way, the first two steps will be working with function-type coefficients instead of numerical coefficients and that will prevent efficient usage of the available codes in numerical analysis. Solving for each set of coefficients at each step for different grid points would also inevitably increase the number of equations drastically, so any better suggestion to solve for this problem?
You can still ask an expert for help

• Questions are typically answered in as fast as 30 minutes

Solve your problem for the price of one coffee

• Math expert for every subject
• Pay only if we can solve it

gudstrufy47j
Ok, if the function that is to be interpolated can be factorized independently in different variables, like $f\left(x,y,z\right)={f}_{x}\left(x\right){f}_{y}\left(y\right){f}_{z}\left(z\right)$, the solution would be easy as the triple sum would be easily written as the multiplication of three sums and the equations can be written in the form that follows:

and finally we would have: ${C}_{abc}={{C}_{x}}_{a}{{C}_{y}}_{b}{{C}_{y}}_{b}$. This is good and efficient, however it works for only a limited class of functions, any idea to extend it to a more general case?