Legend: fo(t), fr(t), f(t)

M-matrix

-0.2  1/h 0  1/h
3.0  1/h -3.0  1/h
0.2  1/h
3.0  1/h