Computing the Jacobian matrix for a solar system simulation

[doonzy]

Hello physicists!

I'm a comp sci student and I am trying to graphically model a simplified version of the solar system as part of a programming exercise. In order to apply the gravitational forces to the planets, I need to compute the Jacobian matrix as it relates to two particles (planetary objects). My (limited) understanding is that the resulting matrix will be 9x9, but I am unsure how it is constructed exactly. This is what I have come up with so far:

$$J = \begin{pmatrix}<br /> \frac{\partial F_1}{\partial r_1} & \frac{\partial F_1}{\partial r_2} & \frac{\partial F_1}{\partial r_3}\\ <br /> \frac{\partial F_2}{\partial r_1} & \frac{\partial F_2}{\partial r_2} & \frac{\partial F_2}{\partial r_3}\\ <br /> \frac{\partial F_3}{\partial r_1} & \frac{\partial F_3}{\partial r_2} & \frac{\partial F_3}{\partial r_3}<br /> \end{pmatrix}$$

where $$\frac{\partial F_i}{\partial r_i}$$ is 3x3, $$F$$ is the gravitational force $$F = \frac{Gm_1m_2}{r^2}$$ and $$r$$ is the respective dimension component $$(x, y, z)$$.

Some clarification would be much appreciated :-)

Thanks.

 

[fresh_42]

Well, $r^2=\sqrt{r_1^2+r_2^2+r_3^2}$, hence $F\, : \,\mathbb{R}^3 \longrightarrow \mathbb{R}$ which means it has only one component and your Jacobian matrix is a gradient, i.e. of size $1 \times 3$. The derivative is the same as at school: differentiate with only one variable $r_i$ and consider all other values as constant.