Помогите написать код в Matlab-e или подскажите с чего начать. Задание:
Цитата:
1. Напишите программу в MatLab для поиска приближенного решения краевой задачи для одномерного уравнения гиперболического типа:


Стройте график решения для каждого момента времени

.
2. Постройте график функции

:

Для численного вычисления интеграла используйте метод прямоугольников.
3. Постройте в координатах

трехмерную поверхность

.
Примечание. Графики необходимо выводить в одно графическое окно, но в разные его части. Не допускается строить различные графики в одной плоскости. Графики должны быть оформлены должным образом: оси координат должны быть подписаны, должна быть нанесена сетка, у каждого графики должен быть заголовок.
Нашел пример реализации явной разностной схемы для уравнениях
Код:
function result=yavn(M)
pause;
N = 3*M^2
h=1/M;
t=1/N;
U=zeros(M+1,N+1);
for X=1:M+1,
U(X,1)=sin(pi*(X-1)*h);
end
for J=1:N,
y=(J-1)*t;
U(1,J+1)=2*t/h^2*(U(2,J)-U(1,J)*(1+h/2-h^2/(2*t)));
U(M,J+1)=2*t/h^2*(U(M+1,J)+U(M,J)*(1+h/2-h^2/(2*t)))+t*8*h^2*y^4;
for X=2:M,
x=(X-1)*h;
U(X,J+1)=U(X,J)+t/h^2*(U(X-1,J)-2*U(X,J)+U(X+1,J))+t*(8*x^2*y^4);
end
end
mesh(0:(1/N):1,0:(1/M):1,U),grid;
result=U;