{X_2,Y_2} = panel_234(curves_Ramsay2, curves_KR2, XRamsay, xk2, 2, 0.5);
// curve_Ramsay2 has 3 cols, for 95,99 and 999
// curve_KR2 has also 3 cols.
{X_3,Y_3} = panel_234(curves_Ramsay3, curves_KR3, XRamsay, xk3, 4.5, 0.5);
// curve_Ramsay3 has 3 cols, for 95,99 and 999
// curve_KR3 has also 3 cols.
{X_4,Y_4} = panel_234(curves_Ramsay4, curves_KR4, XRamsay, xk4, 6.5, 0.5);
// curve_Ramsay4 has 3 cols, for 95,99 and 999
// curve_KR4 has also 3 cols.

Line Zhao1:

new_dir = ChangeDir("E:\\Timo1\\GAUSS_Figure");

Line Zhao2:
_ptek="Egarch_n_Isoquants_cloud.tkf";

begwind; window(2,2,0); setwind(1);
_plwidth = 7.5;
_psymsiz ={ 6 5 5 5}; // symbol size
_plctrl = {-1 0 0 0}; // 0 for line, -1 for disconnected points
_pltype = { 6 6 6 6};
_pstype = { 4 1 1 1}; // 4 for a "plus", 1 for a circle
title("(a)");
xy(X_1,Y_1);
// 1st col of X_1 is the empirical values
// 2nd col of X_1 is the simulated value
nextwind;
title("(b)");
xy(X_2,Y_2);
nextwind;
title("(c)");
xy(X_3,Y_3);
nextwind;
title("(d)");
xy(X_4,Y_4);
endwind;

Line Zhao3:

default_dir=cdir("C:\\gauss6.0"); "Save and Click"; waitc;

standard way is given above via the three key lines below"Line Zhaoj= where j=1,2,3

try it. Egarch_n_Isoquants_cloud.tkf can be then opened easiely by double clicking and you can also convert it to another format after open it.

Hope this useful!