It is unclear to me what the expected solution of this problem should look like, and how you tell if the solution is correct or not for the fractional power.
See: Error while trying to solve coupled PDEs - #4 by nate
and
Default absolute tolerance and relative tolerance - #4 by nate
for tips and tricks on very non-linear problems.