i=0;
for w=0:0.1:1,
	i=i+1;
	w1(i)=w;
	m=fsolve(inline('tanh(w*m+w)-m','m','w'),0,[],w);
	c=[1/(1-m^2), -w;-w,1/(1-m^2)];
	klad=inv(c);
	chi_lr(i)=klad(1,2);
	m_ex=(-exp(-w)+exp(3*w))/(3*exp(-w)+exp(3*w));
	chi_ex(i)=m_ex-m_ex^2;
end;
subplot(2,2,1)
plot(w1,chi_lr,w1,chi_ex)
legend('\chi_{lr}','\chi_{ex}')
xlabel('w')
ylabel('\chi')

