Сформировать массив. Ошибка: Invalid floating point operation. Не могу понять в чем дело

[Sv.]

Дана функция y=3sin^2(x)+cos(6x+1)
Сформировать массив значений аргумента и функции в точках экстремумов (минимумы и максимумы): x1*, y1*, x2*, y2 *, … .

//---------------------------------------------------------------------------

#include <vcl.h>
#include <conio.h>
#include <stdio.h>
#include <math.h>
#pragma hdrstop

//---------------------------------------------------------------------------

#pragma argsused
int massiv_y(double x1,double x2, double h);
void main()
{ double x1, x2, h;
printf("Vvedite interval i shag funkcii: \n");
scanf("%lf%lf%lf", &x1, &x2, &h);
massiv_y(x1,x2,h);
getch();
}

int massiv_y(double x1,double x2, double h)
{
int i=0;
double y[100], y_min, x_min, y_max, x_max;
for(x1;x1<=x2;x1+=h)
{y[i++]=3*sin(x1)*sin(x1)+cos(6*x1+1);
if(y[i]<y[i+1])
{y_min=y[i];
x_min=x1;
printf("y_min=%lf x_min=%lf", y_min, x_min);} //Здесь ошибка
else
{y_max=y[i];
x_max=x1;
printf("y_max=%lf x_max=%lf", y_max, x_max);}}
return 0;
}

[Stilet]

На какой строчке?

[Sv.]

На какой строчке?

где комментарий стоит

с этой ошибкой разобралась. увеличила размер массива.
но выводит все x, но не при всех x считает значение y


значит проблема в условии?

и еще вот такую строку вывода сделала:

printf("\n y_min=%.3lf  x_min=%lf", y_min, x_min);

но выводит много цифр

[Abstraction]

//Пусть i=3
y[i++]=3*sin(x1)*sin(x1)+cos(6*x1+1);//Записали что-то в y[3], i=4
if(y[i]<y[i+1]) //Если y[4]<y[5] ...

Короче: осторожней с инкрементом.

[Sv.]

спасибо

[Sv.]

//---------------------------------------------------------------------------

#include <vcl.h>
#include <conio.h>
#include <stdio.h>
#include <math.h>
#pragma hdrstop

//---------------------------------------------------------------------------

#pragma argsused
int massiv_y(double x1,double x2, double h);
void main()
{ double x1, x2, h;
printf("Vvedite interval i shag funkcii: \n");
scanf("%lf%lf%lf", &x1, &x2, &h);
massiv_y(x1,x2,h);
getch();
}

int massiv_y(double x1,double x2, double h)
{
int i=0;
double y[400], y_min, x_min, y_max, x_max;
for(x1;x1<x2;x1+=h)
{y[++i]=3*sin(x1)*sin(x1)+cos(6*x1+1);
if(y[i-1]<y[i] && y[i]>y[i+1])
{y_max=y[i];
x_max=x1;
printf("\n x_max=%.3lf  y_max=%lf", x_max, y_max);}
else if(y[i-1]>y[i] && y[i]<y[i+1])
{y_min=y[i];
x_min=x1;
printf("\n x_min=%.3lf  y_min=%lf", x_min, y_min);}}
return 0;
}

почему при изменении размера массива, вычисляются другие значения:


и не все точки совпадают с графиком, хотя условие вроде правильное.

[Abstraction]

Дубль два:

//Пусть i=3
y[++i]=3*sin(x1)*sin(x1)+cos(6*x1+1);//В y[4] что-то записано, i=4
if(y[i-1]<y[i] && y[i]>y[i+1]) // Если y[3]<y[4] && y[4]>y[5]...
//Только вот засада в том, что в y[5] ещё никто ничего не писал, там мусор

Лучше сначала заполнить массив y. И настоятельно не советую менять значения x1, x2 и h: это может открыть двери достаточно неприятным ошибкам:

double y[400], y_min, x_min, y_max, x_max;
i=0;
for(double x = x1; x<x2; x+=h, ++i){ //Неплохо бы вначале удостовериться, что x1+400h>=x2
  y[i] = 3*sin(x)*sin(x)+cos(6*x+1); //Сначала заполнили y
}
i=0; //Вернули i "в начало"
for(double x = x1; x<x2; x+=h, ++i){
  //Пошёл основной цикл. Следует быть внимательными, чтобы ненароком не обратиться к y[-1] или y[i+1] при "последнем" i

[Sv.]

будет ли вот так правильно?

//---------------------------------------------------------------------------

#include <vcl.h>
#include <conio.h>
#include <stdio.h>
#include <math.h>
#pragma hdrstop

//---------------------------------------------------------------------------

#pragma argsused
int massiv_y(double x1,double x2, double h);
void main()
{ double x1, x2, h;
printf("Vvedite interval i shag funkcii: \n");
scanf("%lf%lf%lf", &x1, &x2, &h);
massiv_y(x1,x2,h);
getch();
}

int massiv_y(double x1,double x2, double h)
{
int i;
double y[400], y_min, x_min, y_max, x_max;
i=0;
for(double x=x1; x<=x2; x+=h, ++i)
{y[i]=3*sin(x)*sin(x)+cos(6*x+1);}
i=0;
for(double x=x1; x<x2; x+=h, ++i)
{if(x==x1)
{continue;}
y[i]=3*sin(x)*sin(x)+cos(6*x+1);
if(y[i-1]<y[i] && y[i]>y[i+1])
{y_max=y[i];
x_max=x;
printf("\n x_max=%.3lf  y_max=%lf", x_max, y_max);}
else if(y[i-1]>y[i] && y[i]<y[i+1])
{y_min=y[i];
x_min=x;
printf("\n x_min=%.3lf  y_min=%lf", x_min, y_min);}}
return 0;
}

[Abstraction]

Если гарантируется, что в x1-x2 укладывается целое число h, то вроде нормально; чуть более обще будет условие выхода во втором цикле поменять с x<x2 на x+h<=x2.

[Sv.]

Abstraction, спасибо большое))