|
| |
 精品推荐 |
 |
|
| |
|
|
|
|
数值计算程序大放送-线性代数方程组
|
日期：2008年3月31日 作者： 查看:[大字体
中字体 小字体]
|
q=-ai[u]*ai[v];
s=(ar[v]-ai[v])*(ar[u]+ai[u]);
ar[u]=(p-q)/d;
ai[u]=(s-p-q)/d;
}
p=br[k]*ar[v];
q=-bi[k]*ai[v];
s=(ar[v]-ai[v])*(br[k]+bi[k]);
br[k]=(p-q)/d;
bi[k]=(s-p-q)/d;
for (i=k+1;i<=n-1;i++)
{
u=i*n+k;
for (j=k+1;j<=n-1;j++)
{
v=k*n+j;
l=i*n+j;
p=ar[u]*ar[v];
q=ai[u]*ai[v];
s=(ar[u]+ai[u])*(ar[v]+ai[v]);
ar[l]=ar[l]-p+q;
ai[l]=ai[l]-s+p+q;
}
p=ar[u]*br[k];
q=ai[u]*bi[k];
s=(ar[u]+ai[u])*(br[k]+bi[k]);
br[i]=br[i]-p+q;
bi[i]=bi[i]-s+p+q;
}
}
u=(n-1)*n+n-1;
d=ar[u]*ar[u]+ai[u]*ai[u];
if (d+1.0==1.0)
{
free(js);
printf("err**fail\n");
return(0);
}
p=ar[u]*br[n-1];
q=-ai[u]*bi[n-1];
s=(ar[u]-ai[u])*(br[n-1]+bi[n-1]);
br[n-1]=(p-q)/d;
bi[n-1]=(s-p-q)/d;
for (i=n-2;i>=0;i--)
{
for (j=i+1;j<=n-1;j++)
{
u=i*n+j;
p=ar[u]*br[j];
q=ai[u]*bi[j];
s=(ar[u]+ai[u])*(br[j]+bi[j]);
上一篇:IIS 6.0配置HTTP压缩
下一篇:微软怎么设计Exchange 2003的体系结构
|
| 相关文章： |
|
|
|
| 相关软件： |
|
| |
|
