Создадим M-функцию, определяющую правую часть уравнения
function f=func44(y,x)
f(1)=sin(y*x)/(y+x);
end
Воспользуемся функцией решения диф. уравнений ode45
>> [T,Y]=ode45(@func44,[0,20],[1])
>> plot(T,Y)
| Рис. 7. Решение дифференциального уравнения
| |
Проведем сплайновую интерполяцию и построим точечный график
>> pp=csape(T,Y(:,1))
pp =
form: 'pp'
breaks: [1x53 double]
coefs: [52x4 double]
pieces: 52
order: 4
dim: 1
>> fnplt(pp)
>> hold on
>> plot(T,Y(:,1),'o')
Рис.8. График численного решения и полученной интерполяции
4. Найти численное (на промежутке 
) и символьное решение дифференциального уравнения 
с помощью функций ode45 и dsolve. Проверить символьное решение путем подстановки в исходное дифференциальное уравнение. Построить графики символьного и численного решения. Провести сравнение.
