Операция правого деления (/) не определена для векторов (s=1:0.01:9 -- вектор); сравните Yt=1/s; plot(Yt) и Yt=(s)^(-1); plot(Yt) или (равносильно) Yt=ones(s)./s // (покомпонентное деление) plot(Yt) Об этом, к сожалению, не сказано ясно в книге