请教各位一个关于流量计算的问题，谢谢各位了

天然河道，有详细断面资料，断面有三条测流垂线，有水位数据，要计算不同水位下的流量
三条垂线把河道断面分成四份，计算公式为q=0.7*V1*A1+A2*(V1+V2)/2+A3*(V3+V2)/2+0.7*V3*A3
A1,A2,A3,A4 是不同水位下四部分的面积
小弟现在需要这样一个程序，不知大家能否提供一点线索
我们论坛有的下么？