Задача на обтимізацію.

Складіть програму мовою С++ для знаходження функції y(x), де x Є [a, b] з кроком h. Дані для y(x), a, b, h вказані в таблиці.

y = x + sin(3x);

[a, b], h   [-1, 3], h = 0.3

Програма:

С++:

#include <iostream>

#include <cmath>

using namespace std;

int main()

{

   float y, a, b, h;

   cout << "a = ";  cin >> a;

   cout << "b = ";  cin >> b;

   cout << "h = ";  cin >> h;

   while (a <= b) {

       y = a + sin(3 * a);

       cout << "y = " << y << " при x = " << a << endl;

       a += h;

   }

}

Результат:

a = -1

b = 3

h = 0.3

-1.14112 при x = -1

-1.56321 при x = -0.7

-1.33204 при x = -0.4

-0.39552 при x = -0.1

0.764643 при x = 0.2

1.4975 при x = 0.5

1.47546 при x = 0.8

0.942254 при x = 1.1

0.528424 при x = 1.4

0.774186 при x = 1.7

1.72058 при x = 2

2.87844 при x = 2.3

3.59854 при x = 2.6

3.56297 при x = 2.9

Важное замечание: Для типа long long int требуется от компилятора поддержка стандарта C99.

Приложил файл расширения .cpp, чтобы обойти ограничение сайта. При работе с кодом измените на .c

#include <stdio.h>

long long int NOD(long long int var1, long long int var2) //Алгоритм Евклида

{

   while(var1 != var2)

   {

       if(var1 > var2)

           var1 = var1 - var2;

       else var2 = var2 - var1;

   }

   return var1;

}

long long int NOK(long long int var1, long long int var2)

{

   return (var1 * var2) / NOD(var1, var2);

}

int main()

{

   long long int N, K;

   scanf("%lli%lli", &N, &K);

   printf("Количество оборотов: %lli\n", NOK(N, K));

}