Search All of the Math Forum:
Views expressed in these public forums are not endorsed by
NCTM or The Math Forum.


naser
Posts:
61
Registered:
3/20/09


Error, (in int) integer too large in context
Posted:
Mar 23, 2010 3:22 AM


Hello everyone, I wrote a program for calculate of natural frequency of a cantilever beam with RayleighRitz method. This Program works with Chebychev Polynomial. But when I run this program, it has an error. Would you please help me how can I solve this problem? Error: ??? Error using ==> sym.maple Error, (in int) integer too large in context Error in ==> sym.int at 39 r = reshape(maple('map','int',f(:),x),size(f)); Error in ==> eulerbeam at 61 y2= subs(int(dW(i)*dW(j)),x,1)subs(int(dW(i)*dW(j)),x,1); %%%%%%%%%%%%%% Main Program %%%%%%%%%%%%%%%% run constant_cantiliver1 n=20; syms x z T T(1)=1; T(2)=x; for i=3:n T(i)=2*x*T(i1)T(i2); end %%%%%%%%%%%%%%%% FixedFree %%%%%%%%%%%%%%%%%%% FFw=x^44*x^3+6*x^2+28*x+17; W=FFw*T; dW=diff(diff((W))); dks=L/2; %%%dx=L/2*dks M11=zeros(n,n); K11=zeros(n,n); for i=1:n for j=1:n %%Mass Matrix M11(i,j)=subs(int((W(i))*W(j)),x,1)subs(int((W(i))*W(j)),x,1); %%Stiffness Matrix K11(i,j)=subs(int(dW(i)*dW(j)),x,1)subs(int(dW(i)*dW(j)),x,1); end end M11=M11*dks; K11=K11*(1/dks)^3; M=[p*A*M11]; K=[E*I*K11]; y=eig(M^(1)*K); y1=sort(abs(sqrt(y)))/(2*pi);



