追赶法解边值条件微分方程(C++)

一、作业

二、程序

#include

using namespace std;

int main(){

double h(0.1),b(-h*h-2);

double a[10],c[10],d[10],x[11],y[11],u[10],v[10];

int i;

for(i=0;i<11;i++)

x[i]=i*h;

y[0]=0;y[10]=1;

a[1]=0;c[1]=1;d[1]=h*h*h;

for(i=2;i<9;i++){

a[i]=1;

c[i]=1;

d[i]=i*h*h*h;

}

a[9]=1;c[9]=0;d[9]=x[9]*h*h-1;//为各参数赋值

u[1]=c[1]/b;v[1]=d[1]/b;

for(i=2;i<9;i++){

u[i]=c[i]/(b-a[i]*u[i-1]);

v[i]=(d[i]-a[i]*v[i-1])/(b-a[i]*u[i-1]);

}//追

y[9]=(d[9]-a[9]*v[8])/(b-a[9]*u[8]);

for(i=8;i>0;i--)

y[i]=v[i]-u[i]*y[i+1];//赶

for(i=0;i<11;i++)

cout<

return 0;

}

三、结果