Сумма ряда и факториал

[rw61]

Сообщение от Аватар

Как насчет начального значения s? Да и k в формуле суммы от 0, а не от 1

первичная сумма равна члену ряда.
при k=0 прога выдаёт ошибку.

решил я вопрос. косяк был в формуле и в прописке факториала

итоговый код, если вдруг, кому будет интересно, выкладываю.

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

#include <vcl.h>
#pragma hdrstop
#include <math>
#include "Primer3.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
: TForm(Owner)
{

}
//---------------------------------------------------------------------------
void __fastcall TForm1::FormCreate(TObject *Sender)
{
Edit1->Text="0,1";
Edit2->Text="1,0";
Edit3->Text="30";
Edit4->Text="0,1";
Memo1->Lines->Add("Лабораторная работа 3");
}

int fact (int num)
{
int result = 1, i;
for (i = 2; i <= num; i++)
result *= i;
return result;
}
//---------------------------------------------------------------------------
void __fastcall TForm1::Button1Click(TObject *Sender)
{ double a = 0.1, b = 1.0, h = 0.1, e = exp(1), y, s, r, t = cos(M_PI / 4), v = sin(M_PI / 4), m;
int k, n = 30;
for (double x = a; x <= b; x += h)
{
y = pow(e, x * t) * cos(x * v);
k = 0;
//Первый член ряда, при k = 0, равен 1, можно сразу присвоить 1 и не нужна функция fact
m = cos(k * M_PI / 4) * pow(x ,k) / fact(k);
s = m;
for (k = 1; k <= n; k++)
{
// Вычисление члена ряда,в зависимости от k
m *= cos(k * M_PI / 4) / cos((k - 1) * M_PI / 4) * x / k;
s += m; // Суммируем
}
r = fabs(y - s);
Memo1->Lines->Add("При x= "+FloatToStrF(x,ffFixed,8,2)
+"\t Y(x)= "+FloatToStrF(y,ffFixed,8,5)
+"\t S(x)= "+FloatToStrF(s,ffFixed,8,5)
+"\t |Y(x)-S(x)|= "+FloatToStrF(r,ffFixed,8,5));
}
}
void __fastcall TForm1::Button2Click(TObject *Sender)
{
Memo1->Clear();

}