syms y(t) t Y a b x(t)
[V,Subs] = odeToVectorField(diff(y, 2) == -1*(a-b*cos(2*t))*y);
M = matlabFunction(V,'Vars',{t,Y,a, b});
tspan = linspace(0, pi, 75);
av = 0:10;
bv = 0:5;
nav = numel(av);
nbv = numel(bv);
yc = cell(nav, nbv);
xc = cell(nava, nbv);
for k1 = 1:nav
for k2 = 1:nbv
[t,y] = ode45(@(t,Y)M(t,Y,av(k1),bv(k2)),tspan,[0 1]);
[t,x] = ode45(@(t,Y)M(t,Y,av(k1),bv(k2)),tspan,[1 0]);
yc{k1,k2} = y;
xc{k1,k2} = x;
end
end
figure
hold on
for k1 = nav
for k2 = nbv
hold on
plot(t, xc{k1,k2}(:,1), 'displayname', sprintf('x{%d,%d}', k1, k2));
plot(t, yc{k1,k2}(:,2), 'displayname', sprintf('y{%d,%d}', k1, k2));
hold off
end
end
hold off
grid
legend show