function [ rgb ] = jmol_color( Z )
%JMOL_COLOR Get color assigned to an element.
%   rgb = JMOL_COLOR(Z) returns a vector containing the RGB color assigned
%   to the element with atomic number Z. See 
%   http://jmol.sourceforge.net/jscolors/
%
%   rgb = JMOL_COLOR(symbol) returns a vector containing the RGB color 
%   assigned to the element with the given chemical symbol 

    if ischar(Z)
        Z = chemsym2number(Z);
    end

    colors = [...
     [255,255,255];... H
     [217,255,255];... He
     [204,128,255];... Li
     [194,255,0];... Be
     [255,181,181];... B
     [144,144,144];... C
     [48,80,248];... N
     [255,13,13];... O
     [144,224,80];... F
     [179,227,245];... Ne
     [171,92,242];... Na
     [138,255,0];... Mg
     [191,166,166];... Al
     [240,200,160];... Si
     [255,128,0];... P
     [255,255,48];... S
     [31,240,31];... Cl
     [128,209,227];... Ar
     [143,64,212];... K
     [61,255,0];... Ca
     [230,230,230];... Sc
     [191,194,199];... Ti
     [166,166,171];... V
     [138,153,199];... Cr
     [156,122,199];... Mn
     [224,102,51];... Fe
     [240,144,160];... Co
     [80,208,80];... Ni
     [200,128,51];... Cu
     [125,128,176];... Zn
     [194,143,143];... Ga
     [102,143,143];... Ge
     [189,128,227];... As
     [255,161,0];... Se
     [166,41,41];... Br
     [92,184,209];... Kr
     [112,46,176];... Rb
     [0,255,0];... Sr
     [148,255,255];... Y
     [148,224,224];... Zr
     [115,194,201];... Nb
     [84,181,181];... Mo
     [59,158,158];... Tc
     [36,143,143];... Ru
     [10,125,140];... Rh
     [0,105,133];... Pd
     [192,192,192];... Ag
     [255,217,143];... Cd
     [166,117,115];... In
     [102,128,128];... Sn
     [158,99,181];... Sb
     [212,122,0];... Te
     [148,0,148];... I
     [66,158,176];... Xe
     [87,23,143];... Cs
     [0,201,0];... Ba
     [112,212,255];... La
     [255,255,199];... Ce
     [217,255,199];... Pr
     [199,255,199];... Nd
     [163,255,199];... Pm
     [143,255,199];... Sm
     [97,255,199];... Eu
     [69,255,199];... Gd
     [48,255,199];... Tb
     [31,255,199];... Dy
     [0,255,156];... Ho
     [0,230,117];... Er
     [0,212,82];... Tm
     [0,191,56];... Yb
     [0,171,36];... Lu
     [77,194,255];... Hf
     [77,166,255];... Ta
     [33,148,214];... W 
     [38,125,171];... Re
     [38,102,150];... Os
     [23,84,135];... Ir
     [208,208,224];... Pt
     [255,209,35];... Au
     [184,184,208];... Hg
     [166,84,77];... Tl
     [87,89,97];... Pb
     [158,79,181];... Bi
     [171,92,0];... Po
     [117,79,69];... At
     [66,130,150];... Rn
     [66,0,102];... Fr
     [0,125,0];... Ra
     [112,171,250];... Ac
     [0,186,255];... Th
     [0,161,255];... Pa
     [0,143,255];... U
     [0,128,255];... Np
     [0,107,255];... Pu
     [84,92,242];... Am
     [120,92,227];... Cm
     [138,79,227];... Bk
     [161,54,212];... Cf
     [179,31,212];... Es
     [179,31,186];... Fm
     [179,13,166];... Md
     [189,13,135];... No
     [199,0,102];... Lr
     [204,0,89];... Rf
     [209,0,79];... Db
     [217,0,69];... Sg
     [224,0,56];... Bh
     [230,0,46];... Hs
     [235,0,38];... Mt
     ]; 
     
    rgb = colors(Z,:);
            

end

