Need help looping an ode45 function

clear all
xo=[0.05; 0];
ts=[0 3];
[t,x]=ode45(@myfun,ts,xo);
figure(1)
plot(t,x(:,1))

function v = myfun(t,x);
zeta=0.9;
m=5; k=1000;
c=2.*zeta.*sqrt(m*k); 
v=[x(2); x(1).*-k/m+x(2).*-c./m];
end

This is what I have so far. How would I be able to have multiple zeta values (below) loop through the ode45?

zeta=[0, 0.1, 0.25, 0.5, 0.75, 0.9, 1];

[Steven Lord]

Please don't start a new question when you've already received some assistance in your original question.

[Erik Sharrer]

sorry