0001 function [y] = fthetastar_dke_en(theta,rho,equil_fit,B0,ksi0,gamma) 0002 % 0003 0004 % 0005 y = NaN(1,length(theta)); 0006 % 0007 for itheta = 1:length(theta), 0008 % 0009 equil_val = equilval_yp(equil_fit,rho,theta(itheta)); 0010 % 0011 y(itheta) = equil_val.B/B0 - 2.0/(1 - ksi0*ksi0)/(gamma + 1.0); 0012 % 0013 end 0014