solveTridiag(sub.clone(), diag.clone(), sup.clone(), a[dim], np - 2);
solveTridiag(sub.clone(), diag.clone(), sup.clone(), a[dim], np - 2);
solveTridiag(sub.clone(), diag.clone(), sup.clone(), a[dim], np - 2);
solveTridiag(sub.clone(), diag.clone(), sup.clone(), a[dim], np - 2);