竖曲线要素的编程可以通过以下步骤进行：
确定已知要素
变坡点里程桩号
变坡点高程
竖曲线半径
变坡点前坡度（上坡为正，下坡为负）
变坡点后坡度（上坡为正，下坡为负）
计算公式
凹凸型判断：当前坡度 - 后坡度为正，则为凸型，反之为凹型。
转坡角（曲折角）：前坡度 – 后坡度。
竖曲线长：半径 * 转坡角。
切线长：竖曲线长 / 2。
外矢距：切线长的平方 / （2 * 半径）。
待求点到变坡点距离：待求点桩号 – 变坡点桩号（取绝对值）。
曲线起终点桩号：
起点：变坡点的桩号 – 切线长。
终点：变坡点的桩号 + 切线长。
任意点切线标高：变坡点的标高 ± 测点与变坡点里程距离 * 该里程对应坡度。
任意点设计标高：
凸型：该桩号在切线上的设计标高 – 修正值。
凹型：该桩号在切线上的设计标高 + 修正值。
编程实现
根据上述公式，可以使用编程语言（如Python、C++、Java等）编写程序来计算和输出竖曲线要素。
程序可以设计成接收用户输入的已知要素，然后计算并输出相应的竖曲线要素。