Question: Could not get the solution

Hi everyone, I am trying to dsolve a ode but could not get the answer, could anyone please help me or guide me what's wrong with my worksheet


 

restart

u := -3*beta*(2*p^3*sigma^4*(1/3)+(2*p^3*y^2*(1/3)-4*q*(1/3))*sigma^2-4*p^2*y*(k+1)*sigma*(1/3)+p*(k+1)^2)*(y+sigma)*(y-sigma)/(8*sigma^2)+(-p*sigma^3+(p*y^2-k+1)*sigma-(k+1)*y)/(2*sigma)

-(3/8)*beta*((2/3)*p^3*sigma^4+((2/3)*p^3*y^2-(4/3)*q)*sigma^2-(4/3)*p^2*y*(k+1)*sigma+p*(k+1)^2)*(y+sigma)*(y-sigma)/sigma^2+(1/2)*(-p*sigma^3+(p*y^2-k+1)*sigma-(k+1)*y)/sigma

(1)

ode := diff(theta(y), y, y)+G*(diff(u(y), y))^2+G*beta*(diff(u(y), y))^4

diff(diff(theta(y), y), y)+G*(-(3/8)*(diff(beta(y), y))*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(y(y)+sigma(y))*(y(y)-sigma(y))/sigma(y)^2-(3/8)*beta(y)*(2*p(y)^2*sigma(y)^4*(diff(p(y), y))+(8/3)*p(y)^3*sigma(y)^3*(diff(sigma(y), y))+(2*p(y)^2*y(y)^2*(diff(p(y), y))+(4/3)*p(y)^3*y(y)*(diff(y(y), y))-(4/3)*(diff(q(y), y)))*sigma(y)^2+2*((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)*(diff(sigma(y), y))-(8/3)*p(y)*y(y)*(k(y)+1)*sigma(y)*(diff(p(y), y))-(4/3)*p(y)^2*(diff(y(y), y))*(k(y)+1)*sigma(y)-(4/3)*p(y)^2*y(y)*(diff(k(y), y))*sigma(y)-(4/3)*p(y)^2*y(y)*(k(y)+1)*(diff(sigma(y), y))+(diff(p(y), y))*(k(y)+1)^2+2*p(y)*(k(y)+1)*(diff(k(y), y)))*(y(y)+sigma(y))*(y(y)-sigma(y))/sigma(y)^2-(3/8)*beta(y)*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(diff(y(y), y)+diff(sigma(y), y))*(y(y)-sigma(y))/sigma(y)^2-(3/8)*beta(y)*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(y(y)+sigma(y))*(diff(y(y), y)-(diff(sigma(y), y)))/sigma(y)^2+(3/4)*beta(y)*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(y(y)+sigma(y))*(y(y)-sigma(y))*(diff(sigma(y), y))/sigma(y)^3+(1/2)*(-(diff(p(y), y))*sigma(y)^3-3*p(y)*sigma(y)^2*(diff(sigma(y), y))+((diff(p(y), y))*y(y)^2+2*p(y)*y(y)*(diff(y(y), y))-(diff(k(y), y)))*sigma(y)+(p(y)*y(y)^2-k(y)+1)*(diff(sigma(y), y))-(diff(k(y), y))*y(y)-(k(y)+1)*(diff(y(y), y)))/sigma(y)-(1/2)*(-p(y)*sigma(y)^3+(p(y)*y(y)^2-k(y)+1)*sigma(y)-(k(y)+1)*y(y))*(diff(sigma(y), y))/sigma(y)^2)^2+G*beta*(-(3/8)*(diff(beta(y), y))*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(y(y)+sigma(y))*(y(y)-sigma(y))/sigma(y)^2-(3/8)*beta(y)*(2*p(y)^2*sigma(y)^4*(diff(p(y), y))+(8/3)*p(y)^3*sigma(y)^3*(diff(sigma(y), y))+(2*p(y)^2*y(y)^2*(diff(p(y), y))+(4/3)*p(y)^3*y(y)*(diff(y(y), y))-(4/3)*(diff(q(y), y)))*sigma(y)^2+2*((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)*(diff(sigma(y), y))-(8/3)*p(y)*y(y)*(k(y)+1)*sigma(y)*(diff(p(y), y))-(4/3)*p(y)^2*(diff(y(y), y))*(k(y)+1)*sigma(y)-(4/3)*p(y)^2*y(y)*(diff(k(y), y))*sigma(y)-(4/3)*p(y)^2*y(y)*(k(y)+1)*(diff(sigma(y), y))+(diff(p(y), y))*(k(y)+1)^2+2*p(y)*(k(y)+1)*(diff(k(y), y)))*(y(y)+sigma(y))*(y(y)-sigma(y))/sigma(y)^2-(3/8)*beta(y)*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(diff(y(y), y)+diff(sigma(y), y))*(y(y)-sigma(y))/sigma(y)^2-(3/8)*beta(y)*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(y(y)+sigma(y))*(diff(y(y), y)-(diff(sigma(y), y)))/sigma(y)^2+(3/4)*beta(y)*((2/3)*p(y)^3*sigma(y)^4+((2/3)*p(y)^3*y(y)^2-(4/3)*q(y))*sigma(y)^2-(4/3)*p(y)^2*y(y)*(k(y)+1)*sigma(y)+p(y)*(k(y)+1)^2)*(y(y)+sigma(y))*(y(y)-sigma(y))*(diff(sigma(y), y))/sigma(y)^3+(1/2)*(-(diff(p(y), y))*sigma(y)^3-3*p(y)*sigma(y)^2*(diff(sigma(y), y))+((diff(p(y), y))*y(y)^2+2*p(y)*y(y)*(diff(y(y), y))-(diff(k(y), y)))*sigma(y)+(p(y)*y(y)^2-k(y)+1)*(diff(sigma(y), y))-(diff(k(y), y))*y(y)-(k(y)+1)*(diff(y(y), y)))/sigma(y)-(1/2)*(-p(y)*sigma(y)^3+(p(y)*y(y)^2-k(y)+1)*sigma(y)-(k(y)+1)*y(y))*(diff(sigma(y), y))/sigma(y)^2)^4

(2)

dsolve(ode)

Error, (in dsolve) Required a specification of the indeterminate function

 

bc := theta(-sigma) = 0, theta(sigma) = 1

theta(-sigma) = 0, theta(sigma) = 1

(3)

sol := dsolve({bc, ode})

(4)

``


 

Download help_ode.mw

Please Wait...