Итак, я пытаюсь построить вывод этой функции интеграции Эйлера:
typedef double F(double,double);
using std::vector;
void euler(F f, double y0, double a, double b, double h,vector<POINT> Points)
{
POINT Pt;
double y_n = y0;
double t = a;
for (double t = a; t != b; t += h )
{
y_n += h * f(t, y_n);
Pt.x = t; // assign the x value of the point to t.
Pt.y = y_n; // assign the y value of the point to y_n.
Points.push_back(Pt);
}
}
// Example: Newton's cooling law
double newtonCoolingLaw(double, double t)
{
return t; // return statement ends the function; here, it gives the time derivative y' = -0.07 * (t - 20)
}
Я пытаюсь использовать функцию Polyline() в приложении Win32, поэтому я делаю это в случае WM_PAINT
:
case WM_PAINT:
{
hdc = BeginPaint(hWnd, &ps);
//Draw lines to screen.
hPen = CreatePen(PS_SOLID, 1, RGB(255, 25, 5));
SelectObject(hdc, hPen);
using std::vector;
vector<POINT> Points(0);
euler(newtonCoolingLaw, 1, 0, 20, 1,Points);
POINT tmp = Points.at(0);
const POINT* elementPoints[1] = { &tmp };
int numberpoints = (int) Points.size() - 1 ;
Polyline(hdc,elementPoints[1],numberpoints);
Когда я перенаправляю свой ввод-вывод на консоль, вот выходные данные для переменных:
Я могу нарисовать ожидаемые линии на экране с помощью MovetoEx(hdc,0,0,NULL)
и LineTo(hdc,20,20)
, но по какой-то причине ни одна из этих функций не работает с моим vector<POINT> Points
. Какие-либо предложения?