s[e_,v_]:=3/2 Log[e]+ Log[v];

Plot3D[s[e,v],{e,0,2},{v,0,2},ColorFunction->"DarkRainbow",AxesLabel->Automatic]

Export["entropy_idealgas.jpg",%13]

s[e_,h_]:=-(1+e/h)/2 Log[(1+e/h)/2] -(1-e/h)/2 Log[(1-e/h)/2]

Plot3D[s[e,h],{e,0,2},{h,-2,2},ColorFunction->"DarkRainbow",AxesLabel->Automatic]

Export["entropy_paramagnet.jpg",%17]

s[m_]:=-(1+m)/2 Log[(1+m)/2] -(1-m)/2 Log[(1-m)/2]

Plot[s[m],{m,-1,1},AxesLabel->Automatic]

Export["entropym_paramagnet.jpg",%22]

GraphicsColumn[{
Plot[Exp[1s[m]],{m,-1,1},PlotRange->All],
Plot[Exp[10s[m]],{m,-1,1},PlotRange->All],
Plot[Exp[50s[m]],{m,-1,1},PlotRange->All],
Plot[Exp[100s[m]],{m,-1,1},PlotRange->All],
Plot[Exp[200s[m]],{m,-1,1},PlotRange->All]
},ImageSize->Full]

Export["entropym_paramagnetN.jpg",%15]

Show[%7,ImageSize->Full]