[求助]请高手帮忙用MATLAB编写一个常微分方程的程序 [复制链接]

?i\V^3S n$  
先谢谢各位看我的问题哦！ 我现在遇到这样一个方程： 5 wN)N~JE  
z'=(1.49D-B)/(A-1.49C) w}t}Sh  
A=(9000-z)/sqrt((xd-x)^2+(9000-z)^2) d(R8^v/L  
B=(xd-x)/sqrt((xd-x)^2+(9000-z)^2) h4MBw=Tz~  
C=z/sqrt(x^2+z^2) @~N"MsF3  
D=x/sqrt(x^2+z^2) )1R[X!KQ7  
xd=4000x/(x^2+z^2)其中变量只有x，z  而A,B,C,D，xd均是xz变量 a|[f%T<<  
z'=(1.49D-B)/(A-1.49C)最终一定的边界条件去解决这个方程 请问怎么编写啊 谢谢大家的帮忙啊 [;F!\B-  

各位高手给个思路提示哦  我现在脑袋都大了

用龙哥库塔法，直接就能求出来