#include "stdafx.h"
#include "orsciJWVCL.h"
#include "orsciVM.h"
using namespace orsci;
using namespace orsci::vmt;
int main()
{
cout << " orsci: 方程求解,常用solve(...)函数,也可用求逆方法 ... http://www.orsci.cn" << endl;
cout << endl;
mdouble A;
A = "6, 2, 3; 2, 2, 4; 5, 2, 3"; //字符串表达式构造向量,也可构造索引。
cout << A << endl;
cout << "A的秩=" << endl;
cout << vmt::rank(A) << endl;
cout << "A的行列是值=" << endl;
cout << vmt::det(A) << endl;
coldouble B = "1, 2, 3";
cout << "B = " << endl;
cout << B << endl;
mdouble x;
vmt::solve(A, B, x);
cout << "第1种解法:调用solve(...)求解。方程求解 x = " << endl;
cout << x << endl;
cout << "A*x = " << endl;
cout << A * x << endl;
cout << "inv(A) = " << endl;
cout << inv(A) << endl;
cout << "inv(A) * A = " << endl;
cout << inv(A) * A << endl;
cout << "第2种解法:inv(A) * B = " << endl;
cout << inv(A) * B << endl;
cout << endl;
cout << "press any key to stop..." << endl;
char pp;
cin >> pp;
return 0;
} |