http://micro.stanford.edu/mediawiki/index.php?action=history&feed=atom&title=Fit_a0B.m 2026-07-05T13:08:34Z Revision history for this page on the wiki MediaWiki 1.39.7 http://micro.stanford.edu/mediawiki/index.php?title=Fit_a0B.m&diff=3141&oldid=prev 2009-02-06T07:17:19Z

<p>New page: function [a0, B] = fit_a0B( datafile ) %data = load(&#039;platt.B.dat&#039;); data = load( datafile ); a = data(:,1); p = data(:,2); V = a.^3; [P1, S1] = polyfit(a, p, 1); [P2, S2] = polyfit(V, p...</p> <p><b>New page</b></p><div>function [a0, B] = fit_a0B( datafile )<br /> <br /> %data = load(&#039;platt.B.dat&#039;);<br /> data = load( datafile );<br /> <br /> a = data(:,1);<br /> p = data(:,2);<br /> V = a.^3;<br /> <br /> [P1, S1] = polyfit(a, p, 1);<br /> [P2, S2] = polyfit(V, p, 1);<br /> <br /> a0 = -P1(2)/P1(1);<br /> V0 = a0^3;<br /> B = -V0*P2(1) * 0.1; % conversion from kiloBar to GPa<br /> <br /> disp(sprintf(&#039;a0 = %.6f (A)&#039;,a0));<br /> disp(sprintf(&#039;B = %.6f (GPa)&#039;,B));<br /> <br /> aaxis = [0:0.01:1]*(max(a)-min(a))+min(a);<br /> pfit = polyval(P1,aaxis);<br /> vaxis = aaxis.^3;<br /> <br /> figure(2)<br /> plot(a,p,&#039;o&#039;,aaxis,pfit,&#039;-&#039;);</div>

Euniv