#include <math.h>
  double acos( double arg );
    函数返回arg的反余弦值，arg的值应该在-1到1之间

  #include <math.h>
  double asin( double arg );
    函数返回arg的反正弦值，arg的值应该在-1到1之间

  #include <math.h>
  double atan( double arg );
    函数返回arg的反正切值

  #include <math.h>
  double atan2( double y, double x );
    函数返回y/x的反正切值，并且它可以通过x，y的符号判断
    （x，y）所表示的象限，其返回的也是对应象限的角度值

  #include <math.h>
  double ceil( double num );
  double floor( double arg );
    ceil函数返回不小于num的最小整数，如num = 6.04, 则返回7.0
    floor函数返回不大于num的最大的数，如num = 6.04, 则返回6.0

  #include <math.h>
  double cos( double arg );
  double sin( double arg );
  double tan( double arg );
    函数分别返回arg的余弦，正弦，正切值，arg都是用弧度表示

  #include <math.h>
  double cosh( double arg );
  double sinh( double arg );
  double tanh( double arg );
    函数分别返回arg的双曲余弦，双曲正弦，双曲正切，arg都是用弧度表示的

  #include <math.h>
  double fmod( double x, double y );
    函数返回x/y的余数

  #include <stdlib.h>
  div_t div( int numerator, int denominator );
  ldiv_t ldiv( long numerator, long denominator );
    函数返回numerator/demominator操作的，返回一个结构体div_t(ldiv_t)
    div_t(ldiv_t)结构体中定义了quot(商),rem(余数)

  #include <math.h>
  double exp( double arg );
    函数返回e(自然底数)的arg次

  #include <math.h>
  double log( double num );
    函数返回num的自然对数值num应为大于0的数

  #include <math.h>
  double log10( double num );
    函数返回num以10为底的对数值，num也应该为大于0的数

  #include <math.h>
  double pow( double base, double exp );
    函数返回以base为底的exp次，不允许的取值范围:
    当base 为 0 且exp 小于或等于 0
    当base 为 负数 且 exp 不为整数

  #include <math.h>
  double sqrt( double num );
    函数返回num的开方值，num应该为不小于0的值

  #include <math.h>
  double frexp( double num, int* exp );
    函数可以获取科学计数法的参数
    函数返回在0.5到1.0之间的值， 传入exp的参数用于返回num的指数
    (以2为底数计算，即：num = mantissa * (2 ^ exp))

  #include <math.h>
  double ldexp( double num, int exp );
    函数返回num*(2^exp)的值，如果结果溢出，返回HUGE_VAL

  #include <math.h>
  double modf( double num, double *i );
    函数分割num，将整数部分填入i的值，小数部分返回

文章出处：http://www.diybl.com/course/3_program/c/c_js/2007925/73590.html

int a[X][Y];
int b[X][Y];
int c[X][Y];

void matrix(int b[][X],int c[][Y]);
main()
{
int i,j,temp;
//clrsc();
printf("Please input int matrix b[%d][%d]\n",X,Y);
for(i=0;i<Y;i++)
for(j=0;j<Y;j++){
scanf("%d",&temp);
b[i][j]=temp;
}
printf("Please input int matrix c[%d][%d]\n",X,Y);
for(i=0;i<X;i++)
for(j=0;j<Y;j++){
scanf("%d",&temp);
c[i][j]=temp;
}
matrix(b,c);
printf("Now print resource matrix b[%d][%d]=",X,Y);
for(i=0;i<X;i++){
printf("\n");
for(j=0;j<Y;j++)
printf("%d ",b[i][j]);
}
printf("\n");
printf("Now print resource matrix c[%d][%d]=",X,Y);
for(i=0;i<X;i++){
printf("\n");
for(j=0;j<Y;j++)
printf("%d ",c[i][j]);
}
printf("\n");
printf("Now printm multiply results matrix a[%d][%d]=B*C:",X,Y);
for(i=0;i<X;i++){
printf("\n");
for(j=0;j<Y;j++)
printf("%d ",a[i][j]);
}
getch();
return 0;
}
/********************************************************************/
void matrix(int b[][X],int c[][Y])
{
int i,j,k,temp;
for(i=0;i<X;i++)
for(j=0;j<Y;j++){
for(k=0;k<Y;k++)
a[i][j]+=b[i][k]*c[k][j];
}
}

Runge-Kutta非线性高阶单步法，p阶R-K法的整体阶段误差为O(h^p)

R-K四阶算法为：
u(i+1)=u(i)+h*(k1+3*k2+3*k3+k4)/8
k1=f(t(i),u(i))
k2=f(t(i+h/3),u(i+h*k1/3))
k3=f(t(i+h/3),u(i+h*k2/3))
k4=f(t(i+h),u(i+h*k3)) */

#include <iostream>
#include <cmath>
using namespace std;

class RK
{
  private:
     double k1,k2,k3,k4;
     double h,b,u,a;
  public:
     void seth(double l=0){h=l;} //设步长
     void setf(double xa=0,double xb=0,double y=0) //设初值和范围（xa,xb）
     {
      b=xb;
      a=xa;
      u=y;
     }
    double f(double t,double u) //函数值，修改它以适应各自需要
    {
  //函数设定
     double f=u-2*t/u; 
     return f;
    }
    /**//*---------------------------*/
    void dork() //R-K 主函数
    {
     for(int count=0;count<(b-a)/h;count++)
     {
      k1=f(a+count*h,u);
      k2=f(a+count*h+h/3,u+h*k1/3);
      k3=f(a+count*h+2*h/3,u-h*k1/3+h*k2);
      k4=f(a+count*h+h,u+h*k1-h*k2+h*k3);
      u=u+h*(k1+3*k2+3*k3+k4)/8;
       cout<<u<<endl;
     }

    }
}; 

void main()
{
RK my;
my.seth(0.1);
my.setf(0,1,1);
my.dork();
}
 
 

该程序对数据直接进行显示，如要画图，可以加入几行，输出数据，然后进行画图。