close all
K=0;

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


F1y=((1+K)./abs(abs(y)-0.07))- (1./0.07);
F1y=F1y./100;

for i=1:max(size(F1y))

    F1y(:,i)=F1y(:,i)+0.5 -0.2*K;

end

for i=1:max(size(F1y))
if y(:,i)>-0.015
    F1y(:,i)=0
end
end
  
for i=1:max(size(F1y))
if F1y(:,i)>2.5
    F1y(:,i)=2.5
end
end

K=1;
F11y=((K+1)./abs(abs(y)-0.07))- (1./0.07);
F11y=F11y./100;

for i=1:max(size(F11y))

    F11y(:,i)=F11y(:,i)+0.5 -0.2*K;

end

for i=1:max(size(F11y))
if y(:,i)>-0.015
    F11y(:,i)=0
end
end
  
for i=1:max(size(F11y))
if F11y(:,i)>2.5
    F11y(:,i)=2.5
end
end
%F2y=alpha.^(yy*zexp)-1;

figure (3)


plot(y,-F1y)
hold on
plot(y,-F11y,'r')
hold off


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
K=0;

F1yy=((K+1)./abs(abs(yy)-0.07))- (1./0.07);
F1yy=F1yy./100;

for i=1:max(size(F1yy))

    F1yy(:,i)=F1yy(:,i)+0.5 - 0.2*K;

end

for i=1:max(size(F1yy))
if yy(:,i)<0.015
    F1yy(:,i)=0
end
end
  
for i=1:max(size(F1yy))
if F1yy(:,i)>2.5
    F1yy(:,i)=2.5
end
end

K=1;
%F2y=alpha.^(yy*zexp)-1;
F11yy=((K+1)./abs(abs(yy)-0.07))- (1./0.07);
F11yy=F11yy./100;

for i=1:max(size(F11yy))

    F11yy(:,i)=F11yy(:,i)+0.5 -0.2*K;

end

for i=1:max(size(F11yy))
if yy(:,i)<0.015
    F11yy(:,i)=0
end
end
  
for i=1:max(size(F11yy))
if F11yy(:,i)>2.5
    F11yy(:,i)=2.5
end
end


K=2;
%F2y=alpha.^(yy*zexp)-1;
F111yy=((K+1)./abs(abs(yy)-0.07))- (1./0.07);
F111yy=F111yy./100;

for i=1:max(size(F111yy))

    F111yy(:,i)=F111yy(:,i)+0.5-0.2*K;

end

for i=1:max(size(F111yy))
if yy(:,i)<0.015
    F111yy(:,i)=0
end
end
  
for i=1:max(size(F111yy))
if F111yy(:,i)>2.5
    F111yy(:,i)=2.5
end
end

K=2.2
%F2y=alpha.^(yy*zexp)-1;
F1111yy=((K+1)./abs(abs(yy)-0.07))- (1./0.07);
F1111yy=F1111yy./100;

for i=1:max(size(F1111yy))

    F1111yy(:,i)=F1111yy(:,i)+0.5-0.2*K;

end

for i=1:max(size(F1111yy))
if yy(:,i)<0.015
    F1111yy(:,i)=0
end
end
  
for i=1:max(size(F1111yy))
if F1111yy(:,i)>2.5
    F1111yy(:,i)=2.5
end
end
figure (4)

plot(yy,F1yy)
hold on
plot(yy,F11yy,'r')
hold on
plot(yy,F111yy,'g')
hold on
plot(yy,F1111yy,'y')
hold off
%%%%%%%%%%%%%%%%%%^Boa%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 
% close all
% K=0;
% 
% x=-0.046:0.0001:0;
% xx=0:0.0001:0.046;
% 
% 
% F1x=((K+1)./abs(abs(x)-0.046))- (1./0.046);
% F1x=F1x./100;
% 
% for i=1:max(size(F1x))
% 
%     F1x(:,i)=F1x(:,i)+0.5-0.25*K;
% 
% end
% 
% for i=1:max(size(F1x))
% if x(:,i)>-0.008
%     F1x(:,i)=0
% end
% end
%   
% for i=1:max(size(F1x))
% if F1x(:,i)>2.5
%     F1x(:,i)=2.5
% end
% end
% 
% K=1;
% F11x=((K+1)./abs(abs(x)-0.046))- (1./0.046);
% F11x=F11x./100;
% 
% for i=1:max(size(F11x))
% 
%     F11x(:,i)=F11x(:,i)+0.5-0.25*K;
% 
% end
% 
% for i=1:max(size(F11x))
% if x(:,i)>-0.008
%     F11x(:,i)=0
% end
% end
%   
% for i=1:max(size(F11x))
% if F11x(:,i)>2.5
%     F11x(:,i)=2.5
% end
% end
% %F2y=alpha.^(yy*zexp)-1;
% 
% figure (5)
% 
% 
% plot(x,-F1x)
% hold on
% plot(x,-F11x,'r')
% hold off
% 
% 
% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 
% K=0;
% F1xx=((K+1)./abs(abs(xx)-0.046))- (1./0.046);
% F1xx=F1xx./100;
% 
% for i=1:max(size(F1xx))
% 
%     F1xx(:,i)=F1xx(:,i)+0.5-0.25*K;
% 
% end
% 
% for i=1:max(size(F1xx))
% if xx(:,i)<0.008
%     F1xx(:,i)=0
% end
% end
%   
% for i=1:max(size(F1xx))
% if F1xx(:,i)>2.5
%     F1xx(:,i)=2.5
% end
% end
% K=1;
% %F2y=alpha.^(yy*zexp)-1;
% F11xx=((K+1)./abs(abs(xx)-0.046))- (1./0.046);
% F11xx=F11xx./100;
% 
% for i=1:max(size(F11xx))
% 
%     F11xx(:,i)=F11xx(:,i)+0.5-0.25*K;
% 
% end
% 
% for i=1:max(size(F11xx))
% if xx(:,i)<0.008
%     F11xx(:,i)=0
% end
% end
%   
% for i=1:max(size(F11xx))
% if F11xx(:,i)>2.5
%     F11xx(:,i)=2.5
% end
% end
% 
% 
% K=2;
% %F2y=alpha.^(yy*zexp)-1;
% F111xx=((K+1)./abs(abs(xx)-0.046))- (1./0.046);
% F111xx=F111xx./100;
% 
% for i=1:max(size(F111xx))
% 
%     F111xx(:,i)=F111xx(:,i)+0.5-0.25*K;
% 
% end
% 
% for i=1:max(size(F111xx))
% if xx(:,i)<0.008
%     F111xx(:,i)=0
% end
% end
%   
% for i=1:max(size(F111xx))
% if F111xx(:,i)>2.5
%     F111xx(:,i)=2.5
% end
% end
% 
% K=3;
% %F2y=alpha.^(yy*zexp)-1;
% F1111xx=((K+1)./abs(abs(xx)-0.046))- (1./0.046);
% F1111xx=F1111xx./100;
% 
% for i=1:max(size(F1111xx))
% 
%     F1111xx(:,i)=F1111xx(:,i)+0.5-0.25*K;
% 
% end
% 
% for i=1:max(size(F1111xx))
% if xx(:,i)<0.008
%     F1111xx(:,i)=0
% end
% end
%   
% for i=1:max(size(F1111xx))
% if F1111xx(:,i)>2.5
%     F1111xx(:,i)=2.5
% end
% end
% figure (6)
% 
% plot(xx,F1xx)
% hold on
% plot(xx,F11xx,'r')
% hold on
% plot(xx,F111xx,'g')
% hold on
% plot(xx,F1111xx,'y')
% hold off
% 
% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%