diff --git a/Sheet6/Ex6A.m b/Sheet6/Ex6A.m index 2f8a8d0..693c2db 100644 --- a/Sheet6/Ex6A.m +++ b/Sheet6/Ex6A.m @@ -1,6 +1,6 @@ clear; clc; close all; pList = [5, 10, 20, 100]; -nIter = 9; +nIter = 5; figure; tiledlayout(2,2) @@ -22,6 +22,10 @@ for ip = 1:length(pList) K(1,:) = 0; K(1,1) = 1; F(1) = -atan(p); + % GH: Dirichlet boundary conditions at x=+1 + K(end,:) = 0; + K(end,end) = 1; + F(end) = atan(p); u = K \ F; diff --git a/Sheet6/Sheet6.pdf b/Sheet6/Sheet6.pdf new file mode 100644 index 0000000..e5c5e09 Binary files /dev/null and b/Sheet6/Sheet6.pdf differ diff --git a/Sheet6/gh_response.txt b/Sheet6/gh_response.txt new file mode 100644 index 0000000..c728563 --- /dev/null +++ b/Sheet6/gh_response.txt @@ -0,0 +1,16 @@ +h- and r-adaptivity. + +A: +* h-adaptivity gets worse with increasing p (not correct) +* GH: choose nIter=1 ==> all solutions are incorrect +* GH: choose nIter=20 ==> all solutions are correct +* GH: nIter = 5 and Dirchlet bc at x=1.5 ==> aweful for p=20, 100 + Integration of f(x) is not accurate enough. ==> better integration formula needed, see integral(...) + +B: +* h-adaptivity is correct + +C: +* r-adaptivity + looks suspicious, but similar to Mandl +