//Find critical current
//Ex:16.3
clc;
clear;
close;
r=1*10^-3;//in m
hc=7.9*10^3;//in A/m
ic=2*3.14*r*hc;//in m
disp(ic,"Critical current in superconducting state (in A) = ");