IN MATLAB

Consider the system of differential equations:

{? y=Ay, a?t?b,

y(a) = ya.

In the above, y is a (column) vector function in t as y(t) = (y1(t),y2(t),...,ym(t))tp, A is a constant coefficient matrix of size m m, and ya is the given initial vector. Here the subscript tp means transpose. Use both forward Eulers method and backward Eulers method to solve the system.

y = feuler ivp sys(n,a,b,A,ya). y = beuler ivp sys(n,a,b,A,ya). ya is the initial column vector. y is the column vector as the numerical approximations.