请教一个算法问题

疯长的野草

大家帮看下这个问题怎么处理：
我要编程序解微分代数方程：
　　　　　　　　　　　dx/dt=f(t,x,y)
                      g(t,x,y)=0
我的问题是: 怎样在程序里识别哪些是代数方程g(t,x,y)，哪些是微分方程dx/dt=f(t,x,y).因为这两种方程的解的表达式不一样。
现在的问题是我只想让用户指定总的方程数N，而让程序进行分开识别微分和代数方程。这怎么现实啊？
我现在这里谢过大家乐。

MMMIX

原帖由 疯长的野草 于 2007-5-23 20:14 发表
大家帮看下这个问题怎么处理：
我要编程序解微分代数方程：
　　　　　　　　　　　dx/dt=f(t,x,y)
                      g(t,x,y)=0
我的问题是: 怎样在程序里识别哪些是代数方程g(t,x,y)，哪些是微分方程d ...

建议还是找些专门/专业点的资料自己先看看。

freecutelei

没看明白。。。

doctorjxd

建议楼主找一本《数值分析》的教材来看一看。

对常微分方程组用程序求数值解，可以考虑尤拉算法，或龙格-库塔算法。

疯长的野草

回七楼的，首先谢谢你提供的建议。我举个例子：dy/dt=5x+3y
                                 0=x+y
如果按你的说法，给定一个x0,然后计算出由微分方程计算出y，然后再由y带入代数方程计算x,那么此时的x只满足代数方程的解，而不满足微分方程的解。所以用纯解微分方程的算法是不行的。

cugb_cat

那就要看楼主怎么写这个方程组了，如果是写在文本文件中，可以根据微分方程和代数方程的各自的特点来进行区分

疯长的野草

回９楼的，我是在vc6.0里面写着程序。