"??? Undefined function or method 'int' for input arguments of type 'double'."

Hi all,

I have been trying to solve this problem. To set up a e array[200 by 200] to create some mesh beam diagram. ' However,I keep on receiving this error message:

"??? Undefined function or method 'int' for input arguments of type 'double'." on the integration (int) part no matter how many times I tried.

i tried using quadgk function as searched on the net, but the same message appears.

I wish to obtain numeric integration answer, so I didn't use syms x.

Really appreciate your help !

e=zeros(200); %store the calculated value in e array 

double x;

for i=1:1:200       

for j=1:1:200

          x=i-100;

          y=j-100;

          rho=sqrt(x^2+y^2); %To generate rho values to sub into e equation

e(i,j)=double(int(cos(x)^0.5*sin(2*x)*((exp(1)^(-(9/4)*(sin(x)/sin(1.125))^2))*(besselj(1,3*(sin(x)/sin(1.125)))))*(besselj(1,2*3.142*rho*sin(x))),0,1.125));

       end
  end

% %

Answer 1

[Walter Roberson]

You appear to be attempting to do numeric integration. int() is only for symbolic integration. There are a number of numeric integration routines such as quadgk()

Side note:

double x;

would be the same thing as

double('x');

which would convert the letter 'x' to double precision and then throw away the result.

MATLAB does not use declarations of type. Dataypes are determined by assignment.

[Shing]

Hi, yes. I'm trying to do numeric integration and tried using quadgk() previously. However, when I did so, I got this error: ??? Undefined function or method 'quadqk' for input arguments of type 'function_handle'. I'm not sure where it went wrong :( Please help ! Thank you in advance!

e=zeros(200);

for i=1:1:200 %range will be -100 to 100 for x column for j=1:1:200 x=i-100; y=j-100; rho=sqrt(x^2+y^2); e(i,j)= quadqk(@myfun,0,1.125);

-------------functionfile------------------------------------------------------------------------------------------------------------ function z = myfun(x)

z= (cos(x)^0.5.*sin(2*(x)).*((exp(1)^(-(9/4)*(sin(x)./sin(1.125))^2))*(besselj(1,3*(sin(x)./sin(1.125))))).*(besselj(1,2.*3.142.*rho.*sin(x))));

[Walter Roberson]

You used quadQk rather than quadGk

[Shing]

oh terrible mistake, thanks for the help :)