麻烦各位大神帮忙改一下，谢谢大神

clear all
u0 = 1; f = 1; r = 0.17; lamda = 5*10^(-4); z = 1;
rho= linspace(0,.1,20);
phi = linspace(0,2*pi,20);
[phi,rho]=meshgrid(phi,rho);
phip=@(rp)(pi/2)*(1+(-1)^(ceil((rp^2)/(f.*lamda))));
f1=@(rp,theta,rho, phi，phip)(...
u0/(2*I*lamda)*z*rp/(rp^2 + rho^2 + z^2 - 2*rp*rh*...
cos(theta - phi)* e^(I*(2 pi/lamda*sqrt(rp^2 + rho^2 + z^2 - 2 rp*...
rho*cos(theta- phi)) - phip))
);
X= @(rho,phi)(rho .*cos(phi));
Y=@(rho,phi)(rho.*sin(phi));
Q=real(dblquad(@(rp，theta）f1(rp,theta,rho, phi，phip)，0, r ,0, 2*pi ));
[X,Y,Q]=pol2cart(phi,rho,Q);
surf(X,Y,Q);
麻烦各位大神帮忙改一下，谢谢大神