Question: Error, (in thaw) numeric exception: division by zero

Please sir,

I am trying to solve this problem and is giving me the error code "Error, (in thaw) numeric exception: division by zero
". Anyone with useful informations please.

Thanking you in anticipations for your response. 

restart

h := proc (x) options operator, arrow; 1+x*tan(theta)/delta+phi*sin(2*Pi*x) end proc:

A := M*(1+K):

theta[1] := sqrt((-B+sqrt(-4*A*C+B^2))/(2*A)):

xi := H[a]^2+(1+K)/K[p]:

L := xi[1]*cosh(theta[1]*h(x))+beta*xi[1]*theta[1]*sinh(theta[1]*h(x))*sinh(theta[2]*h(x))-xi[2]*cosh(theta[2]*h(x))-beta*xi[2]*theta[2]*sinh(theta[2]*h(x))*sinh(theta[1]*h(x)):

F := proc (x) options operator, arrow; (2*xi[1]/theta[1]-2*xi[2]/theta[2])*sinh(theta[1]*h(x))*sinh(theta[2]*h(x))-2*L*h(x) end proc:

w := S*(sinh(theta[1]*y)*sinh(theta[2]*h(x))-sinh(theta[2]*y)*sinh(theta[1]*h(x)))/(L*xi):

u := S*(xi[1]*cosh(theta[1]*y)*sinh(theta[2]*h(x))-xi[2]*cosh(theta[2]*y)*sinh(theta[1]*h(x))-L)/(L*xi)-1:

psi := S*(xi[1]*sinh(theta[1]*y)*sinh(theta[2]*h(x))/theta[1]-xi[2]*sinh(theta[2]*y)*sinh(theta[1]*h(x))/theta[2]-L*y)/(L*xi)-y:

q := int(u, y = -h(x) .. h(x)):

"(∂P)/(∂x):=((q+2*h(x))*L*xi)/(F(x)):"

Error, recursive assignment

 

"Deltap:=(∫)[0]^(1)(∂P)/(∂x) ⅆx"


 

Download 123.mw

Please Wait...