exm020804.m

matlab电子书籍

syms t x
ft=1/log(t);
sx=int(ft,t,0,x)  

x=0.1:0.05:0.9;							
sx_n=-mfun('Ei',1,-log(x));		 

clf
ezplot('1/log(t)',[0.1,0.9])		
grid on
hold on
plot(x,sx_n,'LineWidth',3)			
Char1='1/ln(t)';
Char2='{\int_0^x} 1/ln(t) dt';	
title([Char1,'   and    ',Char2])		
legend(Char1,Char2,'Location','SouthWest')