Now my plot works, thanks for help:
LG(x,p):=Round(log(x)*p)/p;
Round(v):=floor(v+0.5);
Td(t,h):=-((1188500*t+282506450)*log(h)+20301957*t)/((5000*t+1188500)*log(h)-20301957);
TD(t,h,p):=-Round(((1188500*t+282506450)*LG(h,p)+20301957*t)/(Round(((5000*t+1188500)*LG(h,p)-20301957)/P)));
P:1000;
E(t,h,p):=(Td(t,h)*P-TD(t,h,p))/P;
plot3d(E(t,h,2^15), [t,-30,50], [h,0.005,1], [plot_format,gnuplot]);
Next I would like to compute min and max on some 2D interval
instead of plot3d.
How to do that?
Thanks,
Albrecht