2
h_i=t_(i+1)-t_i for i=1,...,n-1 where n is a positive integer.
矩阵Q
是n
通过(n-2)
矩阵通过给出i=1,...,n
和j=2,...n-1
项q_(i,j)
:
q_(j-1,j)=1/h_(j-1)
q_(j,j)=-(1/h_(j-1)+1/h_j)
q_(j+1,j)=1/h_j
q_(i,j)=0 for |i-j|>=2
我希望得到一个矩阵Q
。我如何在R中编写这个矩阵的程序?提前多了。