Solving Diffusion Equation with Modified Euler Algorithm

[squenshl]

Consider the diffusion equation:
u_t = u_xx, 0 < x < 1, t > 0,
u(0,t) = p(t), u(1,t) = q(t), t > 0,
u(x,0) = f(x), 0 < x < 1,
for t < t_final

Suppose that the boundary condition at x = 0 is now replaced by
$$\partial$$u/$$\partial$$x(0,t) = 0

Using the forward difference formula
$$\partial$$u/$$\partial$$x(x,t) = ((-u(x+2$$\Delta$$x,t) + 4u(x+$$\Delta$$x,t) - 3u(x,t))/2$$\Delta$$x) + O($$\Delta$$x)²

show that the Euler algorithm needs to be modified to include the formula for U_0,k+1:
U_j,k+1 = rU_j-1,k + (1-2r)U_j,k + rU_j+1,k, j = 1,2,...,N-1,
U_0,k+1 = (4U_1,k+1 - U_2,k+1)/3.

I don't have a clue where to start. Please help.

 

[HallsofIvy]

If you are, in fact, taking a course in the numerical solution of differential equations and have no "clue where to start" with all of the problems you have posted, then you have far worse problems that we can help you with! Talk to your teacher about this.