CTP5.m

function [z, v] = CTP5(x)

    n = numel(x);
    
    z1 = x(1);
   
    g = 1 + 9 * (sum(x(2 : end)) / (n - 1))^0.25;
    
    z2 = g * (1 - z1 / g);
    
    z = [z1 z2]';
    
    theta = -0.2 * pi; a = 0.2; 
    
    b = 10; c = 2; d = 6; e = 1;
    
    v = cos(theta) * (z2 - e) - sin(theta) * z1 - ...
        a * abs(sin(b * pi * (sin(theta) * (z2 - e) + cos(theta) * z1)^c))^d;
    
    v = -v;
    
    v = max(v, 0);
    
end