%LINK_BA_CU_FU.M

function [link_bwd_fwd, variables_ba1, variables_ba1p] = link_ba_cu_fu

allsyms    

variables_ba1 = [iv_ba1; la_ba1; ptil_ba1; r_ba1; xi4_ba1]; 

variables_ba1p = [iv_ba1p; la_ba1p; ptil_ba1p; r_ba1p; xi4_ba1p];

variables_ba1p_shiftedfwd = [iv_cu; la_cu; ptil_cu; r_cu; xi4_cu];

link_ba1_cu = -variables_ba1p_shiftedfwd  + variables_ba1p;


link_bwd_fwd = [link_ba1_cu];