Dateien hochladen nach „“

euler_cauchy.c

@@ -0,0 +1,61 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+
+
+void dglrhs(int neq, double w0, double gamma, double *x, double *f) // t muss nicht übergeben werden, nicht explizit zeitabh.
+{ //für unsere DGL mit neq = 2
+    //sicherstellen wie in Vorlesung:(Brauch man das wirklich?
+    if (neq!= 2)
+    {
+        printf ("neq passt nicht.\n");
+        abort(); //Bedeutung?->Abbruch
+    }
+    //Speicherplatz muss zuvor bereitgestellt werden
+    f[0] = x[1];
+    f[1] = -w0*w0*x[0]-2*gamma*x[1];
+}
+
+void euler (int neq, double h, double w0, double gamma, double *x, double *f, void (*dglrhs) (int, double,
+            double, double*, double*))
+{
+    int i;
+    
+    (*dglrhs)(neq,w0,gamma,x,f); // Iformationen über rechte Seite holen
+    
+    //Euler-Verfahren anwenden und x_n+1 berechnen
+    
+    for (i = 0; i < neq; i++)
+    {
+        x[i] += h*f[i];
+    }
+}
+    
+    
+int main()
+{
+    double h,t0,x0,v0,tend;
+    int neq = 2;
+    double exact, diff;
+    double *x, *f;
+    double t;
+    double w0, gamma;
+    
+    printf("Bitte geben Sie h, t0, x0, v0 und tend, sowie omega0 und gamma ein:\n");
+    scanf("%le %le %le %le %le %le %le", &h, &t0, &x0, &v0, &tend, &w0, &gamma);
+    
+    x = malloc(sizeof(double)*neq);
+    f = malloc(sizeof(double)*neq);
+    
+    x[0] = x0;
+    x[1] = v0;
+    
+    for (t = t0; t <= tend; t+=h)
+    {
+        printf("t= %20.5le ,dgl = %20.5le\n", t, x[0]);
+        
+        euler (neq,h,w0,gamma,x,f,&(dglrhs));
+    }
+    free(x); free(f);
+}
+