load '/your_path/local_defs.ijs' load '/your_path/MTsi.ijs' load '/your_path/s_and_p.ijs' load '/your_path/IQMsi.ijs' load '/your_path/IQs.ijs' NB. define tau grid: 1st point (after 0) is 0.05, last NB. point is 16, logarithmic spacing, 8 points per decade: tau=: tau_grid 0.05 8 16 NB. 23 points: 0 0.05 0.0666761 ... 11.8569 15.8114 16 L0=: 1 MTsi tau NB. get Lam, M & N transforms M0=: 3 MTsi tau N0=: 5 MTsi tau B=: (2%3) + tau NB. linear source term B(tau) lam=: 0.5 NB. constant lambda 's p'=: (L0;M0;N0) s_and_p (lam;B) NB. solve for s & p load 'plot' NB. load J plotting package plot tau;>s;10*p NB. plot s and 10 times p vs. tau NB. the above plot may be overwritten too quickly to see ... NB. define a grid of mu (= cos(theta)) values: mu=: 0.01 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.5 0.6 0.7 0.8 0.9 1 I0=: IQMsi t;mu NB. get matrix for I(0,mu) 'Im Qm pol'=: (I0;mu) IQs s;p NB. evaluate emergent radiation NB. plot I(mu) and -100 times the polarization as function of mu: plot mu;>Im;_100*pol