Colormap of the difference between a function Q(x,y) and a constant C in the (x,y)-plane

Here we have

So, I want the colormap for the difference C-Q(x,y) in the x-y plane and map out the regions where Q(x,y)=C, C<Q(x,y), and C>Q(x,y)

[KSSV]

What have you tried?

Answer 1

[KSSV]

I am helping you for a major part. Rest you need to follow up .

m = 100 ; n = 100 ;
x = linspace(0,0.01,m) ;
y = linspace(0,0.01,n) ;
[X,Y] = meshgrid(x,y) ;
L = 0.81*10^-3 ;
C = 16.3 ;
N = 10000 ;

Q = zeros(m,n) ;
for i = 1:N
    vn = ((2*i-1)/2*L)^2*pi^2 ;
    A = -vn*X./(1+vn*Y) ;
    Q = Q+1/vn*(1-exp(A))  ;
end
Q = 2./X.*Q ;

pcolor(X,Y,Q) ;
shading interp ;
colorbar

Check the code thoroughly, before using.

[Micheal oguntola]

Thank you @KSSV I was able to modify the code you provided to suit my need.

Answer 2

[Micheal oguntola]

@KSSV, i have been able to write a code that plot Q(x,y) in 3D