第五次数值分析作业

Author: lxb1226

Diff for: simpson.c

@@ -0,0 +1,27 @@
+//
+// Created by heyjude on 18-10-25.
+//
+
+#include<stdio.h>
+#include<math.h>
+
+//float f(float x);
+//
+//int main(void){
+//    float b=2.0,a=0.0;
+//    int m=5;
+//    float h=(b-a)/(2*m);
+//    int i;
+//    float sum=0;
+//    for(i=1;i<=m;i++){
+//        sum+=4*f(a+h*(2*i-1))+2*f(a+2*i*h);
+//
+//    }
+//    sum=h/6*(f(a)+f(b)+sum);
+//    printf("sum = %f\n",sum);
+//    return 0;
+//}
+//
+//float f(float x){
+//    return sin(x)*exp(-x*x);
+//}

Diff for: tixing.c

@@ -0,0 +1,25 @@
+#include<stdio.h>
+#include<math.h>
+
+double f(float x);
+
+int main(void){
+    int M=10;
+    int i;
+    float b=2.0,a=0.0;
+    float h;
+    double sum=0;
+    h=(b-a)/M;
+    printf("h = %f\n",h);
+    for(i=1;i<=M-1;i++){
+        sum+=f(a+i*h);
+    }
+    sum=h/2*(f(a)+f(b)+2*sum);
+    printf("sum = %f",sum);
+    return 0;
+}
+double f(float x){
+
+    return sin(x)*exp(-x*x);
+
+}