matlab编的龙格库塔法解微分方程，求改错

@[TOmatlab编的龙格库塔法解微分方程，求改错 数学模型： 主程序数据： R1=0.019; R2=0.022; R3=0.044; R4=0.062;% 一级杆直径，二级杆直径，泵径,油管径 L1=1100; L2=500; %L1碳纤维杆的长度，L2钢杆长度 L=5; %单元节长 h1=1600;h2=1047; %泵挂深度 ，动液面 hc=h1-h2; %泵沉没度 S1=1900; S

>>阅读原文<<