function diffw=highersys(x,w) % a y" + by '+ c y = f(x) % Only state variable part goes here without ICs A=[-3/2 -5/2; 1 0]; B=[11*exp(-x)/2 ; 0]; diffw=A*w+B; end