помагите!!!!

[xxx-100]

[code][/code]Проблемма в следующем! Есть функция по имени fac :
function fac=fac(x)
global R R0;
z=sqrt(x);
c1=-besselj(0,z*R);
c2=z*bessely(1,z*R0);
c3=bessely(0,z*R);
c4=z*besselj(1,z*R0);
fac=c3*c4-c1*c2;
Решаем записанную задачу следующем образом.
global R R0
R=1;
R0=2;
k=1;
q=fzero('fac',2);
lamb(1)=q;
for k=1:4
for j=1:50
p=fzero('fac',j);
if p>q+1
k=k+1;
lamb(k)=p;
q=p;
end
end
end
в чем ошибка? Почему for k=1:4 не работает?

[Prorab]

В карантин - неинформативный заголовок (уж одно слово можно было без ошибок написать?!), невнятный вопрос, неудачное оформление (хотя бы тег code нужно было добавить)