function [wprime]=f54(x,w)
% wprime must be a col!
wprime=[w(2); (w(2)+8*x^3*w(1)^3)/x];