function y=f_2(x)

%%
%% Function y=(x-1)^3 + 2x
%%
y=2*x;
I=find(x >= 1);
y(I)=y(I)+(x(I) - 1).^3;

