I am writing a code for a system of equations for which the variables vary in time and space. I have written an implicit code to solve a system of equations before but in that case I could write the equations in matrix form and used Newton Raphson's method to solve the equations. This time I am just using for loops and I'm not assembling the matrices. Would it be possible to use Newton Raphson's method and Implicit Euler for this case. I mean I was having problems in terms of calculating the Jacobian (the gradients of the function). Any hints would be appreciated.
2026-04-12 22:48:09.1776034089
Solving a system of equations that vary in time and space implicitly
385 Views Asked by Bumbble Comm https://math.techqa.club/user/bumbble-comm/detail At
1
There are 1 best solutions below
Related Questions in NUMERICAL-METHODS
- The Runge-Kutta method for a system of equations
- How to solve the exponential equation $e^{a+bx}+e^{c+dx}=1$?
- Is the calculated solution, if it exists, unique?
- Modified conjugate gradient method to minimise quadratic functional restricted to positive solutions
- Minimum of the 2-norm
- Is method of exhaustion the same as numerical integration?
- Prove that Newton's Method is invariant under invertible linear transformations
- Initial Value Problem into Euler and Runge-Kutta scheme
- What are the possible ways to write an equation in $x=\phi(x)$ form for Iteration method?
- Numerical solution for a two dimensional third order nonlinear differential equation
Related Questions in PARTIAL-DIFFERENTIAL-EQUATIONS
- PDE Separation of Variables Generality
- Partial Derivative vs Total Derivative: Function depending Implicitly and Explicitly on Variable
- Transition from theory of PDEs to applied analysis and industrial problems and models with PDEs
- Harmonic Functions are Analytic Evan’s Proof
- If $A$ generates the $C_0$-semigroup $\{T_t;t\ge0\}$, then $Au=f \Rightarrow u=-\int_0^\infty T_t f dt$?
- Regular surfaces with boundary and $C^1$ domains
- How might we express a second order PDE as a system of first order PDE's?
- Inhomogeneous biharmonic equation on $\mathbb{R}^d$
- PDE: Determine the region above the $x$-axis for which there is a classical solution.
- Division in differential equations when the dividing function is equal to $0$
Related Questions in MATLAB
- Taking snapshots of an animation in PDE toolbox in Matlab
- Including a time delay term for a differential equation
- Dealing with a large Kronecker product in Matlab
- Apply affine heat equation on images
- How to construct a B-spline from nodal point in Matlab?
- How to solve an algebraic Riccati equation when the Hamiltonian spectrum is too close to the imaginary axis.
- Error calculating diffusion equation solution by fft
- How to simulate a random unitary matrix with the condition that each entry is a complex number with the absolute value 1 in matlab
- Implementation help for Extended Euclidean Algorithm
- Optimization problem in Matlab
Trending Questions
- Induction on the number of equations
- How to convince a math teacher of this simple and obvious fact?
- Find $E[XY|Y+Z=1 ]$
- Refuting the Anti-Cantor Cranks
- What are imaginary numbers?
- Determine the adjoint of $\tilde Q(x)$ for $\tilde Q(x)u:=(Qu)(x)$ where $Q:U→L^2(Ω,ℝ^d$ is a Hilbert-Schmidt operator and $U$ is a Hilbert space
- Why does this innovative method of subtraction from a third grader always work?
- How do we know that the number $1$ is not equal to the number $-1$?
- What are the Implications of having VΩ as a model for a theory?
- Defining a Galois Field based on primitive element versus polynomial?
- Can't find the relationship between two columns of numbers. Please Help
- Is computer science a branch of mathematics?
- Is there a bijection of $\mathbb{R}^n$ with itself such that the forward map is connected but the inverse is not?
- Identification of a quadrilateral as a trapezoid, rectangle, or square
- Generator of inertia group in function field extension
Popular # Hahtags
second-order-logic
numerical-methods
puzzle
logic
probability
number-theory
winding-number
real-analysis
integration
calculus
complex-analysis
sequences-and-series
proof-writing
set-theory
functions
homotopy-theory
elementary-number-theory
ordinary-differential-equations
circles
derivatives
game-theory
definite-integrals
elementary-set-theory
limits
multivariable-calculus
geometry
algebraic-number-theory
proof-verification
partial-derivative
algebra-precalculus
Popular Questions
- What is the integral of 1/x?
- How many squares actually ARE in this picture? Is this a trick question with no right answer?
- Is a matrix multiplied with its transpose something special?
- What is the difference between independent and mutually exclusive events?
- Visually stunning math concepts which are easy to explain
- taylor series of $\ln(1+x)$?
- How to tell if a set of vectors spans a space?
- Calculus question taking derivative to find horizontal tangent line
- How to determine if a function is one-to-one?
- Determine if vectors are linearly independent
- What does it mean to have a determinant equal to zero?
- Is this Batman equation for real?
- How to find perpendicular vector to another vector?
- How to find mean and median from histogram
- How many sides does a circle have?
It is difficult to use the Newton-Raphson method without assembling the matrix. The reason is that Newton-Raphson requires you to solve a matrix equation with the Jacobian, and the standard method for solving matrix equations is Gaussian elimination / LU-factorization, and for that you do need to assemble the matrix. There are methods for solving matrix equations that do not require you to assemble the matrix but these are more advanced.
If you cannot calculate the Jacobian (your question is a bit vague here), then you cannot use Newton-Raphson.
Assuming that you want to solve your equation without using any built-in solver, perhaps you should first try fixed-point iteration to solve the nonlinear equation in the implicit Euler method as this is the easiest method. There is a very short description on Wikipedia at http://en.wikipedia.org/wiki/Backward_Euler_method . However, fixed-point iteration destroys the stability of the implicit Euler method and that is probably the reason why you are using an implicit method in the first place.