f=[x(2);-Re*x(2)/J-G*x(1)/J+(1.5*10^(-3)*e*v.*v./2/J)*(10^(-3)/(d0-10^(-3)*x(1))/x(1)+log(d0-10^(-3)*x(1))/(x(1)^2)-log(d0)/(x(1^2)))];矩阵输入 有问题 检查下 比如最后一个
怎么样，看不懂吧，求直了需要姐，所以叫求解
