Changes between Version 7 and Version 8 of u/erica/truncerror


Ignore:
Timestamp:
08/28/13 14:47:42 (11 years ago)
Author:
Erica Kaminski
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • u/erica/truncerror

    v7 v8  
    6767= A check that numerically solving Laplace's Eqn with f(x) = 6x gives the exact solution =
    6868
     69I wanted to check if numerically solving
    6970
     71[[latex($\triangledown ^2 U = 6x$)]]
     72
     73gives the exact solution. This is the same as saying
     74
     75[[latex($A\vec{U} = 6\vec{x}$)]],
     76
     77that is the LTE is zero (i.e. discretization of the equation leads to the exact solution).
     78
     79To check this I wrote out the form for the LTE, which is the same as given above but with f(xi) replaced with 6x(i). Using the exact solution in this formula gives,
     80
     81[[latex($\tau_i = U''_i + \frac{1}{12}h^2 U''''_i + O(h^4) - 6x_i$)]]
     82
     83The ODE is easy to solve by hand, and gives a general solution of,
     84
     85[[latex($U(x) = x^3 + cx + d$)]]
     86
     87As you can see, the truncation error requires 2nd derivatives and ''higher'', and so the LTE is exactly zero as I wanted to verify.
     88
     89
     90
     91
     92
     93