Recent content by kriss2

I have explicite method and I need something similar for implicite metod k=1:T Pold=P; Uold=U; i=2:X-1; j=2:Y-1; P(i,j)=Pold(i,j) ... + ((Pold(i-1,j)-2*Pold(i,j)+Pold(i+1,j))/X^2 ... + (Pold(i,j-1)-2*Pold(i,j)+Pold(i,j+1))/Y^2)*tau ...

Hi a want to code implicite scheme for PDEs u_t =a laplace u + L p_t p_t = laplace p + u pls help I made a classical FEM discretization, but i don't know how to programe it. (u(n+1)ij - un)/t = a ((u(n+1)(i-1j) -2 u(n+1)(ij) +u(n+1)(i+1 j)) /x2 + (u(n+1)(i j-1)-2u(n+1)(ij) +...