34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
|
{
Ly.push_back(r);
}
}
// write header
out << "set data style lines" << endl
<< "set xlabel \"time\"" << endl
<< "set grid" << endl
<< endl;
// write states (X11)
out << "set term X11" << endl
<< "plot\\" << endl;
for (list<record_t>::iterator i = Lx.begin(); i != Lx.end() ; i++)
{
if (Lx.begin() != i)
{
out << ",\\" << endl;
}
out << "\t\"MTT_work/" << sys << "_odes.dat2\" using 1:" << 2 + Ly.size() + i->num
<< " axes x1y1"
<< " title \"" << i->mod << "_" << i->cmp << "\"";
}
out << endl << endl << "pause(-1)" << endl << endl;
// write states (postscript)
out << "set term postscript eps color" << endl
<< "set output \"" << sys << "_states.eps\"" << endl
<< "replot" << endl << endl;
// write outputs (X11)
out << "set term X11" << endl
<< "plot\\" << endl;
for (list<record_t>::iterator i = Ly.begin(); i != Ly.end() ; i++)
{
if (Ly.begin() != i)
{
out << ",\\" << endl;
}
out << "\t\"MTT_work/" << sys << "_odes.dat2\" using 1:" << 1 + i->num
<< " axes x1y1"
<< " title \"" << i->mod << "_" << i->cmp << "\"";
}
out << endl << endl << "pause(-1)" << endl << endl;
// write outputs (postscript)
out << "set term postscript eps color" << endl
<< "set output \"" << sys << "_outputs.eps\"" << endl
<< "replot" << endl << endl;
return out;
}
#ifdef MAIN
int main(int argc, char *argv[])
{
|
>
|
>
<
<
<
<
<
<
|
<
|
<
<
|
<
|
<
<
<
<
<
<
<
<
<
|
<
|
<
<
|
<
|
<
<
<
|
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
|
{
Ly.push_back(r);
}
}
// write header
out << "wait=-1" << endl
<< "set data style lines" << endl
<< "set xlabel \"time\"" << endl
<< "set grid" << endl
<< "set term X11" << endl
<< endl;
// write states (X11)
for (list<record_t>::iterator i = Lx.begin(); i != Lx.end() ; i++)
{
out << "plot " << "\"MTT_work/" << sys << "_odes.dat2\" using 1:" << 2 + Ly.size() + i->num
<< " axes x1y1 title \"" << i->mod << "_" << i->cmp << ";" << endl
<< "pause(wait);" << endl;
}
// write outputs (X11)
for (list<record_t>::iterator i = Ly.begin(); i != Ly.end() ; i++)
{
out << "plot " << "\"MTT_work/" << sys << "_odes.dat2\" using 1:" << 1 + i->num
<< " axes x1y1 title \"" << i->mod << "_" << i->cmp << endl
<< "; pause(wait);" << endl;
}
return out;
}
#ifdef MAIN
int main(int argc, char *argv[])
{
|