//Example number 3.6, Page number 3.34


clc;clear;close

// Variable declaration
H0=64*10**3;    // initial field(ampere/m)
T=5;    // temperature(K)
Tc=7.26;   // transition temperature(K)

// Calculation
H=H0*(1-(T/Tc)**2);     // critical field(ampere/m)

// Result
printf("critical field is : %0.3e ampere/m",H)