File:Pu241heatPlo800x517.jpg

This is low resolution duplicate of image http://mizugadro.mydns.jp/t/index.php/File:Pu241heatPlot.jpg

The original is horribly slow to load..

Decay heat by mixture of Pu-241 and Am-241, used in the Century heater, described in utopia Tartaria.

Abscise axis is time, measured in years. The initial concentration of Pu-241 is 100%.

C++ generator of curves

 * 1) include
 * 2) include
 * 3) include
 * 4) define DB double

void ado(FILE *O, int X, int Y) {      fprintf(O,"%c!PS-Adobe-2.0 EPSF-2.0\n",'%'); fprintf(O,"%c%cBoundingBox: 0 0 %d %d\n",'%','%',X,Y); fprintf(O,"/M {moveto} bind def\n"); fprintf(O,"/L {lineto} bind def\n"); fprintf(O,"/S {stroke} bind def\n"); fprintf(O,"/s {show newpath} bind def\n"); fprintf(O,"/C {closepath} bind def\n"); fprintf(O,"/F {fill} bind def\n"); fprintf(O,"/o {.1 0 360 arc C S} bind def\n"); fprintf(O,"/times-Roman findfont 20 scalefont setfont\n"); fprintf(O,"/W {setlinewidth} bind def\n"); fprintf(O,"/RGB {setrgbcolor} bind def\n");} //DB casde(DB k, DB p, DB t){ return k*p/(k-p)*(exp(-p*t)-exp(-k*t));}

int main{int n,m; DB year=365.25*24*60*60; //second

DB PuT, Puk, Pue; PuT=14.4 ; // Years Puk=log(2.)/PuT; Pue=0.021; //MeV

DB Am, AmT, Amk, Ame; AmT=432; Ame=5.638; Amk=log(2.)/AmT;

DB MeV=1.60218e-13; // Joule DB Year=365.242374*24*3600; //Second DB amu=1.660539040e-27 ; //kg

DB k=Puk/Year; //Herz DB p=Amk/Year; //Herz DB kp=k*p/(k-p); //Herz DB W=241.1*amu; //kg, mass of Pu atom DB E=Ame*MeV; //Joule DB ko=kp*E; //Watt DB mo=ko/W; //Watt/kg printf("kp=%12.4le W=%12.6le E=%12.6le mo=%12.6le\n",kp,W,E,mo);

FILE *o; o=fopen("11.eps","w"); ado(o,1620,1040); fprintf(o,"10 10 translate 10 10 scale 2 setlinecap 1 setlinejoin\n"); M(0,104) L(0,0) L(162,0) fprintf(o,".5 W S\n");
 * 1) define M(x,y) fprintf(o,"%6.4lf %6.4lf M\n",0.+x,0.+y);
 * 2) define L(x,y) fprintf(o,"%6.4lf %6.4lf L\n",0.+x,0.+y);

for(n=10; n<101; n+=10){ M(0,n) L(160,n)} for(m=10; m<161; m+=10){ M(m,0) L(m,100)} fprintf(o,".3 W S\n");

DB t; for(n=0;n<164;n+=2){t=n; DB c=exp(-Amk*t)-exp(-Puk*t); DB d=c*mo; //Watt/kg DB f=d/E; //Bq/kg printf("%3.0f %6.4lf %9.3le %9.4lf\n",t,c,f,d); if(n==0) M(t,d) else L(t,d) } fprintf(o,"1 0 0 RGB .5 W S\n");

fprintf(o,"showpage\n"); fclose(o);

system("epstopdf 11.eps"); system( "open 11.pdf");

}

Latex generator of labels
\documentclass[12pt]{article} \usepackage{geometry} \topmargin -96pt \oddsidemargin -64pt \paperwidth 1680pt \paperheight 1086pt \textwidth 1600pt \textheight 1200pt \usepackage{amsmath} \usepackage{graphicx} \usepackage{rotating} \parindent 0pt \newcommand{\sx}{\scalebox} \newcommand \ing {\includegraphics} \newcommand \rot {\begin{rotate}} \newcommand \ero {\end{rotate}} \begin{document} \begin{picture}(1700,1054) \put(50,30){\ing{11}} \put(-9,1026){\sx{5}{$P,$}} \put(-5,958){\sx{5}{$\rm \frac{W}{kg}$}} %\put(-8,924){\sx{5}{$90$}} \put(-8,824){\sx{5}{$80$}} %\put(-8,724){\sx{5}{$70$}} \put(-8,624){\sx{5}{$60$}} %\put(-8,524){\sx{5}{$50$}} \put(-8,424){\sx{5}{$40$}} %\put(-8,324){\sx{5}{$30$}} \put(-8,224){\sx{5}{$20$}} %\put(-8,124){\sx{5}{$10$}} \put(4, 24){\sx{5}{$0$}} \put(46,-9){\sx{5}{$0$}} %\put(128,-9){\sx{5}{$10$}} \put(228,-9){\sx{5}{$20$}} \put(428,-9){\sx{5}{$40$}} \put(628,-9){\sx{5}{$60$}} \put(828,-9){\sx{5}{$80$}} \put(1020,-9){\sx{5}{$100$}} \put(1220,-9){\sx{5}{$120$}} \put(1420,-9){\sx{5}{$140$}} \put(1584,-6){\sx{5}{$t,$ y}} \end{picture} \end{document}