alpha=5;
k=11;

y=-0.07:0.0001:0;
yy=0:0.0001:0.07;

F1y=-alpha.^(-y*k)+1 ;
F2y=alpha.^(yy*k)-1;
F3y=-(alpha+1).^(-y*k)+1 ;
F4y=(alpha+1).^(yy*k)-1;
F5y=-(alpha+2).^(-y*k)+1 ;
F6y=(alpha+2).^(yy*k)-1;
F7y=-(alpha+3).^(-y*k)+1 ;
F8y=(alpha+3).^(yy*k)-1;
F9y=-(alpha+4).^(-y*k)+1 ;
F10y=(alpha+4).^(yy*k)-1;


for i=1:size(y,2)
if (y(:,i)>-0.015)
    F1y(:,i)=0;
    F3y(:,i)=0;
    F5y(:,i)=0;
    F7y(:,i)=0;
    F9y(:,i)=0;
end

if F1y(:,i)<-2.5
    F1y(:,i)=-2.5
end

if F3y(:,i)<-2.5
    F3y(:,i)=-2.5
end
if F5y(:,i)<-2.5
    F5y(:,i)=-2.5
end

if F7y(:,i)<-2.5
    F7y(:,i)=-2.5
end

if F9y(:,i)<-2.5
    F9y(:,i)=-2.5
end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:size(y,2)
if (yy(:,i)<0.015)
    F2y(:,i)=0;
    F4y(:,i)=0;
    F6y(:,i)=0;
    F8y(:,i)=0;
    F10y(:,i)=0;
end

if F2y(:,i)>2.5
    F2y(:,i)=2.5
end

if F4y(:,i)>2.5
    F4y(:,i)=2.5
end
if F6y(:,i)>2.5
    F6y(:,i)=2.5
end

if F8y(:,i)>2.5
    F8y(:,i)=2.5
end

if F10y(:,i)>2.5
    F10y(:,i)=2.5
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure (1)

%subplot(121)

plot(y,F1y,'LineWidth',2)
hold on
plot(y,F3y,'r','LineWidth',2)
hold on
plot(y,F5y,'g','LineWidth',2)
hold on
plot(y,F7y,'y','LineWidth',2)
hold on
plot(y,F9y,'k','LineWidth',2)
%axis([-0.08 0 -3 0])
legend('Alpha=5','Alpha=6','Alpha=7','Alpha=8','Alpha=9')
hold on
%xlabel('COPy','FontSize',14)
%ylabel('Fz','FontSize',14)
%title('M=-0.015 L=-2.5 K=11','FontSize',14)
hold on 
plot(yy,F2y,'LineWidth',2)
hold on
plot(yy,F4y,'r','LineWidth',2)
hold on
plot(yy,F6y,'g','LineWidth',2)
hold on
plot(yy,F8y,'y','LineWidth',2)
hold on
plot(yy,F10y,'k','LineWidth',2)
axis([-0.075 0.075 -2.6 2.6])
%legend('Alpha=5','Alpha=6','Alpha=7','Alpha=8','Alpha=9')
title('|M|=0.015 |L|=2.5 K=11','FontSize',14)
xlabel('COPy','FontSize',14)
ylabel('Fz','FontSize',14)
hold off
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
alpha=5;
k=17;

x=-0.046:0.0001:0;
xx=0:0.0001:0.046;

F1x=-alpha.^(-x*k)+1 ;
F2x=alpha.^(xx*k)-1;
F3x=-(alpha+1).^(-x*k)+1 ;
F4x=(alpha+1).^(xx*k)-1;
F5x=-(alpha+2).^(-x*k)+1 ;
F6x=(alpha+2).^(xx*k)-1;
F7x=-(alpha+3).^(-x*k)+1 ;
F8x=(alpha+3).^(xx*k)-1;
F9x=-(alpha+4).^(-x*k)+1 ;
F10x=(alpha+4).^(xx*k)-1;

for i=1:size(xx,2)
if (xx(:,i)<0.008)
    F2x(:,i)=0;
    F4x(:,i)=0;
    F6x(:,i)=0;
    F8x(:,i)=0;
    F10x(:,i)=0;
end

if F2x(:,i)>2.5
    F2x(:,i)=2.5
end

if F4x(:,i)>2.5
    F4x(:,i)=2.5
end
if F6x(:,i)>2.5
    F6x(:,i)=2.5
end

if F8x(:,i)>2.5
    F8x(:,i)=2.5
end

if F10x(:,i)>2.5
    F10x(:,i)=2.5
end
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:size(xx,2)
if (x(:,i)>-0.008)
    F1x(:,i)=0;
    F3x(:,i)=0;
    F5x(:,i)=0;
    F7x(:,i)=0;
    F9x(:,i)=0;
end

if F1x(:,i)<-2.5
    F1x(:,i)=-2.5
end

if F3x(:,i)<-2.5
    F3x(:,i)=-2.5
end
if F5x(:,i)<-2.5
    F5x(:,i)=-2.5
end

if F7x(:,i)<-2.5
    F7x(:,i)=-2.5
end

if F9x(:,i)<-2.5
    F9x(:,i)=-2.5
end
end

figure (2)

%subplot(121)

plot(x,F1x,'LineWidth',2)
hold on
plot(x,F3x,'r','LineWidth',2)
hold on
plot(x,F5x,'g','LineWidth',2)
hold on
plot(x,F7x,'y','LineWidth',2)
hold on
plot(x,F9x,'k','LineWidth',2)
axis([-0.05 0.05 -2.6 2.6])
legend('Alpha=5','Alpha=6','Alpha=7','Alpha=8','Alpha=9')
hold on
xlabel('COPx','FontSize',14)
ylabel('Fx','FontSize',14)
title('|M|=0.008 |L|=2.5 K=17','FontSize',14)
%subplot(122)
plot(xx,F2x,'LineWidth',2)
hold on
plot(xx,F4x,'r','LineWidth',2)
hold on
plot(xx,F6x,'g','LineWidth',2)
hold on
plot(xx,F8x,'y','LineWidth',2)
hold on
plot(xx,F10x,'k','LineWidth',2)
%axis([0 0.06 0 3])
%legend('Alpha=5','Alpha=6','Alpha=7','Alpha=8','Alpha=9')
hold off
%title('M=0.008 L=2.5 K=17','FontSize',14)
%xlabel('COPx','FontSize',14)
%ylabel('Fx','FontSize',14)

