数值分析上机题

(t,k)/a[k][k]);\r)\rvoidret()\r(\rintt,z;floatsum;\rx[m-l]=(float)a[m-l][m]/a[m-l][m-l];\rfor(t=(m-2);t>-l;t-)\r(\rsum=0;\rfor(z=(t+l);z<m;z++)\rsum+=a[t][z]*x[z];\rx[t]=(float)(a[t][m]-sum)/a[t][t];\r})\rintmain()\r(\rInput();//readdata\rintr,i\rfor(k=0;k<(m-1);k++)\r(\rintselect();//selectmainelement\ri-select();\rvoidexchange(intg);\rexchange(r);//exchange\rvoidanalyze();\ranalyze();//analyze\r)\rvoidret();\rret();//replaceback\rcout9HThesolutionvectorisbelow:H9endl;\rfor(i=0;i<m;i++)\rcout9,'x[,,9i9,,J=u9xliJ9endl;\rreturn0;\r)\r÷øù™šl==NumericalSolutionofPartialDifferentialEquations===\rNumericalResultAccurateValueEiTors\r(0.2,1.0)3.319473.320120.00045\r(0.4,1.0)4.054324.055200.00088\r(0.6,1.0)4.952384.953030.00065\r(0.8,1.0)4.048976.049650.00067\r™§:Crank-NicolsonÉ[WŽúìW›ÐûÔü—ýŒ¡þÿ