Function with matrix input arguments

Hi, I want to evaluate the following equation

F=((sqrt(2)/(81*sqrt(pi)))*(6*R-R^2)*exp(-R/3)*cos(TH))^2;

Both R and TH are matrices. How can I make sure that F can handle matrix inputs?

I later want to plot this function using:

[X,Y] = pol2cart(TH,R)

surf(X,Y,F)

Answer 1

[KSSV]

r = linspace(0,1) ;
th = linspace(0,2*pi) ;
[R,TH] = meshgrid(r,th) ;
F=((sqrt(2)/(81*sqrt(pi)))*(6*R-R^2)*exp(-R/3)*cos(TH))^2;
[X,Y] = pol2cart(TH,R) ;
surf(X,Y,F)

Answer 2

[Walter Roberson]

Vectorize.

F = ((sqrt(2) ./ (81*sqrt(pi))) .* (6*R-R.^2) .* exp(-R/3) .* cos(TH)).^2;