matlab dy,高手，请问用matlab如何解下面方程：y*Dy=a+b*y;我的计算结果里面含有wrightOmega ，怎样解出一般解?...

满意答案

weiliyao77

2013.07.20

采纳率：49%    等级：12

已帮助：11659人

#include

#include

using namespace std;

#define N 20

double a[N][N];

double x[N+1];

double b[N+1];

int n;//n方程个数，n未知数个数

int set( )

{

cout<

cin>>n;

int i,j;

cout<

for(i = 1;i <= n;i++)

{

for(j = 1;j<=n;j++)

{

cin>>a[i][j];//系数

}

cin>>b[i];//结果

}

return 0;

}

int find(int k)//寻找第k列主元

{

int i,tag = k;

double maxv=0;

for(i = k;i <= n;i++)

{

if(a[i][k] > maxv)

{

maxv = a[i][k];

tag = i;

}

}

return tag;

}

void exchange(int i1,int i2)//将第 i1 i2行互换

{

int j;

for(j = 1;j <= n;j++)

{

swap(a[i1][j],a[i2][j]);

}

swap(b[i1],b[i2]);

}

void div(int k)//让第k个方程的首项系数为1

{

double temp = a[k][k];

for(int j = k; j <= n;j++)

{

a[k][j]/=temp;

}

b[k]/=temp;

}

void disME(int k)

{

int i,j;

for(i =1 ;i<=n;i++)

{

for(j = i;j<= n;j++)

{

if(a[i][j])

{

if(a[i][j]==1)

{if(j==n)

cout<

else

cout<

}

else

{

if(j==n)

cout<

else

cout<