Добавил:
Maestro
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:LAB8
.H#include <stdio.h>
#include <math.h>
#include <conio.h>
void main()
{
//FILE *f1;
int i=0;
float x,x0=1,x2=2,n=2000,y0=-1;
float ye,y1e;
float k1,k2,k3,k4;
float y;
float y2,y1;
float c;
float h=(x2-x0)/n;
x=x0;
ye=y1=y0;
c=-exp(-3.0/2.0);
//f1=fopen("lab67.txt","r+");
clrscr();
do
{
y1e=ye+h*((x+1)*ye);
y=c*exp((x*x)/2+x);
k1=h*((x+1)*y1);
k2=h*((x+1+h/2.0)*(y1+k1/2.0));
k3=h*((x+1+h/2.0)*(y1+k2/2.0));
k4=h*((x+1+h)*(y1+k3));
y2=y1+((1.0/6.0)*(k1+2*k2+2*k3+k4));
printf("%4d %1.2f %5.2f %5.3f %9.5f %5.3f %12.8f\n",i++,x,y,ye,fabs(y-ye),y1,fabs(y-y1));
//fprintf(f1,"%2d %5.2f %5.1f %5.1f %7.3f %7.3f %7.3f %7.3f %7.3f %7.3f %7.3f %7.3f %7.3f\n",n++,x,t,t1,f,l2,fabs(f-l2),N1,fabs(f-N1),N2,fabs(f-N2),P,fabs(f-P));
ye=y1e;
y1=y2;
x+=h;
}
while(x<x2);
getch();
}
Соседние файлы в предмете Вычислительная математика