function [gx, gpx, fpx] = bindingfuncs(X)

h=1;
g1 = 0.0756;
gp1 = 0.0709;
fp1 = 0.2838;
g2=20.*g1;
gx=g2.*(X<0)+g1.*(X).*(X>=0 & X<=h)+ (g1).*(X>h);
gpx= 3*(fp1+gp1).*(X<0)+gp1.*(X).*(X>=0 & X<=h)+(4*gp1).*(X).*(X>h);
fpx=0.*(X<0)+fp1.*(X).*(X>=0 & X<=h)+0.*(X>h);
    
