commit 3027250a9192ce0aeab5eaae04e19a208253b9b2 Author: schrom01 Date: Thu Dec 15 09:20:37 2022 +0100 initial commit diff --git a/HM1_Serie11.pdf b/HM1_Serie11.pdf new file mode 100644 index 0000000..ec2ab39 Binary files /dev/null and b/HM1_Serie11.pdf differ diff --git a/Serie11_Aufg3_Gerüst.py b/Serie11_Aufg3_Gerüst.py new file mode 100644 index 0000000..e4f23e3 --- /dev/null +++ b/Serie11_Aufg3_Gerüst.py @@ -0,0 +1,39 @@ +# -*- coding: utf-8 -*- +""" +Created on Sun Nov 29 14:43:17 2020 + +Höhere Mathematik 1, Serie 11, Aufgabe 3 (Gerüst) + +@author: knaa +""" +import numpy as np +import matplotlib.pyplot as plt + + +detail = 1000 #number of pixels in x and y direction +maxit = 70 #maximum n for iterations +x_min = """???""" #minimim value of x-interval +x_max = """???""" #maximum value of x-interval +y_min = """???""" #minimim vale of y-interval +y_max ="""???""" #minimim vale of y-interval + +a = np.linspace("""???""",detail,dtype=np.float64) #define real axis [x_min,x_max] +b = np.linspace("""???""",detail,dtype=np.float64) #define imaginary axis [y_min,y_max] + +B = np.zeros((detail,detail)) #for color valzues n + +[x,y] = np.meshgrid("""???""") #to create the complex plane with the axes defined by a and b +C = x+y*1j #creating the plane +Z = np.zeros("""???""", np.complex64) #initial conditions (first iteration), Z has same dimension as C +for n in np.arange(1,maxit+1): #start iteration + Z = """???""" #calculating Z + expl = np.where("""???""") #finding exploded values (i.e. with an absolute value > 2) + Z[expl] = 0 #removing from iteration + C[expl] = 0 #removing from plane + B[expl] = """???""" #saving color value n + +plt.figure(1) +B = B/np.max(np.max(B)) #deviding by max value for correct color +plt.imshow(B,extent=[x_min,x_max,y_min,y_max],origin='lower',interpolation='bilinear') #display image + +