Y___Y的专栏

#include #include ////////////////////////////////////////////////////////////////////////////////////////// //矩阵结构体 struct Matrix {  int m,n;//m为行数,n为列数  double **pm;//指向矩阵二维数组的指针 }; //初始化矩阵mt，并

//高斯积分程序，高斯积分具有计算速度快，精度高，能计算反常积分等优点#include #include #define PI 3.14159265static const double gp3[3]={-0.77459666924148338, 0.0                , 0.77459666924148338};//高斯点static const double gc3[3]={

//追赶法解三对角阵方程组Ax=b。  //n---方程组的阶数。  //a1,a2,a3依次为三对角阵的三对角。  //成功返回1，失败返回0。  int RunEq(int n,double *a1,double *a2,double *a3,double *b,double *x)  {   int i;   double l;   for(i=1;i   {    if(fabs(a2[i

#include #include #include double cs(double f[],double x[],int n) { double s=0.0,t=0.0; int i,j; for(i=0;i  {   t=1.0;   for(j=0;j   for(j=i+1;j   t=f[i]/t;   s=s+t; if(i>n)break;  }  return s; }d

///////////////快速傅立叶变换程序/////////////////////#include "stdafx.h" #include #define pi 3.14159265359 #define M  8////////////////取反函数////////////////////////// int rebit(int num,int p) {  int i,rb=

#include #include #include #define N 3 void Guass(double a[][N+1],double x[]) {  double l[N];  int i,j,k;  l[0]=1.0;  for(i=0;i  {   if(fabs(a[i][i])   {    printf("不能求解！/n");    return;   }   for(j=i

采用方法：先用二分法求一实根（三次方程必有一实根），然后化为二次方程，用解析法求得另外两根（包括实根和虚根）#include #include #include #define eps 1.0e-8 //三次函数 inline double pow3(double a3,double a2,double a1,double a0,double x) {  return ((a3*x+a