Hi!!!

I want to write a code in Matlab for the Backward Euler Method for 2x2 systems, using the fixed point iteration to find the y^{n+1}.

y1^{n+1}=y1^{n}+h*f(t^{n+1},y1^{n+1},y2^{n+1}) (1)

y2^{n+1}=y2^{n}+h*g(t^{n+1},y1^{n+1},y2^{n+1}) (2)

Could you tell how I use the fixed point iteration??

At (1) the fixed point iteration will calculate y1^{n+1}, y2^{n+1}will be calculated at (2) but it is already used in the equation (1) ...