Добавил:
okley
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Matlab / f23
.m clc; %clean common window
%(y(k-1)-2*y(k)+y(k+1))/(n-1);
n=101;
h=pi/(n-1);
x=0:h:pi;
f=[0,h^2*cos(x(2:n-1)),1]';
v=[-2-h^2,1,zeros(1,n-2)];
M=toeplitz(v);
M(1,1)=1;
M(1,2)=0;
M(n,n-1)=0;
M(n,n)=1;
y=M\f;
plot(x,y),grid
%(y(k-1)-2*y(k)+y(k+1))/(n-1);
n=101;
h=pi/(n-1);
x=0:h:pi;
f=[0,h^2*cos(x(2:n-1)),1]';
v=[-2-h^2,1,zeros(1,n-2)];
M=toeplitz(v);
M(1,1)=1;
M(1,2)=0;
M(n,n-1)=0;
M(n,n)=1;
y=M\f;
plot(x,y),grid
Соседние файлы в папке Matlab