水仙花数(matlab编程求解)