MHB How to Solve Non-Linear Equations of 3 Variables Using Newton-Raphson Method?

AI Thread Summary
The discussion focuses on solving a system of three non-linear equations using the Newton-Raphson method, specifically for the variables c, s, and q. The user expresses difficulty in calculating the derivatives of the functions involved and applying the Newton-Raphson iteration due to the interdependence of the variables. They seek guidance on how to derive the necessary derivatives and proceed with the iteration process. The conversation emphasizes the importance of forming a vector equation and utilizing the Jacobian matrix for the Newton-Raphson method. Overall, the thread highlights the challenges of applying this numerical method to multi-variable non-linear equations.
ra_forever8
Messages
106
Reaction score
0
The three non-linear equations are given by
\begin{equation}
c[(6.7 * 10^8) + (1.2 * 10^8)s+(1-q)(2.6*10^8)]-0.00114532=0
\end{equation}
\begin{equation}
s[2.001 *c + 835(1-q)]-2.001*c =0
\end{equation}
\begin{equation}
q[2.73 + (5.98*10^{10})c]-(5.98 *10^{10})c =0
\end{equation}
Using the Newton-Raphson Method solve these equations in terms of $c$,$s$ and $q$.=> It is really difficult question for me because i don't know very much about the Newton-Raphson Method and also these non-linear equations contain 3 variables.I have try by applying the Newton-Raphson method to each equations:-
\begin{equation}
f(c,s,q)=0= c[(6.7 * 10^8) + (1.2 * 10^8)s+(1-q)(2.6*10^8)]-0.00114532
\end{equation}
\begin{equation}
g(c,s,q)=0= s[2.001 *c + 835(1-q)]-2.001*c
\end{equation}
\begin{equation}
h(c,s,q)=0= q[2.73 + (5.98*10^{10})c]-(5.98 *10^{10})c
\end{equation}
now i guess i need to work out $f'(c,s,q), g'(c,s,q), h'(c,s,q)$ but i don't know how?and after working out $f'(c,s,q), g'(c,s,q), h'(c,s,q)$ . After that i think i need to use Newton-raphson iteration:$c_{n+1}= c_n - \frac{f(c,s,q)}{f'(c,s,q)}$ but the $f(c,s,q)$ and $f'(c,s,q)$ contains the $s$ and $q$.

Similarly, for $s_{n+1}= s_n - \frac{g(c,s,q)}{g'(c,s,q)}$ will have $g(c,s,q)$ and $g'(c,s,q)$ containing the $c$ and $q$.$q_{n+1}= q_n - \frac{h(c,s,q)}{h'(c,s,q)}$ will have $h(c,s,q)$ and $h'(c,s,q)$ containing the $c$.so am i not sure what to do please help me. to find the values of $c,s,q$.
 
Mathematics news on Phys.org
grandy said:
The three non-linear equations are given by
\begin{equation}
c[(6.7 * 10^8) + (1.2 * 10^8)s+(1-q)(2.6*10^8)]-0.00114532=0
\end{equation}
\begin{equation}
s[2.001 *c + 835(1-q)]-2.001*c =0
\end{equation}
\begin{equation}
q[2.73 + (5.98*10^{10})c]-(5.98 *10^{10})c =0
\end{equation}
Using the Newton-Raphson Method solve these equations in terms of $c$,$s$ and $q$.=> It is really difficult question for me because i don't know very much about the Newton-Raphson Method and also these non-linear equations contain 3 variables.I have try by applying the Newton-Raphson method to each equations:-
\begin{equation}
f(c,s,q)=0= c[(6.7 * 10^8) + (1.2 * 10^8)s+(1-q)(2.6*10^8)]-0.00114532
\end{equation}
\begin{equation}
g(c,s,q)=0= s[2.001 *c + 835(1-q)]-2.001*c
\end{equation}
\begin{equation}
h(c,s,q)=0= q[2.73 + (5.98*10^{10})c]-(5.98 *10^{10})c
\end{equation}
now i guess i need to work out $f'(c,s,q), g'(c,s,q), h'(c,s,q)$ but i don't know how?and after working out $f'(c,s,q), g'(c,s,q), h'(c,s,q)$ . After that i think i need to use Newton-raphson iteration:$c_{n+1}= c_n - \frac{f(c,s,q)}{f'(c,s,q)}$ but the $f(c,s,q)$ and $f'(c,s,q)$ contains the $s$ and $q$.

Similarly, for $s_{n+1}= s_n - \frac{g(c,s,q)}{g'(c,s,q)}$ will have $g(c,s,q)$ and $g'(c,s,q)$ containing the $c$ and $q$.$q_{n+1}= q_n - \frac{h(c,s,q)}{h'(c,s,q)}$ will have $h(c,s,q)$ and $h'(c,s,q)$ containing the $c$.so am i not sure what to do please help me. to find the values of $c,s,q$.

Write your system as a vector equation in a vector variable, you want a solution of the system:
$${\bf{f}}({\bf{x}})={\bf{0}}$$

The NR iteration for this is:

$${\bf{x}}_{n+1}={\bf{x}}_n-[{\rm{J}}({\bf{x}}_n)]^{-1}{\bf{f}}({\bf{x_n}})$$

where ${\rm{J}}({\bf{x}})$ is the matrix with $i.j$ th element equal to $ \dfrac{\partial {\bf{f}}_i}{\partial{\bf{x}}_j}$
 
Thread 'Video on imaginary numbers and some queries'
Hi, I was watching the following video. I found some points confusing. Could you please help me to understand the gaps? Thanks, in advance! Question 1: Around 4:22, the video says the following. So for those mathematicians, negative numbers didn't exist. You could subtract, that is find the difference between two positive quantities, but you couldn't have a negative answer or negative coefficients. Mathematicians were so averse to negative numbers that there was no single quadratic...
Insights auto threads is broken atm, so I'm manually creating these for new Insight articles. In Dirac’s Principles of Quantum Mechanics published in 1930 he introduced a “convenient notation” he referred to as a “delta function” which he treated as a continuum analog to the discrete Kronecker delta. The Kronecker delta is simply the indexed components of the identity operator in matrix algebra Source: https://www.physicsforums.com/insights/what-exactly-is-diracs-delta-function/ by...
Thread 'Unit Circle Double Angle Derivations'
Here I made a terrible mistake of assuming this to be an equilateral triangle and set 2sinx=1 => x=pi/6. Although this did derive the double angle formulas it also led into a terrible mess trying to find all the combinations of sides. I must have been tired and just assumed 6x=180 and 2sinx=1. By that time, I was so mindset that I nearly scolded a person for even saying 90-x. I wonder if this is a case of biased observation that seeks to dis credit me like Jesus of Nazareth since in reality...
Back
Top