function[x2]=p35zz(x0)
b=0.1;
m(1)=0;
E(1)=0;
x2(1)=x0(1);
x2(2)=x0(1);
a(1)=3093;
b(1)=0;
for i=2:13
e(i)=x0(i)-x2(i);
E(i)=b*e(i)+(1-b)*E(i-1);
m(i)=b*abs(e(i))+(1-b)*m(i-1);
a(i)=abs(e(i))/m(i);
s1(1)=x0(1);
s2(1)=x0(1);
s1(i)=a(i)*x0(i)+(1-a(i))*s1(i-1);
s2(i)=a(i)*s1(i)+(1-a(i))*s2(i-1);
A(i)=2*s1(i)-s2(i);
B(i)=a(i)/(1-a(i))*(s1(i)-s2(i));
x2(i+1)=A(i)+B(i);
end