Решение системы уравнений с изменяющейся переменной в Mathlab

[reptype]

Необходимо решить систему уравнений:

5*da/dt+2*x(3)+x(4)=0
5*da/dt+x(4)=0
2*x(3)+x(4)=-3

a(t) - синусоидальная функция времени.
Пусть a(t)=sin(t)
Нужно найти решение системы в каждый момент времени.

clear, clc
h=5
t=0:h:15;
n=length(t);
a=t.*sin(t);
for i=1: (n-1)
da(i)=(sin(t(i+1))-sin(t(i)))/h;
end
A=[ da*5 2 1; da*5 0 1; 0 2 1];
b=[0; 0; -3];
x = A\b;

Не знаю, где ошибка. Не могу добиться получения массива решений.