Your only option is likely to iterate over the vectors in nested loops:
v0=0.01;
y =[0.094 0.15];
yv = linspace(min(y), max(y), 25);
c0=[1 2 5 6];
c1=[5 6 8 7];
for k1 = 1:numel(c0)
for k2 = 1:numel(c1)
[y,v(k1, k2, :)] = ode45(@(y,v)((c0(k1)/v^2)-(c1(k2)/v)),yv ,v0);
% disp(v);
end
end
Plotting them is going to be something of a challenge. This could be an option:
figure
hold all
for k = 1:size(v,3)
surf(c0, c1, v(:,:,k))
end
hold off
grid on
xlabel('c0')
ylabel('c1')
zlabel('zv')
view(-45, 20)
