close all;clear all;clc;
tao=-10:0.05:10;
miu=tao;
T=1;
[TAO,MIU]=meshgrid(tao,miu);
R=sqrt(tao.^2+miu.^2)+eps;%tao+miu;%.*(T-abs(tao)))+eps;

%KAI_pcsf=abs(R/pi./miu);
surf(TAO,MIU,R);