ex2

ex2/task_c.py

@@ -1,8 +1,19 @@
-import numpy as numpy
+import numpy as np
 import matplotlib.pyplot as plt
 
+x = np.linspace(0,1,200)
+p_list = [-100, -10, -1, 0, 1, 10, 100]
 
-x = np.linspace(0,1,100)
-p = [10**i for i in range(-3,4)]
+for p in p_list:
+    if p != 0:
+        u = (np.exp(p*x)-1)/(np.exp(p)-1)
+    else:
+        u = x
 
-print(p)
+    plt.plot(x,u, label=f"p={p}")
+
+plt.xlabel("x")
+plt.ylabel("u(x)")
+plt.title("Solution of -u''+pu'=0 with u(0)=0 u(1)=1")
+plt.legend()
+plt.show()