(Subscript[I, n])=(1/(Δ t))(((Subscript[ϕ, n])((Subscript[t, 0]),Δ t))^T)p