x=[0.1:0.001:0.8];
n=size(x,2);
y1=zeros(n,1);
y2=zeros(n,1);
for i=1:n
  y2(i)=0;
  for j=1:2:701
    y2(i)=y2(i)+sin(j*(x(i)-0.4))/j;
  end;
  y2(i)=y2(i)*4/(5*pi)+0.2;
end;
for i=1:n
  y1(i)=-0.1*x(i)^5-x(i)^4-8*x(i)^3+12*x(i)^2-3*x(i)+y2(i);
end;
plot(x,y1,'k');
