Fixed small error in alpha. Made sure all returned vectors have length n+1. crankNicolson() now has a supermysterious bug which causes the return value to be a little wrong.