计算方法与实习上机题答案

实习题11用两种不容的顺序计算，分析误差的变化6483102N（1）顺序计算源代码INCLUDEINCLUDEVOIDMAINDOUBLESUM0INTN1WHILE1SUMSUM1/POWN,2IFN10000PRINTF“SUND30F“,N,SUMIFN10000BREAKNPRINTF“SUMDFN“,N,SUM结果（2）逆序计算源代码INCLUDEINCLUDEVOIDMAINDOUBLESUM0INTN10000WHILE1SUMSUM1/POWN,2IFN10000PRINTF“SUMD30F“,N,SUMIFNINCLUDEVOIDMAININTI0,NFLOATA1024,B1024,D1024PRINTF“PLEASEINPUTN,N“SCANF“D“,PRINTF“NPLEASEINPUTA1TOANN“FORI1I0IDIBIAI1/DI1PRINTF“NFFN“,D0结果3给出一个有效的算法和一个无效的算法计算积分1010,4NDXYN源代码INCLUDEINCLUDEMAINDOUBLEY_01/40LOG5,Y_1DOUBLEY_210/55010/110/2,Y_3INTN1,M10PRINTF“有效算法输出结果N“PRINTF“Y020F“,Y_0WHILE1Y_110/4NY_0/40PRINTF“YD20F“,N,Y_1IFN10BREAKY_0Y_1NIFN30PRINTF“N“PRINTF“N无效算法的输出结果N“PRINTF“Y1020F“,Y_2WHILE1Y_310/N40Y_2PRINTF“YD20F“,M1,Y_3IFMMAININTNDOUBLESN30000SN3000030/2010/3000001/300010/20FORN30000N2NSNN1SNN10/NN1PRINTF“从大到小顺序计算NSN1000FNSN10000FNSN30000FN“,SN1000,SN10000,SN30000SN230/210/201/30/20FORN3NINCLUDETYPEDEFFLOATPFLOATFLOATFF1FLOATXRETURNXXEXPXFLOATFF2FLOATXRETURNXEXPX1FLOATFF3FLOATXRETURNLOGXX2FLOATANSWERFLOATPFLOATINTK2FLOATM1,N1,X2,A,B,CIFPFF3N2PRINTF“X04F,X14F,“,M,NWHILE1IFFABSMNINCLUDEFLOATGEXIANFLOAT,FLOATFLOATFFLOATMAININTI,JFLOATX122FLOATX22,X3SCANF“D“,IFI1PRINTF“F“,X1ELSEIFI2PRINTF“F“,X2ELSEFORJ3JINCLUDEVOIDCOLPIVOTFLOAT,INT,FLOATVOIDCOLPIVOTFLOATC,INTN,FLOATXINTI,J,T,KFLOATPFORI0IFABSCKN1IKJIFKIFORJIJ0IFORJN1JI1JCIN1NXJCIN1JXICIN1N/CIN1IVOIDMAININTIFLOATX4FLOATC451,1,0,3,4,2,1,1,1,1,3,1,1,3,3,1,2,3,1,4COLPIVOTC0,4,XFORI0IVOIDMAINFLOATX4INTIFLOATA4548,24,0,12,4,24,24,12,12,4,0,6,20,2,2,6,6,2,16,2VOIDDIRECTLUFLOAT,INT,FLOATDIRECTLUA0,4,XFORI0I0IFORRN1RI1RUIN1NUIN1RXRXIUIN1N/UIN1I实习题41、源代码INCLUDEFLOATLAGRANGEFLOATX,FLOATY,FLOATXX,INTN/N为（N1）次插值；INTI,JFLOATA,YY0ANEWFLOATNFORI0IFLOATLAGRANGEFLOATX,FLOATY,FLOATXX,INTN/N为（N1）次插值；INTI,JFLOATA,YY0ANEWFLOATNFORI0IDEFINEN3VOIDDIFFERENCEFLOATY,FLOATF44,INTNINTK,IF00Y0F10Y1F20Y2F30Y3FORK1KINCLUDEVOIDMAININTI,NFLOATA2FLOATX151,15,2,25,3,35,4,45,5,55,6,65,7,75,8,Z15FLOATY15334,7950,12265,15905,18915,21415,23865,25250,26755,28050,29665,30140,31040,31815,32515FORN0NFABSCKN1IKJIFKIFORJIJ0IFORJN1JI1JCIN1NXJCIN1JXICIN1N/CIN1IFLOATPOWERINTI,FLOATVFLOATA1WHILEIAVRETURNA结果实习题61、源代码（1）INCLUDEINCLUDEFLOATCOTESFLOATFFLOAT,FLOATA,FLOATB,INTNINTKFLOATC,C10,C2,C3,C4FLOATHBA/NC2FAH/4C3FAH/2C4FA3H/4FORK1KINCLUDEFLOATCOTESFLOATFFLOAT,FLOATA,FLOATB,INTNINTKFLOATC,C10,C2,C3,C4FLOATHBA/NC2FAH/4C3FAH/2C4FA3H/4FORK1KINCLUDEDOUBLEADAMSDOUBLEFDOUBLEX,DOUBLEY,DOUBLEX0,DOUBLEY0,DOUBLEH,INTNFORINTN0NINCLUDEUSINGNAMESPACESTD/F为函数的入口地址，X0、Y0为初值，XN为所求点，STEP为计算次数DOUBLERUNGE_KUTADOUBLEFDOUBLEX,DOUBLEY,DOUBLEX0,DOUBLEY0,DOUBLEXN,LONGSTEPDOUBLEK1,K2,K3,K4,RESULTDOUBLEHXNX0/STEPIFSTEPVOIDRUNGE_KUTTAFLOATFFLOATX,FLOATY,FLOATA,FLOATB,FLOATY0,INTN,FLOATYYFLOATXA,YY0,K1,K2,K3,K4FLOATHBA/NINTIFORI1I3IK1FX,YK2FXH/2,YHK1/2K3FXH/2,YHK2/2K4FXH/2,YHK3YYHK12K22K22K3K4/6XAIHYYI1YVOIDADAMSFLOATA,FLOATB,INTN,FLOATFFLOATX,FLOATY,FLOATY0INTIFLOATY1,Y2,Y,YP,YC,YY3,H,XPRINTF“X0FTY0FN“,A,Y0RUNGE_KUTTAF,A,B,Y0,N,YYY1YY0Y2YY1YYY2HBA/NFORI1I3IPRINTF“XDFTYDFN“,I,AIH,I,YYI1FORI3INIXAIHYPYH55FX,Y59FXH,Y237FX2H,Y19FX3H,Y0/24YCYH9FXH,YP19FX,Y5FXH,Y2FX2H,Y1/24PRINTF“XDFTYDFN“,I1,XH,I1,YCY0Y1Y1Y2Y2YYYCFLOATFFLOATX,FLOATYRETURN01XXXYYVOIDRUNGE_KUTTAFLOATFFLOATX,F