surface colour bar, octave-compatible

This commit is contained in:
andy 2021-03-10 19:59:34 +00:00
parent 393c8605a2
commit 2a38009d4f
2 changed files with 5 additions and 3 deletions

View File

@ -36,7 +36,8 @@ figure(1)
surf(f_vals, y_vals, transpose(real(cond))); surf(f_vals, y_vals, transpose(real(cond)));
h = gca; h = gca;
rotate3d on rotate3d on
grid(); grid;
colorbar;
set(h, 'xscale', 'log') set(h, 'xscale', 'log')
set(h, 'yscale', 'log') set(h, 'yscale', 'log')
title('2D Sheet Real Conductivity'); title('2D Sheet Real Conductivity');
@ -52,7 +53,8 @@ figure(2)
surf(f_vals, y_vals, transpose(imag(cond))); surf(f_vals, y_vals, transpose(imag(cond)));
h = gca; h = gca;
rotate3d on rotate3d on
grid(); grid;
colorbar;
set(h, 'xscale', 'log') set(h, 'xscale', 'log')
set(h, 'yscale', 'log') set(h, 'yscale', 'log')
title('2D Sheet Imaginary Conductivity'); title('2D Sheet Imaginary Conductivity');

View File

@ -3,7 +3,7 @@ function conductivity = sheet_conductivity(omega, fermi_level, temp, scatter_lif
%% CONSTANTS %% CONSTANTS
e = 1.602e-19; % coulombs e = 1.602e-19; % coulombs
kb = physconst('Boltzmann'); % J/K kb = 1.380e-23; % J/K
hbar = 6.626e-34 / (2*pi); % Js hbar = 6.626e-34 / (2*pi); % Js
%% TERM 1 %% TERM 1