for i=-1:1:1
    for j=-1:1:1

ode45(@(t,x,a,b) te41(t,x,2,1),[-1.5,100],[i;j],opt1);
hold on
    end 
end
