"--------------------------------------------"
" How to make lin-log10 paper. (intro)       "
"(log in Maxima,base  '%e'=2.718281828459045)"
"--------------------------------------------"
peter.vlasschaert@gmail.com,17/01/2017

we need conversion factor:
maxima: log(x) =  y <=> x = (exp(1))^y

kill(all)$

How to find : e

float((%e)^1);

float(exp(1));

float(log(9))/float(log(10));

principle 'for' => log ('scale' part)

don't use (....)$

for i:1 thru 6 step 1 do
 ( conv:float(1/log(10)),
   log10:float(log(i)*conv),
   print("number = ",i,",","log10(",i,") = ",log10)
 )     ;

fill a array:without know the dimension

use of: makelist

-----------------------------------------------------------

v:makelist(i^2,i,1,10);
v[3];

----------------------------------------------------------

principle 'makelist' => log ('scale' part)

1e) definition

w(j):=float(log(j)*conv);

check : definition

w(2);

2e) make list,vv

vv:makelist(w(j),j,1,10);

some manupulations on the list command

vv[3];

ww:append(vv,[w(3)]);

zz:delete(w(3),ww);

length(zz);

write_data(zz,"zzfilelist.dat")$

printfile("zzfilelist.dat")$

zz:append([zz],[1,2,3,4,5,6,7,8,9]);

write_data(zz,"zzfilearray.dat")$

printfile("zzfilearray.dat")$

data:read_matrix("zzfilearray.dat");

element :from matrix 'data'

data[1,3];

row : 'first row' => data

data[1];

row : 'second row' => data

data[2];

datav:transpose(data);

data[1,1];

plot2d ' discrete data'

x_cor:[1,2,3,4,5]$
y_cor:[0.2,0.5,0.7,0.9,1.2]$
xy:[[1,0.2],[2,0.5],[3,0.7],[4,.9],[5,1.2]]$

plot2d([discrete,x_cor,y_cor],[style,points])$

data plot from above from the above ('makelist','for')

plot2d([discrete,data[2],data[1]],[style,points],[legend,"data"],[xlabel,"data[2]"],
    [ylabel,"data[1]"])$