《数值分析》上机习题

m++)Рfor(i=0;i<=3;i++)Рf[m][i]=0;Рf[0][0]=h;Рf[1][0]=h*y[2];Рf[2][0]=h*(1000-1000*y[1]-100*y[2]);Рf[0][1]=h;Рf[1][1]=h*(y[2]+(1/3.0)*f[0][0]);Рf[2][1]=h*(1000-1000*(y[1]+(1/3.0)*f[1][0])-100*(y[2]+(1/3.0)*f[1][0]));Рf[0][2]=h;Рf[1][2]=h*(y[2]+(1/3.0)*f[1][0]+f[1][1]);Рf[2][2]=h*(1000-1000*(y[1]+(1/3.0)*f[2][0]+f[2][1])-100*(y[2]+(1/3.0)*f[2][0]+f[2][1]));Рf[0][3]=h;Рf[1][3]=h*(y[2]+f[1][0]-f[1][1]+f[1][2]);Рf[2][3]=h*(1000-1000*(y[1]+f[2][0]-f[2][1]+f[2][2])-100*(y[2]+f[2][0]-f[2][1]+Рf[2][2]));Рfor(i=0;i<3;i++)Р y[i]=y[i]+(f[i][0]+3*f[i][1]+3*f[i][2]+f[i][3])/8;Рif(n==50||n==90||n==170||n==200)Р{Рfor(i=0;i<3;i++) b[j][i]=y[i];Рj++;Рn++;Р}Рwhile(n<=200);Рprintf("The result is::\n");Рfor(j=0;j<4;j++)Р {for(i=0;i<3;i++)Р printf("y%d(%5.4f)=%.7f,",i+1,a[j],b[j][i]);Р printf("\n");Р}Р}Р(3)运行结果