diff --git a/ex2/task_c.py b/ex2/task_c.py new file mode 100644 index 0000000..a12ae74 --- /dev/null +++ b/ex2/task_c.py @@ -0,0 +1,19 @@ +import numpy as np +import matplotlib.pyplot as plt + +x = np.linspace(0,1,200) +p_list = [-100, -10, -1, 0, 1, 10, 100] + +for p in p_list: + if p != 0: + u = (np.exp(p*x)-1)/(np.exp(p)-1) + else: + u = x + + 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() \ No newline at end of file