Question

大家帮忙做一个编程题目，嗯，老出错，不知是何原因。

题目如下： 5、已知微分方程组
x'=x(1-x^2-y^2)+y(x^2-y^2-b)

y'=y(1-x^2-y^2)-x(x^2-y^2-b)
计算，b=0.5 和 1.1的解。
提示：利用四阶龙格——库塔法求解。
我的程序如下：
#include<stdio.h>
#include<math.h>
double X(double b,double x,double y)
{
double f1;
f1=x*(1-x*x-y*y)+y*(x*x-y*y-b);
return f1;
}
double Y(double x,double y,double b)
{
double f2;
f2=y*(1-x*x-y*y)-x*(x*x-y*y-b);
return f2;
}
void mian()
{
double x[6],y[6],b[6],h=0.1;
b[0]=0.5;
double k11,k21,k12,k22,k13,k23,k14,k24;
int i;
for(i=0;i<6;i++)
{
b[i]=0.5+i*h;
k11=X(b[i],x[i],y[i]);
k21=Y(b[i],x[i],y[i]);
k12=X(b[i]+h/2,x[i]+h/2*k11,y[i]+h/2*k21);
k22=Y(b[i]+h/2,x[i]+h/2*k11,y[i]+h/2*k21);
k13=X(b[i]+h/2,x[i]+h/2*k12,y[i]+h/2*k22);
k23=Y(b[i]+h/2,x[i]+h/2*k12,y[i]+h/2*k22);
k14=X(b[i]+h/2,x[i]+h/2*k13,y[i]+h/2*k23);
k24=Y(b[i]+h/2,x[i]+h/2*k13,y[i]+h/2*k23);
x[i+1]=x[i]+h/6(k11+2*k12+2*k13+k14);
y[i+1]=y[i]+h/6(k21+2*k22+2*k23+k24);
printf("%lf %lf %lf\n",b[i],x[i],y[i]);
}
}
提示错误为：
x[i+1]=x[i]+h/6(k11+2*k12+2*k13+k14);
y[i+1]=y[i]+h/6(k21+2*k22+2*k23+k24);
这两行term does not evaluate to a function。
求大侠帮忙！！！

Answer 1

x[i+1]=x[i]+h/6(k11+2*k12+2*k13+k14);

y[i+1]=y[i]+h/6(k21+2*k22+2*k23+k24);

6后面你是要表示乘法么？  如果是， 应该用 *

加上*以后的结果：

Answer 2

void mian() 应为void main()