Second exploration of phase transitions: Transitions due to alternate optima

>> x=-1:.00001:1;
>> coste=(abs(x-.5)).^.5+.7*(abs(x+.6)).^.5;
>> plot(x,coste)

>> [calorespecifico,T,costemedio,varianza]=coste2calorespecifico(coste,[0 .5],1);

Apparently, no effect with just a change in the slope:

>> coste=(abs(x-.5)).^1+.7*(abs(x+.6)).^1;
>> plot(x,coste)

>> [calorespecifico,T,costemedio,varianza]=coste2calorespecifico(coste,[0 .5],1);