- 积分
- 0
- 注册时间
- 2006-12-20
- 仿真币
-
- 最后登录
- 1970-1-1
|
具体程序如下 谢谢!
- #include "stdafx.h"
- //Add C++ Math Library to project
- #pragma comment(lib, "libmx.lib")
- #pragma comment(lib, "libmatlb.lib")
- #pragma comment(lib, "libmat.lib")
- #pragma comment(lib, "libmmfile.lib")
- #pragma comment(lib, "libmatpm.lib")
- mwArray func1(mwArray ,mwArray );
- DECLARE_FEVAL_TABLE
- FEVAL_ENTRY(func1)
- END_FEVAL_TABLE
- int main(int argc, char* argv[])
- {
- double y00[] = {1,0};
- mwArray TT,YY;
- mwArray tspan;
- mwArray y0(1,2,y00);
- tspan = colon(0,0.2,1);
- TT = ode45(&YY,"func1",tspan.DIN,y0.DIN);
- cout<<TT<<endl<<YY<<endl;
- getchar();
- return (EXIT_SUCCESS);
- }
- mwArray func1(mwArray t,mwArray y)
- {
- // One argument test function.
- //return (y(1,1)-y(1,2),y(1,2)+2*y(1,1));
- mwArray dxtemp;
- dxtemp(1,1)= -y(1,1);//-y(1,2);
- dxtemp(1,2)= -y(1,2);//+2*y(1,1);
- return dxtemp;
- }
复制代码 |
|