程式筆記
標籤
C
(1)
maple
(2)
Matlab
(11)
paper
(1)
python
(14)
2011年9月28日 星期三
Center difference matric A
function A=GEA(m)
h=1/m+1;
X=linspace(0,1,m);
x=sin(X(2:m));
A=zeros(m,m);
for i=1:m
A(i,i)=-2;
if i-1 ~= 0
A(i,i-1)=1;
A(i-1,i)=1;
end
end
A=(1/h^2)*A;
沒有留言:
張貼留言
較新的文章
較舊的文章
首頁
訂閱:
張貼留言 (Atom)
沒有留言:
張貼留言