file saves binom.wxmx but will not come back up

[zmth]

again and again it does this file saves ok and then will NOT come back up and i spend hours and hours then to have the file not load says xml parsing error. Then why did it save it ! Can you please fix this file. by giving me somewhere to email it.

[gunterkoenigsmann]

Wow... ...such a thing didn't happen for a long time. Which wxMaxima version are you using? If it is a version from<2015 chances are high that the current nightly build of wxMaxima will automatically fix the issue. If not the process of repairing nearly every possible error is:

• Make a copy of the file,
• rename the copy to something ending in .zip,
• unzip it
• open the content.xml with an XML-aware text editor
• and hope that your text editor instantly informs you where the error lies.

The current nightly build of wxMaxima will be able to open the content.xml, after it has been fixed.

If that doesn't help you can still open the .wxmx file in an text editor. And my mail address is wxMaxima at peterpall . de.

If you are willing to send me the broken file I am grateful as with a bit of luck it will help me resolving the actual problem that caused the file to corrupt.

Unfortunately my computer currently isn't in working order so I will only be able to repair the file by hand in approximately a week.

[gunterkoenigsmann]

What is the exact error message? The things I have ever seen were:

• an illegal character in the file. The current nightly build should fix that kind of error whilst opening the file.
• a truncated file due to a broken hard disk. Normally the content.xml file can be recovered in this case using an text editor as it is at the beginning of the .wxmx file and isn't compressed.
• a missing closing XML tag
• and a completely garbled output tag. Which isn't a problem as if you delete output tags including the matching closing /output using the text editor the content.xml cannot still be opened by the nightly build of wxMaxima.

[zmth]

thanks for your reply: wxMaxima version: 13.4.0 Maxima version: 5.31.1 Maxima build date: 2013-09-24 09:49:12 Host type: i686-pc-mingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8

before years ago same thing hapenned and they fixed it easily. I don't have an xml editor don't even know what XML means. Don't have any way to send it to peterpall . de as is not a valid email address. i renamed it binom.zip and will try to upload here binom.zip thanks

[zmth]

error says " wxMaxima encountered an error loading C:\maxima\binom.wxmx " then i press ok and it replies " XML parsing error: 'not well-formed (invalid token)' at line 828 "

I brought up the file in wordpad but unfortunately there is nothing in wordpad that can show line numbers

[zmth]

well ok here is a copy of the file made into a file 'content.xml': <?xml version="1.0" encoding="UTF-8"?>

(assume(positive,true),assume_pos_true,showtime:true,ratprint:false,load("lrats"),simpsum:true ,load(simplify_sum),file_output_append:true );

Evaluation took 3.6800 seconds (3.6800 elapsed) (%o1) true

(kill(a,a1,b,belis,blis,i,ii,in,lisr,listmp,lisx,ls1m,mlis,mm,n,nps,nt,ntop,ntp,p,pb,plis,pmb,pnps,t,tt,u2lis,u2ls) ,pmb(bblis):=block([fl1:false,jj,tmp],pb:makelist(i,i,1,nps), i0:1,loop, jj:0,if bblis[pb[i0+1]] >bblis[pb[i0]] then(fl1:true,for j thru nps-1-i0 while bblis[pb[i0+1+j]]>bblis[pb[i0+j]] do jj:j,for i:0 thru floor(jj/2) do (tmp:pb[i0+i],pb[i0+i]:pb[i0-i+jj+1],pb[i0-i+jj+1]:tmp)/*end i*/) /*end then*/,if i0+jj<nps-1 then(i0:i0+jj+1,go (loop)) else if fl1 then (fl1:false,i0:1,go (loop)) , blis:makelist(bblis[pb[i]],i,1,nps),for i:0 thru nps-2 while blis[nps-i]=0 do blis:rest(blis,-1) , d2lis:makelist(0,nps),pai:makelist(sum(p[pb[j]]^2/a[pb[j]],j,1,i-1)+sum(p[pb[j]]^2/a[pb[j]],j ,i+1,nps),i,1,nps) )/*end blk*/ /*bblis must be atleast length 2 */ ,tpd2(i,im,mm):=block([tt],if i<ii then tpd2(i+1,im,mm) else u2lis[ta+mm]:u2lis[ta+mm]+im*tb ,for m thru d2lis[i] do(tt:d2lis[i]!/m!/(d2lis[i]-m)!*pai[i]^m,if i<ii then tpd2(i+1,tt*im,mm+m) else u2lis [ta+mm+m]:u2lis[ta+mm+m]+im*tb*tt) )/*end blk*/ ,d2enum(blis,ii,in,ntop):=block([n1:blis[ii] ,n1r,nnt],if ii=length(blis) then(d2lis[ii]:(n1-ntop)/2,ta:n/2 -sum(d2lis[i],i,1,ii),tb:in*prod(n1-i,i,0,ntop-1)*abs(n1-ntop-1)!!*prod(a[pb[i]]^(d2lis[i]-blis[i]),i,1,ii) *(-1)^sum(d2lis[i],i,1,ii), for i thru ii do if blis[i]-2*d2lis[i]>0 then tb:tb*p[pb[i]]^(blis[i]-2*d2lis[i]) ,tpd2(1,1,1)) else(n1r:sum(blis[i],i,ii+1,length(blis)),for ntp:max(0,ntop-n1r) thru min(ntop,n1) do (nnt:prod(n1-i,i,0,ntp-1)*binomial(ntop,ntp), for d2:max((n1+ntop-n1r)/2-ntp,0) thru (n1-ntp)/2 do (d2lis[ii]:d2,d2enum(blis,ii+1,in*abs(2*d2-1)!!*binomial(n1-ntp,2*d2)*nnt,ntop+n1-2*(d2+ntp) ) )/*end d2*/)/*end ntp*/)/*end else*/)/*end blk*/ ,PEM(nps,a,belis,p):=block([a1,u2ls:[],nt:0,t], for k:0 thru 2 do(n:0, for i thru nps do n:n+belis[k*nps+i] , if n>0 then(u2lis:makelist(0,n/2+1),pmb(rest(rest(belis,(k-2)*nps),k*nps)),nt:nt+n,d2enum(blis,1,1,0) ) else u2lis:[1], if k=0 then u2ls:append(u2ls,u2lis) else(u2ls:append (u2ls,makelist(0,n/2)),for i:0 thru nt/2 do u2ls[nt/2+1-i]:sum(u2ls[nt/2+1-n/2-j]*u2lis[n/2+1-i+j],j,max(0 ,i-n/2),min(i,(nt-n)/2)))/*endelse*/)/*end k*/ ,a1:sum(p[i]^2/a[i],i,1,nps) ,t:u2ls[nt/2+1], for i:0 thru nt/2-1 do t:u2ls[nt/2-i]+t*(nt-2*i-1)/a1/2/(i+1) , t:2*%pi^((3*nps-1)/2)*t*(-1)^(nt/2)*(nt/2)!/prod((2*i+1),i ,1,nt/2)/sqrt(a1)/prod(a[i],i,1,nps)^(3/2), return(t))/*end blk*/ ,kill(bdlis,i,iej,ir,j1,j2,j,m1,m2,m,n,np,nps,dout,ii,lisin,ls,mn,nps,ts,tm,t2,zd1,zd2) /*HERE try direct cj of edmonds... ANd find identical result*/ /* ,vc(j1,m1,j2,m2,j):=block([t],if j<abs(j1-j2) or j>j1+j2 or j<abs(m1+m2) then return(0) else(t:(-1)^m1 *sqrt(factor( /*(2*j+1)*(j1+j2-j)!*/(j1-m1)!*(j2-m2)!*(j+m1+m2)!*(j-m1-m2)! /(/*(j+j1+j2+1)!*(j +j1-j2)!*(j-j1+j2)!*/(j1+m1)!*(j2+m2)!)))*sum((-1)^s*(j1+m1+s)!*(j2+j-m1-s)! /( s!*(j1-m1-s)!* (j-m1-m2-s)!*(j2-j+m1+s)! ),s,max(0,j-j2-m1), min(j1-m1,j-m1-m2)),return(t))/*end else*/ )/*end blk*/ */ ,jkz(i,ir,n,np):=block([tm],if i=1 then lisi:makelist(0,2*n-min(3,n)),if i<n+1 then(if i<n then(if i<n-1 then for ii:0 thru min(np,2*np-ir) do(lisi[i]:ii,jkz(i+1,ir+ii,n,np))/*end ii*/ else for ii:max(0,np-ir) thru min(np,2*np-ir) do(lisi[i]:ii,jkz(i+1,ir+ii,n,np))/*end ii*/)/*end i<n then*/ else(lisi[n]:2*np-ir,if n<4 then liss:append(liss, [copylist(lisi)])/*MUSt use copylist*/ else jkz(n+1,ir+lisi[n]-lisi[1]-lisi[2],n,np))/*end else*/)/*end i<n +1 then*/ else (tm:lmax(rest(makelist(lisi[i],i,3,n),i-n-1)),if i>n+1 then for ii:max(2*tm-ir,abs(lisi[i-1]- lisi[i+1-n])) thru min(ir,lisi[i-1]+lisi[i+1-n]) do(lisi[i]:ii,if i<2*n-3 then jkz(i+1,ir-lisi[i+2-n],n,np) else liss:append(liss,[copylist(lisi)]))/*end ii*/ else for ii:max(2*tm-ir,abs(lisi[2]-lisi[1])) thru min(ir,lisi[1]+ lisi[2]) do(lisi[n+1]:ii,if i<2*n-3 then jkz(n+2,ir-lisi[3],n,np) else liss:append(liss,[copylist(lisi)])) /*end ii*/)/*end else*/)/*end blk*/ ,ad3n(ii,nn,lisin):= if ii<nps then for i:0 thru nn do ad3n(ii+1,nn-i,makelist(lisin[j]+(1-abs(signum(ii-j))) *i,j,1,nps)) else(lisin[ii]:lisin[ii]+nn,lisout:append(lisout,[lisin])) /* ,wsd(ii,lisin,mn,ls):=block([j1,j2,j3,ni,liso],if ii=nps+1 then(j1:ls[1],j2:ls[2],ni:2) else(if ii<length(ls)+2 then j1:ls[ii-1],j2:ls[ii+1-nps],ni:ii+1-nps),if ii=length(ls)+2 then(j3:0,j1:j2) else if ii=length(ls)+1 then j3:ls[nps] *signum(nps-2) else j3:ls[ii],if ii> length(ls)+signum(nps-2) then(lout:[],dout:[]), for m:max(-j1,mn-j2) thru min(j1,mn+j2) do(tt:vc(j1,m,j2,mn-m,j3),if tt#0 then(liso:copylist(lisin),liso[ni]:mn-m,liso[nps+1]:radcan (liso[nps+1]*tt),if ii>nps+1 then wsd(ii-1,liso,m,ls) else(liso[1]:m, lout:append(lout,[copylist(liso)])) /*end else*/ )/*endthen*/ )/*endm*/)/*end blk*/ */ ,szp(iej,zd1,zd2):=block([tn:0], for i thru length(zd1) do for j:1+iej*(i-1) thru length(zd2) do tn:tn+(1+iej *signum(j-i))*zd1[i][nps+1]*zd2[j][nps+1]*prod((2*(ls[k]+zd2[j][k]+zd1[i][k])+1)!!/2^(ls[k]+zd1[i][k] +zd2[j][k]+2)/a[k]^(zd1[i][k]+zd2[j][k]+ls[k]+3/2),k,1,nps)*%pi^(nps/2)*(1+sum(km[k]*a[k]*(4*(zd1 [i][k]+zd2[j][k]-(zd2[j][k]-zd1[i][k])^2+ls[k]*(ls[k]+2))+3)/(4*(ls[k]+zd1[i][k]+zd2[j][k])+2),k,1,nps)) ,return(tn))/*end blk*/)$

Evaluation took 0.0100 seconds (0.0100 elapsed)

/*Need before call lm12r(1,1) or lm12nc(1,1) : nps, lisr:[],ipx:[1,2,3], lisr:[], listmp:makelist(0,3*nps) ,lisp1p2:makelist(0,2*nps-2),belis:makelist(0,3*nps) , mlis[1],mlis[2],ls1[1],ls1[2],pnps like irps of fortran 7/Apr/2018 21:52 */ (kill(belis,lisr,listmp,lisx,ipx,plis,pnps,ls1,mlis) ,lm12nc(ii,rr):=block([tt,tk,rp:pnps[ii],tlm1,tlm2,m1:abs(mlis[1][ii]),m2:abs(mlis[2][ii]),l1:ls1[1][ii],l2 :ls1[2][ii],fl,fl1,fl2,p1b,p2b,stp] ,if ii<nps then(stp:1,p1b:0,p2b:0) else (stp:2, p1b:mod(sum(lisp1p2[i] ,i,1,nps-1),2), p2b:mod(sum(lisp1p2[i],i,nps,2*nps-2),2)),for k1:0 thru (l1-m1)/2 do for p1:p1b thru 2*k1 +m1 step stp do(tlm1:sum(binomial(k1,j)*binomial(m1,p1-2*j)*(-1)^j,j,max(floor((p1-m1+1)/2),0) ,min(k1,floor(p1/2))),if tlm1#0 then(if mlis[1][ii]#m1 then tlm1:tlm1:tlm1*(-1)^(m1+p1), if ii<nps then lisp1p2[ii]:p1,for k2:0 thru (l2-m2)/2 do for p2:p2b thru 2*k2+m2 step stp do(tlm2:sum (binomial(k2,j)*binomial(m2,p2-2*j)*(-1)^j,j,max(floor((p2-m2+1)/2),0),min(k2,floor(p2/2))) ,if tlm2#0 then(if ii<nps then lisp1p2[nps-1+ii]:p2,if mlis[2][ii]#m2 then tlm2:tlm2*(-1)^(m2+p2) /*outer scalar loop via i1,i2 follows per above rp:pnps[ii], ordering listmp particle 1's 3 rect coord's next particle 2,..., nps*/ , for i1:0 thru rp do for i2:0 thru rp-i1 do (lisx:[2*(k1+k2+i1)+m1+m2-p1-p2, p1+p2+2*i2, 2*(rp -i1-i2-k1-k2)+ l1+l2-m1-m2], listmp[3*ii-2]:lisx[1],listmp[3*ii-1]:lisx[2], listmp[3*ii]:lisx[3], tk :sqrt((2*l1+1)*(l1+m1)!*(l1-m1)!*(2*l2+1)*(l2+m2)!*(l2-m2)!)*(-1)^(k1+k2+m1+m2)/2^(2*(k1 +k2)+m1+m2)/(m1+k1)!/k1!/(l1-m1-2*k1)!/(m2+k2)!/k2!/(l2-m2-2*k2)!*tlm1*tlm2*rp!/i1!/i2! /(rp-i1-i2)! , if ii<nps then lm12nc(ii+1,rr*tk) else (/*tk:tk/4^nps/%pi^nps*/tk:tk*(-1)^((sum(lisp1p2 [i],i,1,2*nps-2)+p1+p2)/2),jb:0, for i thru 2 do(fl2:true,ip:i, for j:jb thru nps-1 while fl2 do(imx:listmp [3*j+ipx[i]],for ii:i+1 thru 3 do(it:listmp[3*j+ipx[ii]], if fl2 then if it#imx then (fl2:false,jb:j), if it >imx then(imx:it, ip:ii)), if ip>i then(itmp:ipx[i], ipx[i]:ipx[ip],ipx[ip]:itmp),if not fl2 and i=2 and j<nps -1 and listmp[3*j+ipx[1]]=listmp[3*j+ipx[2]] then( j1:j ,for i:j+1 thru nps-1 while listmp[3*i+ipx[1]] =listmp[3*i+ipx[2]] do j1:i, if j1<nps-1 then (j1:j1+1, if listmp[3*j1+ipx[1]]<listmp[3*j1 +ipx[2]] then(itmp:ipx[1],ipx[1]:ipx[2],ipx[2]:itmp))))) ,for i thru 3*nps do belis[floor((i-1)/3)+mod((i-1),3)*nps+1] :listmp[ipx[1+mod((i-1),3)]+3*floor((i-1)/3)] ,fl1:true, for i thru length(lisr) while fl1 do if belis=rest(lisr[i],-1) then(fl1:false ,lisr[i][3*nps+1]:lisr[i][3*nps+1]+float(rr*tk) )/*end then*/ ,if fl1 then lisr:append(lisr,[append(belis,[float(rr*tk)])]) )/*end else*/) /*end i2*/)/*end then*/)/*end p2*/)/*end then*/)/*end p1*/ )/*end blk*/ )/*end NO complex lm12r renamed lm12nc*/ $

Evaluation took 0.0100 seconds (0.0100 elapsed)

float((4*%pi)^4);

Evaluation took 0.0000 seconds (0.0000 elapsed) (%o14) 24936.72730470462

(nps:2, lisr:[], ipx:[1,2,3], listmp:makelist(0,3*nps) ,lisp1p2:makelist(0,2*nps-2),belis:makelist(0,3*nps) ,pnps:[0,14] ,ls1[1]:makelist(0,nps) ,ls1[2]:ls1[1], mlis[1]:ls1[1] , mlis[2]:ls1[1] , lm12nc(1,1),ldisp(["len(lisr)",length(lisr),"lisr",lisr]) );

(%t10) [len(lisr),24,lisr,[[0,28,0,0,0,0,3.0],[0,26,0,2,0,0,84.0],[0,24,0,4,0,0,546.0],[0,22,0,6,0,0,2184.0],[0,20,0,8,0,0,6006.0],[0,18,0,10,0,0,12012.0],[0,16,0,12,0,0,18018.0],[0,14,0,14,0,0,10296.0],[0,24,0,2,0,2,546.0],[0,22,0,4,0,2,6552.0],[0,20,0,6,0,2,24024.0],[0,18,0,8,0,2,60060.0],[0,16,0,10,0,2,108108.0],[0,14,0,12,0,2,144144.0],[0,20,0,4,0,4,18018.0],[0,18,0,6,0,4,120120.0],[0,16,0,8,0,4,270270.0],[0,14,0,10,0,4,432432.0],[0,12,0,12,0,4,252252.0],[0,16,0,6,0,6,180180.0],[0,14,0,8,0,6,720720.0],[0,12,0,10,0,6,1009008.0],[0,12,0,8,0,8,630630.0],[0,10,0,10,0,8,756756.0]]] Evaluation took 0.6100 seconds (0.6100 elapsed) (%o10) [%t10]

(fl:true,kill(a,km),file_output_append:true,lisnp:[],icnt:0,lsa:[] ,nps:2,npmx:7 ,lisin:makelist(0,nps), ad3al:[lisin], for ihz thru 2*npmx do(lisin:makelist(0,nps) ,lisout:[],ad3n(1,ihz,lisin),for iz thru length(lisout) do ad3al :append(ad3al,[copylist(lisout[iz])])) /*,ldisp(["length(ad3al)",length(ad3al)," ad3al",ad3al])*/);

Evaluation took 0.0300 seconds (0.0300 elapsed) (%o3) done

(nps:2,alis:[1.03,1.03], qq:[1,-2,-2], pmat:matrix([-1,1],[1,0],[0,1]) ,rirps:makelist(float(3*10^38),length(ad3al)) ,ii0:float(16*%pi^2) , for i thru length(ad3al) do rirps[i]: -integrate(integrate(r1^(1+2*ad3al[i][1])*exp(-alis[1]*r1^2) *r2^(2+2*ad3al[i][2])*exp(-alis[2]*r2^2),r1,r2,inf),r2,0,inf) -integrate(integrate(r1^(2+2*ad3al[i][1])*exp(-alis[1]*r1^2) *r2^(1+2*ad3al[i][2])*exp(-alis[2]*r2^2),r2,r1,inf),r1,0,inf) +2*integrate(r^(1+2*ad3al[i][1])*exp(-alis[1]*r^2),r,0,inf) *integrate(r^(2+2*ad3al[i][2])*exp(-alis[2]*r^2),r,0,inf) +2*integrate(r^(2+2*ad3al[i][1])*exp(-alis[1]*r^2),r,0,inf) *integrate(r^(1+2*ad3al[i][2])*exp(-alis[2]*r^2),r,0,inf) ,rirps:float(ii0*rirps), for i thru length(ad3al) do stringout("lm12.t",[ad3al[i], i, rirps[i] ]) );

Evaluation took 16.4300 seconds (16.4300 elapsed) (%o384) done

(tp:PEM(2,[1,1],[0,0, 0,0, 0,0],[1,0]), tt:16*%pi^2*integrate(r*exp(-r^2),r,0,inf)*integrate(r^2*exp(-r^2) ,r,0,inf),ldisp(["tp",tp,"ck 0",tp-tt]) );

(%t20) [tp,2*%pi52,ck 0,0] Evaluation took 0.0200 seconds (0.0200 elapsed) (%o20) [%t20]

/*for any positive cst and any positive alis ck's 0. Reason for the .../c is 1/r dependence is c times as far away, that is for the potential of a test charge due to a point charge at origin assuming product of test and point charges being unity*/ (kill(alis),cst:.1,alis:[3,5],for c:cst thru 4 do(tp:PEM(2,alis,[0,2, 0,0, 0,0], [c,0]) ,tt:16*%pi^2*integrate(r*exp(-alis[1]*r^2),r,0,inf)*integrate(r^4*exp(-alis[2]*r^2) ,r,0,inf)/c , ldisp(["c",c,"tp",float(tp),"ck 0",float(radcan(tp-tt/3))])))$ /*for any positive cst and any positive alis*/ (kill(alis),cst:1,alis:[3,5],for c:cst thru 4 do(tp:PEM(2,alis,[0,2, 0,0, 0,0], [0,-c]), tt:16*%pi^2 *integrate(r^2*exp(-alis[1]*r^2),r,0,inf)*integrate(r^3*exp(-alis[2]*r^2),r,0,inf)/c ,ldisp(["c",c,"tp",float(tp),"ck 0",float(radcan(tp-tt/3))])))$

(%t125) [c,1,tp,0.089776264715563,ck 0,0.0] (%t126) [c,2,tp,0.044888132357781,ck 0,0.0] (%t127) [c,3,tp,0.029925421571854,ck 0,0.0] (%t128) [c,4,tp,0.02244406617889,ck 0,0.0] Evaluation took 0.1000 seconds (0.1000 elapsed)

(kill(alis),cst:1,alis:[3,5],for c:cst thru 4 do(tp:PEM(2,alis,[0,2, 0,0, 0,0], [-c,c]), tt:16*%pi^2 *(integrate(integrate(r1*exp(-alis[1]*r1^2)*r2^4*exp(-alis[2]*r2^2),r1,r2,inf),r2,0,inf) +integrate(integrate(r1^2*exp(-alis[1]*r1^2)*r2^3*exp(-alis[2]*r2^2),r2,r1,inf),r1,0,inf))/c ,ldisp(["c",c,"tp",float(tp),"tt",float(tt),"ck 0",float(radcan(tp-tt/3))])))$ (kill(alis),c2:2,cst:1,alis:[3,5],for c:cst thru 4 do(tp:PEM(2,alis,[2,0, 0,0, 0,0], [c,-c]), tt:16*%pi^2 *(integrate(integrate(r1^3*exp(-alis[1]*r1^2)*r2^2*exp(-alis[2]*r2^2),r1,r2,inf),r2,0,inf)/c +integrate(integrate(r1^4*exp(-alis[1]*r1^2)*r2*exp(-alis[2]*r2^2),r2,r1,inf),r1,0,inf)/c) ,ldisp(["c",c,"tp",float(tp),"tt",float(tt),"ck 0",float(radcan(tp-tt/3))])))$

(%t219) [c,1,tp,0.108807675827,tt,0.32642302748102,ck 0,0.0] (%t220) [c,2,tp,0.054403837913503,tt,0.16321151374051,ck 0,0.0] (%t221) [c,3,tp,0.036269225275668,tt,0.108807675827,ck 0,0.0] (%t222) [c,4,tp,0.027201918956751,tt,0.081605756870255,ck 0,0.0] Evaluation took 0.2700 seconds (0.2700 elapsed)

(kill(alis),cst:1,alis:[3,5],for c:cst thru 4 do(tp:PEM(2,alis,[2,2, 0,0, 0,0], [c,-c]) +PEM(2,alis,[2,0, 0,2, 0,0], [c,-c])+PEM(2,alis,[2,0, 0,0, 0,2], [c,-c]), tt:16*%pi^2 *(integrate(integrate(r1^3*exp(-alis[1]*r1^2)*r2^4*exp(-alis[2]*r2^2),r1,r2,inf),r2,0,inf) +integrate(integrate(r1^4*exp(-alis[1]*r1^2)*r2^3*exp(-alis[2]*r2^2),r2,r1,inf),r1,0,inf))/c ,ldisp(["c",c,"tp",float(tp),"tt",float(tt),"ck 0",float(radcan(tp-tt/3))])))$ (kill(alis),cst:1,alis:[3,5],for c:cst thru 4 do(tp:PEM(2,alis,[2,2, 0,0, 0,0], [c,-c]) +PEM(2,alis,[0,2, 2,0, 0,0], [c,-c])+PEM(2,alis,[0,2, 0,0,2,0], [c,-c]), tt:16*%pi^2 *(integrate(integrate(r1^3*exp(-alis[1]*r1^2)*r2^4*exp(-alis[2]*r2^2),r1,r2,inf),r2,0,inf) +integrate(integrate(r1^4*exp(-alis[1]*r1^2)*r2^3*exp(-alis[2]*r2^2),r2,r1,inf),r1,0,inf))/c ,ldisp(["c",c,"tp",float(tp),"tt",float(tt),"ck 0",float(radcan(tp-tt/3))])))$ (kill(alis),cst:1,alis:[3,5],for c:cst thru 4 do(tp:2*PEM(2,alis,[0,2, 0,2, 0,0], [c,-c]) +PEM(2,alis,[0,0, 0,0, 0,4], [c,-c]), tt:16*%pi^2 *(integrate(integrate(r1*exp(-alis[1]*r1^2)*r2^6*exp(-alis[2]*r2^2),r1,r2,inf),r2,0,inf) +integrate(integrate(r1^2*exp(-alis[1]*r1^2)*r2^5*exp(-alis[2]*r2^2),r2,r1,inf),r1,0,inf))/c ,ldisp(["c",c,"tp",float(tp),"tt",float(tt),"ck 0",float(radcan(tp-tt/3))])))$ (nps:2,em:float((4*%pi)^nps));

Evaluation took 0.0000 seconds (0.0000 elapsed) (%o3) 157.9136704174297

Fotran compare /*11/May/2018 10:05*/ (nps:4,npmx:3, alis:[.35, 14, .8, .8],qlis:[1, 1, -1, -1, 1] ,pmat:matrix([-0.5, 1, 0, 0],[-0.5, .333333313,1,0],[-0.5, .333333313,1.81506315E-04,1] ,[1, 0, 0,0],[0, -0.666666627,1, 0],[0, -0.666666627, 1.81506315E-04,1] ,[.5,1,0,0],[0,0,-0.999818504,1],[.333333313,0.333333313,1,0] ,[0.333272815, 0.333272815, 1.81506315E-04,1]) ,ts:0, lisr:[], ipx:[1,2,3], listmp:makelist(0,3*nps) ,lisp1p2:makelist(0,2*nps-2),belis:makelist(0,3*nps) ,pnps:[2,0,0,0] ,ls1[1]:makelist(0,nps) ,ls1[2]:ls1[1], mlis[1]:ls1[1] , mlis[2]:ls1[1] , lm12nc(1,1),ldisp(["len(lisr)",length(lisr),"lisr",lisr]) ,for i thru length(lisr) do(tmp:0,lisin:rest(lisr[i],-1) ,for i thru nps do for j:i+1 thru nps+1 do (rp:pmat[(2*i*nps-2*nps-i**2+i-2)/2+j] t:float(PEM(nps,alis,lisin,rp))/*,ldisp([i,j," t",t])*/ ,tmp:tmp+qlis[i]*qlis[j]*t) , ts:ts-rest(lisr[i],3*nps)*tmp),ts );

(%t59) [len(lisr),2,lisr,[[4,0,0,0,0,0,0,0,0,0,0,0,3.0],[2,0,0,0,2,0,0,0,0,0,0,0,6.0]]] Evaluation took 0.1600 seconds (0.1600 elapsed) (%o59) [12358.68808231463]

/* 16/Apr/2018 15:02R */ (nps:4, alis:[.35 , 14 , .8 , .8 ] ,er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[5.44617767E-04, 4.08463355E-04, 0.500090778, 0.500090778] ,ls:[0,0,0,0],irps1lis:[[0,0,0,0],[0,0,0,1],[0,0,1,0],[0,1,0,0],[1,0,0,0] ] , irps2lis:[[0,0,0,0]] /*irps1lis*/ ,nm1lis:[-9.88848782,9.77482123E-05,9.69184839E-05,1.56217811E-05 ,2.30744338] , nm2lis:[8.07506847] /*nm1lis*/ ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm);

Evaluation took 0.7000 seconds (0.7000 elapsed) (%o5) 3.099127450223449*10-7

/*13/Apr/2018 14:43 ERROR */ (nps:4, alis:[1.44 , 1.44 , .830 , .830 ] ,er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[5.44617767E-04, 4.08463355E-04, 0.500090778, 0.500090778] , ls:[3,0,3, 0],irps1lis:[[0,0,0,0]],irps2lis:[[0,0,0,0]] ,nm1lis:[0.75953704] , nm2lis:nm1lis ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm);

Evaluation took 0.1400 seconds (0.1400 elapsed) (%o4) 4.593492688239066

/*4/Apr/2018 1:07*/ (kill(alis),alis:[1.03,1.03],t0:integrate(r *exp(-alis[2]*r^2),r,0,inf) ,for i thru 7 do(tt:integrate(r^(2*i+1) *exp(-alis[2]*r^2),r,0,inf) ,ldisp(["i",i,"tt/t0",float(radcan(tt/t0))])))$

(%t31) [i,1,tt/t0,0.97087378640776] (%t32) [i,2,tt/t0,1.885191818267508] (%t33) [i,3,tt/t0,5.490849956118957] (%t34) [i,4,tt/t0,21.32368914997655] (%t35) [i,5,tt/t0,103.5130541260996] (%t36) [i,6,tt/t0,602.9886648122311] (%t37) [i,7,tt/t0,4097.981217170502] Evaluation took 0.0800 seconds (0.0800 elapsed)

(nps:3,id:84,mm:ident(id), read_matrix("braketmat.t",mm) ,eij:eigens_by_jacobi(mm),m2:part(eij,2) ,eival:part(eij,1),m1:mm.m2,ms:transpose(m2).m1 /*,ldisp(["mm",mm,"eival",eival,"m2",m2,"m1",m1,"ms",ms])*/ ,[float((lmax(eival)/(4*%pi)^nps)^2),float((lmin(eival)/(4*%pi)^nps)^2)] );

Warning: WITH-COMPILATION-UNIT is being redefined. Evaluation took 60.1600 seconds (60.1600 elapsed) (%o3) [0.80464391760693,0.0073504853167731]

(nps:2 , alis:[1.03,1.03],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[0.500136137, 0.500136137] , ls:[0, 0],irps1lis: [[ 0 , 0 ], [ 0 , 1 ], [ 1 , 0 ], [ 0 , 2 ], [ 1 , 1 ], [ 2 , 0 ], [ 0 , 3 ], [ 1 , 2 ], [ 2 , 1 ], [ 3 , 0 ], [ 0 , 4 ], [ 1 , 3 ], [ 2 , 2 ], [ 3 , 1 ], [ 4 , 0 ], [ 0 , 5 ], [ 1 , 4 ], [ 2 , 3 ], [ 3 , 2 ], [ 4 , 1 ], [ 5 , 0 ], [ 0 , 6 ], [ 1 , 5 ], [ 2 , 4 ], [ 3 , 3 ], [ 4 , 2 ], [ 5 , 1 ], [ 6 , 0 ], [ 0 , 7 ] ] ,irps2lis:irps1lis,nm1lis:[ -0.984453917, 4.98551035 , -2.04678550E-02 , -6.50452614 , 7.05079511E-02 , 3.07178684E-03 , 3.37654257 , -6.03787862E-02 , -6.91111479E-03 , -4.04286227E-04 , -0.820529819 , 1.90783106E-02 , 3.78440111E-03 , 6.11288880E-04 , 3.08473573E-05 , 9.82600600E-02 , -2.41830479E-03 , -6.73635164E-04 , -2.25071999E-04 , -2.23464795E-05 , -1.53048086E-06 , -5.55719808E-03 , 1.03745107E-04 , 3.91982539E-05 , 1.77588245E-05 , 5.51672747E-06 , 1.70829168E-07 , 4.33522267E-08 , 1.17456002E-04] ,nm2lis:nm1lis,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate (r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps) *tket -sum(rkm[i]*(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]) , i,1,nps)),tij:tbra*tij, for i thru nps do tij:integrate (r[i]^2*tij,r[i],0,inf),nm:float(nm+nm1lis[i1] *nm2lis[i2]*(tad+tij))),nm);

Evaluation took 47.3800 seconds (47.3800 elapsed) (%o10) 1.019263295209839

/*ck nps 3 00 scalar diag*/ (nps:3, alis:[.1,.94,.94],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[5.44617767E-04, 0.500136137, 0.500136137] , ls:[4,2,2],irps1lis:[[0,0,0]],irps2lis:[[0,0,0]] ,nm1lis:[8.13202496E-05],nm2lis:nm1lis ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm);

Evaluation took 0.1600 seconds (0.1600 elapsed) (%o15) 0.99999979220816

(nps:3,c11:.353890479,c21 :-0.302185416, c22:0.235252231 , alis:[.1,.94,.94],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[5.44617767E-04, 0.500136137, 0.500136137] , ls:[0, 0,0],irps1lis:[[0,0,0]],irps2lis:[[0,0,0],[0,0,1]] ,nm1lis:[ c11],nm2lis:[ c21,c22] ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm);

Evaluation took 0.2100 seconds (0.2100 elapsed) (%o9) 4.714530010829776*10-8

(nps:2,c11:1.42729914,c21 :-1.20808673, c22:1.00673914 ,alis:[1,1], qq:[1,-2,-2], pmat:matrix([-1,1],[1,0],[0,1]) ,bkmat:zeromatrix(nps,nps) ,tt:sum(-qq[i]*PEM(2,alis,[0,0, 0,0, 0,0],pmat[i]),i,1,nps*(nps+1)/2) ,bkmat[1,1]:float(tt*c11^2)/*agree fortran*/ ,i00p1:16*%pi^2*(integrate(integrate(r1*exp(-alis[1]*r1^2) *r2^2*exp(-alis[2]*r2^2),r1,r2,inf),r2,0,inf) +integrate(integrate(r1^2*exp(-alis[1]*r1^2)*r2*exp (-alis[2]*r2^2),r2,r1,inf),r1,0,inf)) ,i00p2:16*%pi^2*integrate(r^2*exp(-alis[1]*r^2),r,0,inf) *integrate(r*exp(-alis[2]*r^2),r,0,inf),tt00:4*i00p2-i00p1 ,float(c11^2*tt00)/*agrees prior bkmat[1,1]*/ ,i20p1:16*%pi^2*(integrate(integrate(r1*exp(-alis[1]*r1^2) *r2^4*exp(-alis[2]*r2^2),r1,r2,inf),r2,0,inf) +integrate(integrate(r1^2*exp(-alis[1]*r1^2)*r2^3 *exp(-alis[2]*r2^2),r2,r1,inf),r1,0,inf)) ,i20p2:16*%pi^2*integrate(r*exp(-alis[1]*r^2),r,0,inf) *integrate(r^4*exp(-alis[2]*r^2),r,0,inf) ,i20p3:16*%pi^2*integrate(r^2*exp(-alis[1]*r^2),r,0,inf) *integrate(r^3*exp(-alis[2]*r^2),r,0,inf) ,ldisp(["i00p1",float(i00p1)," i00p2",float(i00p2) ," i20p1",float(i20p1)," i20p2",float(i20p2)," i20p3" ,float(i20p3)]) ,tt20:-i20p1+2*(i20p2+i20p3) ,bkmat[2,1]:float(c11*(c21*tt00+c22*tt20)) ,ldisp(["bkmat",bkmat]) ,tt:float(c11*sum(-qq[i]*(c22*3*PEM(2,alis,[0,2, 0,0, 0,0],pmat[i]) +c21*PEM(2,alis,[0,0, 0,0, 0,0],pmat[i])) ,i,1,nps*(nps+1)/2)));

(%t27) [i00p1,24.73942945119314, i00p2,34.98683665524971, i20p1,30.92428681399141, i20p2,52.48025498287457, i20p3,34.98683665524971] (%t28) [bkmat,234.699591319584208.2772551062240710] Evaluation took 0.1900 seconds (0.1900 elapsed) (%o28) 8.277255106224101

(nps:2,c11:1.42729914,c21 :-1.20808673, c22:1.00673914 , alis:[1,1],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[0.5, 0.5 ] , ls:[0, 0],irps1lis:[[0,0]],irps2lis:[[0,0]] ,nm1lis:[ c11],nm2lis:[ c11] ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm); (nps:2,c11:1.42729914,c21 :-1.20808673, c22:1.00673914 , alis:[1,1],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[0.5, 0.5 ] , ls:[0, 0],irps1lis:[[0,0]],irps2lis:[[0,0],[0,1]] ,nm1lis:[ c11],nm2lis:[ c21,c22] ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm); (nps:2,c11:1.42729914,c21 :-1.20808673, c22:1.00673914 , alis:[1,1],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[0.5, 0.5 ] , ls:[0, 0],irps1lis:[[0,0],[0,1]],irps2lis:[[0,0],[0,1]] ,nm1lis:[ c21,c22], nm2lis:[ c21,c22] ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm); /*ck nps 2 00 scalar diag*/ (nps:2, alis:[1,1],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[0.500136137, 0.500136137 ] , ls:[3, 3],irps1lis:[[0,0]] ,irps2lis:irps1lis ,nm1lis:[7.33087808E-02] ,nm2lis:nm1lis ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm);

Evaluation took 0.0600 seconds (0.0600 elapsed) (%o3) 0.99999984836816

(nps:2, alis:[1,1],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[0.500136137, 0.500136137 ] , ls:[0, 0],irps1lis:[[0,0],[0,1]] ,irps2lis:[[0,0],[0,1]] ,nm1lis:[ -1.20794177,1.00664556] ,nm2lis:[ -1.20794177,1.00664556] ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm); (nps:2, alis:[1,1],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[0.500136137, 0.500136137 ] , ls:[0, 0],irps1lis:[[0,0],[0,1]] ,nm1lis:[ -1.20794177,1.00664556] ,irps2lis:[[0,0]] ,nm2lis:[ 1.42718267] ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm); /*finite nuclear mass helium*/ (nps:2,mm:matrix([234.654449,8.28337097],[8.28337097,158.217590]) ,eij:eigens_by_jacobi(mm),m2:part(eij,2) ,eival:part(eij,1),m1:mm.m2,ms:transpose(m2).m1 ,ldisp(["mm",mm,"eival",eival,"m2",m2,"m1",m1,"ms",ms]) ,float((lmax(eival)/(4*%pi)^nps)^2) ); /*20/Mar/2018 20:55 LEAVE diag's at 1 so jacobian for change of coord's is 1, that is NOT orthonormal but still orthog in KE space*/ (nps:4,amat:ident(nps),bmat:ident(nps), for i:2 thru nps do(ri:sum (m[k],k,1,i-1)+m[nps+1],for j thru i-1 do amat[i,j]:-m[j]/ri) ,ldisp(["amat",amat]), for i thru nps-1 do for k:i+1 thru nps do bmat[k,i]:-ratsimp(sum(amat[k,j]*bmat[j,i],j,i,k-1)) ,ldisp(["bmat",bmat,"amat.bmat",ratsimp(amat.bmat)]) ,for i thru nps do for j:i thru nps do(ro:factor(sum(amat[i,k] *amat[j,k]/m[k]/2,k,1,i) +sum(amat[i,k],k,1,i) *sum(amat[j,k],k,1,j)/m[nps+1]/2) ,ldisp([i,j,ro])) ,for i thru nps-1 do for j:i+1 thru nps do(tji:ratsimp(bmat[j]-bmat[i]) ,ldisp(["j",j,"i",i,"tji",tji,"tji.amat",ratsimp(tji.amat)])));

(%t4) [amat,1000-m1m5+m1100-m1m5+m2+m1-m2m5+m2+m110-m1m5+m3+m2+m1-m2m5+m3+m2+m1-m3m5+m3+m2+m11] (%t5) [bmat,1000m1m5+m1100m1m5+m1m2m5+m2+m110m1m5+m1m2m5+m2+m1m3m5+m3+m2+m11,amat.bmat,1000010000100001] (%t6) [1,1,m5+m12*m1*m5] (%t7) [1,2,0] (%t8) [1,3,0] (%t9) [1,4,0] (%t10) [2,2,m5+m2+m12*m2*

m5+m1

] (%t11) [2,3,0] (%t12) [2,4,0] (%t13) [3,3,m5+m3+m2+m12*m3*

m5+m2+m1

] (%t14) [3,4,0] (%t15) [4,4,m5+m4+m3+m2+m12*m4*

m5+m3+m2+m1

] (%t16) [j,2,i,1,tji,[-m5m5+m1,1,0,0],tji.amat,-1100] (%t17) [j,3,i,1,tji,[-m5m5+m1,m2m5+m2+m1,1,0],tji.amat,-1010] (%t18) [j,4,i,1,tji,[-m5m5+m1,m2m5+m2+m1,m3m5+m3+m2+m1,1],tji.amat,-1001] (%t19) [j,3,i,2,tji,[0,-m5+m1m5+m2+m1,1,0],tji.amat,0-110] (%t20) [j,4,i,2,tji,[0,-m5+m1m5+m2+m1,m3m5+m3+m2+m1,1],tji.amat,0-101] (%t21) [j,4,i,3,tji,[0,0,-m5+m2+m1m5+m3+m2+m1,1],tji.amat,00-11] Evaluation took 0.3400 seconds (0.3400 elapsed) (%o21) done

/*8/Mar/2018 5:07 more suitable for fortran let index 0 go to index nps+1 and further in terms of actual masses so mm[0]-> 1/m[nps+1]*/ (nps:4,amat:zeromatrix(nps,nps), amat[1,1]:1, for i:2 thru nps do(amat[i,i]:1 ,ri:sum(m[k],k,1,i-1)+m[nps+1], for j thru i-1 do amat[i,j]:-m[j]/ri) ,ldisp(["amat init",amat]) /* , for i thru nps do (ro:ratsimp(sum(amat[i,k]^2/m[k] ,k,1,i) +sum(amat[i,k],k,1,i)^2/m[nps+1]),for j thru i do amat[i,j]:amat[i,j] /sqrt(ro)), amat:amat*sqrt(2), ldisp(["amat",radcan(amat)])*/ ,for i thru nps do for j:i thru nps do(ro:ratsimp(sum(amat[i,k]*amat[j,k]/m[k]/2,k,1,i) +sum(amat[i,k],k,1,i) *sum(amat[j,k],k,1,j)/m[nps+1]/2) ,ldisp([i,j,ro])) ,disp("agrees fortran uselink... expression for rkm"));

(%t2) [amat init,1000-m1m5+m1100-m1m5+m2+m1-m2m5+m2+m110-m1m5+m3+m2+m1-m2m5+m3+m2+m1-m3m5+m3+m2+m11] (%t3) [1,1,m5+m12*m1*m5] (%t4) [1,2,0] (%t5) [1,3,0] (%t6) [1,4,0] (%t7) [2,2,m5+m2+m12*m2*m5+2*m1*m2] (%t8) [2,3,0] (%t9) [2,4,0] (%t10) [3,3,m5+m3+m2+m12*m3*m5+

2*m2+2*m1

*m3] (%t11) [3,4,0] (%t12) [4,4,m5+m4+m3+m2+m12*m4*m5+

2*m3+2*m2+2*m1

*m4] agrees fortran uselink... expression for rkm Evaluation took 0.0800 seconds (0.0800 elapsed) (%o12) done

 /*13/Mar/2018 12:12 found can simplify and directly have*/ disp("not using this orthonormalized") /* (nps:4,amat:zeromatrix(nps,nps),for i thru nps do(amat[i,i]:sqrt(2*m[i]*(m[nps+1] +sum(m[j],j,1,i-1))/(m[nps+1]+sum(m[j],j,1,i))),for j thru i-1 do amat[i,j]:-m[j] *sqrt(2*m[i]/(m[nps+1]+sum(m[j],j,1,i-1))/(m[nps+1]+sum(m[j],j,1,i))) ) , ldisp(["amat",amat]),for i thru nps do for j:i thru nps do(ro:radcan(sum(amat[i,k] *amat[j,k]/m[k]/2,k,1,i) +sum(amat[i,k],k,1,i)*sum(amat[j,k],k,1,j)/m[nps+1]/2) ,ldisp([i,j,ro]))) */;

not using this orthonormalized Evaluation took 0.0000 seconds (0.0000 elapsed) (%o56) done

(kill(nps),ratsimp((2*i*nps-2*nps-i**2+i-2)/2+nps+1));

Evaluation took 0.0000 seconds (0.0000 elapsed) (%o344) 2*i*nps-i2+i2

/*13/Mar/2018 13:43 direct expression for inverse of amat*/ (nps:4,kill(bmat),bmat:zeromatrix(nps,nps),for i thru nps do(bmat[i,i]:radcan(1/amat[i,i]) ,for k:i+1 thru nps do bmat[k,i]:-factor(radcan(sum(amat[k,j]*bmat[j,i],j,i,k-1)/amat[k,k]))) ,ldisp(["bmat",bmat,"amat.bmat",radcan(amat.bmat)]) ,for i thru nps-1 do for j:i+1 thru nps do(tji:radcan(bmat[j]-bmat[i]) ,ldisp(["j",j,"i",i,"tji",tji,"tji.amat",radcan(tji.amat)]))); (xlis:makelist(x[i]-x[nps+1],i,1,nps),ylis:amat.xlis ,for i thru nps do ldisp([i,ylis[i]]) ,for i thru nps do for j:i thru nps do( tt:ratsimp(sum(diff(ylis[i],x[k])*diff(ylis[j],x[k])/m[k]/2,k,1,nps+1)) ,ldisp([i,j,tt])));

(%t116) [1,[2*

x1-x5

m5+m1m1*m5]] (%t117) [2,[2*

x2-x5

m5+m2+m1m2*m5+m1*m2-2*m1*

x1-x5

m5+m1

*m5+m2+m1m2*m5+m1*m2]] (%t118) [3,[2*

x3-x5

m5+m3+m2+m1m3*m5+

m2+m1

*m3-2*m2*

x2-x5

m5+m2+m1

*m5+m3+m2+m1m3*m5+

m2+m1

*m3-2*m1*

x1-x5

m5+m2+m1

*m5+m3+m2+m1m3*m5+

m2+m1

*m3]] (%t119) [4,[2*

x4-x5

m5+m4+m3+m2+m1m4*m5+

m3+m2+m1

*m4-2*m3*

x3-x5

m5+m3+m2+m1

*m5+m4+m3+m2+m1m4*m5+

m3+m2+m1

*m4-2*m2*

x2-x5

m5+m3+m2+m1

*m5+m4+m3+m2+m1m4*m5+

m3+m2+m1

*m4-2*m1*

x1-x5

m5+m3+m2+m1

*m5+m4+m3+m2+m1m4*m5+

m3+m2+m1

*m4]] (%t120) [1,1,[1]] (%t121) [1,2,[0]] (%t122) [1,3,[0]] (%t123) [1,4,[0]] (%t124) [2,2,[1]] (%t125) [2,3,[0]] (%t126) [2,4,[0]] (%t127) [3,3,[1]] (%t128) [3,4,[0]] (%t129) [4,4,[1]] Evaluation took 0.1100 seconds (0.1100 elapsed) (%o129) done

(nps:4,xlis:makelist(x[i]-x[nps+1],i,1,nps),ylis:amat.xlis /*,for i thru nps do ldisp([i,ylis[i]])*/ ,for i thru nps do for j:i thru nps do( tt:ratsimp(sum(diff(ylis[i],x[k])*diff(ylis[j],x[k])/m[k]/2,k,1,nps+1)) ,ldisp([i,j,tt]))); /*invert lower triangular matrix direct like also easier to apply to fortran*/ (n:3,amat:zeromatrix(n,n),bmat:zeromatrix(n,n),for i thru n do for j thru i do amat[i,j]:a[i,j] ,ldisp(["amat",amat]), for i thru n do(bmat[i,i]:1/amat[i,i] ,for k:i+1 thru n do bmat[k,i]:-sum(amat[k,j]*bmat[j,i],j,i,k-1)/amat[k,k]) ,ldisp(["bmat",bmat,"amat.bmat",ratsimp(amat.bmat)]));

(%t19) [amat,a1,100a2,1a2,20a3,1a3,2a3,3] (%t20) [bmat,1a1,100-a2,1a1,1*a2,21a2,20a2,1*a3,2a1,1*a2,2-a3,1a1,1a3,3-a3,2a2,2*a3,31a3,3,amat.bmat,100010001] Evaluation took 0.0200 seconds (0.0200 elapsed) (%o20) [%t20]

(nps:4, alis:[4,10,8,1],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[7,11,5,1] ,ls:[0,0,0,0], irps1:[0,0,0,0], irps2:[0,0,0,0] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,ldisp(["tbra",tbra,"tket",tket]) ,tad:1, for i thru nps do tad:tad*integrate(r[i]^2*tbra*tket,r[i],0,inf) ,ldisp(["tad",tad]) );

(%t3) [tbra,%e-r42-8*r32-10*r22-4*r122,tket,%e-r42-8*r32-10*r22-4*r122] (%t4) [tad,%pi2*%e-3*r42-24*r32-30*r22-12*r122312*1032] Evaluation took 0.6400 seconds (0.6400 elapsed) (%o4) [%t4]

(nps:4, alis:[4,10,8,1],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[7,11,5,1] ,ls:[0,0,0,0], irps1:[0,0,0,0], irps2:[0,0,0,0] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,ldisp(["tbra",tbra,"tket",tket]) ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,ldisp(["tad",tad]) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i]*(diff(tket,r[i],2) +2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,ldisp(["tij",tij,"tad+tij",float(tad+tij)]) ,nm:23.4721661,ldisp(["(tad+tij)*nm^2",float((tad+tij)*nm^2)]) );

(%t2) [tbra,%e-r42-8*r32-10*r22-4*r122,tket,%e-r42-8*r32-10*r22-4*r122] (%t3) [tad,%pi22312*1032] (%t4) [tij,537*10*%pi225*2372,tad+tij,0.0018150708501672] (%t5) [(tad+tij)*nm^2,0.99999981966221] Evaluation took 0.6400 seconds (0.6400 elapsed) (%o5) [%t5]

(nps:4, alis:[4,10,8,1],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[7,11,5,1] ,ls:[0,0,0,0],irps1lis:[[0,0,0,0]],irps2lis:[[0,0,0,0]] ,nm1lis:[[23.4721661]] ,nm2lis:[[23.4721661]] ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm);

Evaluation took 0.5600 seconds (0.5600 elapsed) (%o2) [0.99999981966221]

(nps:4, alis:[4,10,8,1],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[7,11,5,1] ,ls:[0,0,0,0],irps1lis:[[0,0,0,0]],irps2lis:[[0,0,0,0],[0,0,0,1]] ,nm1lis:[[23.4721661]] ,nm2lis:[[-28.5353470],[19.0944176]] ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm);

Evaluation took 0.3000 seconds (0.3000 elapsed) (%o4) [1.0621665569310323*10-7]

(nps:4, alis:[4,10,8,1],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[7,11,5,1] ,ls:[0,0,0,0],irps1lis:[[0,0,0,0],[0,0,0,1]] ,irps2lis:[[0,0,0,0],[0,0,0,1]] ,nm1lis:[[-28.5353470],[19.0944176]] ,nm2lis:[[-28.5353470],[19.0944176]] ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm);

(%o1) [0.99999962054837]

(nps:4, alis:[4,10,8,1],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[7,11,5,1] , ls:[1, 0, 1, 2],irps1lis:[[0,0,0,0],[0,0,0,1],[0,0,1,0]] ,irps2lis:[[0,0,0,0],[0,0,0,1],[0,0,1,0]] ,nm1lis:[-51.9086838,2.02681106E-02,188.035919] ,nm2lis:[-51.9086838,2.02681106E-02,188.035919] ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm);

(%o2) 0.99999944872644

(nps:4, alis:[4,10,8,1],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[7,11,5,1] , ls:[1, 0, 1, 2],irps1lis:[[0,0,0,0],[0,0,0,1],[0,0,1,0]] ,irps2lis:[[0,0,0,0],[0,0,0,1],[0,0,1,0],[0,1,0,0]] ,nm1lis:[-51.9086838,2.02681106E-02,188.035919] ,nm2lis:[-30.6686954,3.99442762E-02,10.4058046,272.789520] ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,tij: ratsimp(sum(rkm[i]*ls[i]*(ls[i]+1)/r[i]^2,i,1,nps)*tket -sum(rkm[i] *(diff(tket,r[i],2)+2*diff(tket,r[i])/r[i]), i,1,nps)) ,tij:tbra*tij, for i thru nps do tij:integrate(r[i]^2*tij,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*(tad+tij))),nm);

Evaluation took 2.5000 seconds (2.5000 elapsed) (%o5) 2.5731472397261257*10-8

/* 25/Jun/2018 1:00 Next for only normalized wrt tad as in useKEplusPE.F03 */ (nps:3, alis:[1.44,.83,.83],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[5.44617767E-04, 0.500136137, 0.500136137] , ls:[0,0,0],irps1lis:[[0,0,0]],irps2lis:[[0,0,0]] ,nm1lis:[3.36989617],nm2lis:nm1lis ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*tad)),nm); (nps:3, alis:[1.44,.83,.83],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[5.44617767E-04, 0.500136137, 0.500136137] , ls:[0,0,0],irps1lis:[[0,0,0]],irps2lis:[[0,0,0],[0,0,1]] ,nm1lis:[3.36989617],nm2lis:[-4.12726307, 2.28375244] ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*tad)),nm); (nps:3, alis:[1.44,.83,.83],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[5.44617767E-04, 0.500136137, 0.500136137] , ls:[2,2,0],irps1lis:[[0,0,0],[0,0,1]],irps2lis:[[0,0,0],[0,0,1]] ,nm1lis:[-1.31544173, .727877736] ,nm2lis:[-1.31544173, .727877736] , nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*tad)),nm); /* 25/Jun/2018 13:53 tad only */ (nps:4, alis:[.35,11,.8,.8],er:exp(-sum(alis[i]*r[i]^2,i,1,nps)/2) ,rkm:[5.44617767E-04,4.08463355E-04,.500090778,.500090778] , ls:[2,2,0,0] ,irps1lis:[[0,0,0,0],[0,0,0,1],[0,0,1,0],[0,1,0,0],[1,0,0,0]] ,irps2lis:[[0,0,0,0],[0,0,0,1],[0,0,1,0],[0,1,0,0]] ,nm1lis:[-19.2381763,-9.50228923E-07 ,-9.50228923E-07,0,1.92381811] ,nm2lis:[-19.2381897,-9.33261106E-07, 0,60.462883] ,nm:0, for i1 thru length(irps1lis) do for i2 thru length(irps2lis) do ( irps1:irps1lis[i1], irps2:irps2lis[i2] ,tbra:prod(r[i]^(ls[i]+2*irps1[i]),i,1,nps)*er ,tket:prod(r[i]^(ls[i]+2*irps2[i]),i,1,nps)*er ,tad:tbra*tket, for i thru nps do tad:integrate(r[i]^2*tad,r[i],0,inf) ,nm:float(nm+nm1lis[i1]*nm2lis[i2]*tad)),nm);

[zmth]

don't know whether the lines correspond or not but on line 828 of this content.xml :



actually in my editor the funny thing here in the middle shows as a small left pointing arrow

[zmth]

supposed to be on line 828



[zmth]

keeps changing line 828 to nonsense when i press comment above. Try again

(small left arrow)

[gunterkoenigsmann]

As I have said I will have a working computer again in a week. I'll try to find a way to fix the file using my cellphone but I can't guarantee I find one.

But since you are still using the 2013 version of wxMaxima I am confident that a quick fix would be uninstalling wxMaxima and installing the version from https://wolfgang.dautermann.at/maxima/nightlybuild/ will resolve the problem: the current wxMaxima version will automatically fix all kinds of unopenable files I have ever seen from the 2013 version of wxMaxima.

Am 25. August 2018 11:04:23 MESZ schrieb zmth notifications@github.com:

thanks for your reply: wxMaxima version: 13.4.0 Maxima version: 5.31.1 Maxima build date: 2013-09-24 09:49:12 Host type: i686-pc-mingw32 Lisp implementation type: GNU Common Lisp (GCL) Lisp implementation version: GCL 2.6.8

before years ago same thing hapenned and they fixed it easily. I don't have an xml editor don't even know what XML means. Don't have any way to send it to peterpall . de as is not a valid email address. i renamed it binom.zip and will try to upload here binom.zip thanks

-- You are receiving this because you commented. Reply to this email directly or view it on GitHub: https://github.com/andrejv/wxmaxima/issues/1026#issuecomment-415954985

-- Diese Nachricht wurde von meinem Android-Gerät mit K-9 Mail gesendet.

[gunterkoenigsmann]

You are correct: the character is an Escape sign and not allowed in XML so it will cause the error you have encountered.

New enough versions of wxMaxima will automatically repair this error when opening the file.

[zmth]

So now have the .xml file and next if i 'fiddle' with this flle , deleting, changing things etc. then how do i get it back to a .wxmx file to try and load again ? On Saturday, August 25, 2018, 5:37:22 AM CDT, Gunter Königsman notifications@github.com wrote:

You are correct: the character is an Escape sign and not allowed in XML so it will cause the error you have encountered.

New enough versions of wxMaxima will automatically repair this error when opening the file.

— You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub, or mute the thread.

[zmth]

Uninstalling the older version and reinstalling newer. What all does that do ? and what directory does it use? Now can't even find the directories and tree stucture of the old version. Anyway so if i uninstall all prior is gone forever ? What if install the newer without uninstalling ?

[gunterkoenigsmann]

By default it installs itself into the directory c:/maxima_current

The rationale behind this is that this path doesn't contain a whitespace which the default install path in many countries does.

[gunterkoenigsmann]

If you want to transform content.xml back to a valid .wxmx file just create a .zip file containing the content.xml and rename it into something ending in .wxmx

The wxMaxima from dautermann.at is able to open content.xml even if it isn't inside a .zip file.

[gunterkoenigsmann]

My computer works again => opened the file with a new wxMaxima and saved it again. It should be openable now again by any wxMaxima now that was released after 2011. binom.zip

If you have still problems just comment on this thread and I will answer your messages.

[zmth]

ok thanks I did find on another computer did have a newer version 5.38.1 which did bring up the prior damaged file, with a notice, and corrected it also. If in newer version of maxima 5.38.1 bring up a filed saved with older version maxima 5.28.0-2 get a warning notice and also if bring up a file saved with newer 5.38.1 in the older version 5.28.0-2.

Funny seems the older version more efficient because the file saved by the newer version was around 87kb and after bringing it up in older version 5.28.0-2 and saving it then it was only about 13kb.

[gunterkoenigsmann]

There are 3 reasons:

• Newer versions of wxMaxima include higher-resolution plots in the worksheet by default.
• Newer versions of wxMaxima use antialiassing: If only a third of a pixel is hit by a line the pixel is only drawn light gray, not completely black or completely white. Which greatly improves the quality of the images, but also greatly increases the number of colors the image contains => we no more have to store 1..3 bit per pixel but 8...24 bit per pixel.

Additionally the text portion of the file is no more compressed:

The .wxmx file internally is a .zip archive, as is a .docx, a .odt or nearly every other file a modern program creates.

.zip files offer the feature of compressing turning of compression for their contents, which, as you have experienced, makes the file smaller, but comes with a number of drawbacks:

• if the text is compressed and the file is truncated somehow (disk or filesystem error or disk full or...) the information on how to uncompress the file is lost => reconstructing the contents of faulty compressed files is hard
  • many folks use a version control system like subversion or git that records all changes you ever have made to a file: As each individual change is small and the changes are recorded in a compressed form recording the change to uncompressed files doesn't use much memory (normally the memory used for all the changes you did in about 5 years is roughly as big as the actual file). But it provides you with an unlimited undo facility. Additionally if you store your git or subversion repository not on your local disk but on the other side of the world you'll have a complete backup of all the history of all of your files that will even survive a meteorite crash. If you try to do the same to a compressed file any change in the file is likely to change all of the compressed text => storing the history of compressed files tends to use up much memory.
  • and the amount of hard disk space saved often is small as the text rarely is bigger than a few kilobytes and therefore can only be reduced by a fraction of a few kilobytes.

=> Most of the increase of the file size is caused by an improved image quality. The rest (which should be small) helps with recovering broken files or saving the file's history.