%βֵ
function[xp,yp]=fractal_curves(x,y)
alpha=0.4
N=length(x);
xp=[];
yp=[];xp=x;yp=y;
k=1;
for j=1:length(x)-1
a(j)=(x(j+1)-x(j))/(x(N)-x(1));
c(j)=((y(j+1)-y(j))-alpha(y(N)-y(1)))/(x(N)-x(1));
e(j)=(x(N)*x(j)-x(1)*x(j+1))/(x(N)-x(1));
f(j)=(x(N)*y(j)-x(1)*y(j+1)-alpha*(x(N)*y(1)-x(1)*y(N)))/(x(N)-x(1));
 for i=1:N
     xp(k)=a(j)*x(i)+e(j);
     yp(k)=c(j)*x(i)+alpha*y(i)+f(j);
     k=k+1;
 end
end