equation.h

求解方程根的函数和一些微分方程迭代的函数

#include "math.h"
#include "stdio.h"

//求实系数代数方程全部根的牛顿-下山法
//其中a(n+1),存放多项式方程的实系数(从后面开始放)
//n整形变量，多项式方程的次数
//xr返回n个根的实部，xr返回n个根的虚部
int equat_newton(double *a,double *xr,double *xi,int n);
//调试结果与材料提供的结果一致，但是在matlab下检验
//两种方法得到的结果完全不同
//把两组解代入代数方程，都有有些根符合有些根不符合
//复数根更是不相符合
//可以用这个函数吗？？？？？？


//全区间积分的定步长龙格-库塔法
//t为积分的起始点，h为定步长，y存放初始值
//z返回k个积分点上的n个未知函数值，n微分函数个数，k积分的步数（含积分点）
void equat_rungekutta(double t,double h,double *y,double *z,int n,int k);
//每次调用该函数的时候需要自己写一个微分方程的函数

//积分一步的变步长龙格-库塔法(与前面函数重载)
void equat_rungekutta(double t,double h,double *y,int n,double eps);
//t为积分的起始点，h为定步长，y存放初始值,返回时存放结果
//n微分函数个数，eps积分精度要求。
//每次调用该函数的时候同样需要自己写一个微分方程的函数


//积分一步的变步长基尔方法
//t积分的起始点，h积分步长，eps，积分的精度要求
//y存放起始点初值，返回时存放结果
//q，在第一调用时赋值0，循环使用
void equat_gill(double t,double h,double eps,double *y,double *q,int n);


//求实系数代数方程全部根的qr方法
int equate_qr(double *a,double *xr,double *xi,double eps,int n,int jt);
//a，存放多项式方程的n+1个系数
//xr，返回n个根的实部，xi返回n个根的虚部
//eps，控制qr分解的精度要求，jt，qr分解的最大迭待次数
//n，方程的次数