MATLAB

[programming]

Добрый день)) Вопрос по программе в matlab.

y=[1.3 1.7 1.8 1.65 1.35]
x=[0.1 0.2 0.3 0.4 0.5] 
   function [s]=lagrange(x,y)
s=0;
n=length(x);
z=linspace(x(1),x(n),100);

   for k=1:n
   p=1;

         for i=[1:k-1,k+1:n]
p=p.*(z-x(i))/(x(k)-x(i));
         end

s=s+y(k).*p;

    end


plot(x,y,'*r')
hold on
grid on
plot(z,s)
end

Ругается на третью строчку, что Function definitions are not permitted in this context.