#include
#include
#include
#define MAX 20
float equation(float x,float y)
{
return(1+(y*y));
}
void main()
{
FILE *fp;
int i=0,count=-1;
float lower,upper,h,y1,xvalue[MAX],yvalue[MAX],result;
float function[MAX],search,final,temp;
fp=fopen("milne.dat","w"
clrscr();
printf("MILNE'S METHOD "
fprintf(fp,"MILNE'S METHOD "
printf("\n"
fprintf(fp,"\n"
printf("\nEnter the lower bound of x : "
fprintf(fp,"\nEnter the lower bound of x : "
scanf("%f",&lower);
fprintf(fp,"%f",lower);
printf("\nEnter the upper bound of x : "
fprintf(fp,"\nEnter the upper bound of x : "
scanf("%f",&upper);
fprintf(fp,"%f",upper);
printf("\nEnter the value of y(lower) : "
fprintf(fp,"\nEnter the value of y(lower) : "
scanf("%f",&y1);
fprintf(fp,"%f",y1);
printf("\nEnter the value of h : "
fprintf(fp,"\nEnter the value of h : "
scanf("%f",&h);
fprintf(fp,"%f",h);
printf("\nEnter the value of x for which you want to find y :"
fprintf(fp,"\nEnter the value of x for which you want to find y :"
scanf("%f",&search);
fprintf(fp,"%f",search);
xvalue[i]=lower;
yvalue[i]=y1;
for(i=0;xvalue[i]<=upper;i++)
{
xvalue[i+1]=xvalue[i]+h;
}
for(i=0;xvalue[i]<=upper;i++)
{
result=equation(xvalue[i],yvalue[i]);
yvalue[i+1]=yvalue[i]+(h*result);
}
printf("\n\n"
fprintf(fp,"\n\n"
printf("The table is "
fprintf(fp,"The table is "
printf("\n\n"
fprintf(fp,"\n\n"
printf(" i x y f(x,y) "
fprintf(fp," i x y f(x,y) "
printf("\n\n"
fprintf(fp,"\n\n"
for(i=0;xvalue[i]<=upper;i++)
{
function[i]=equation(xvalue[i],yvalue[i]);
printf(" %d. %.4f %.4f %.4f ",i,xvalue[i],yvalue[i],function[i]);
fprintf(fp," %d. %.4f %.4f %.4f ",i,xvalue[i],yvalue[i],function[i]);
count=count+1;
printf("\n"
fprintf(fp,"\n"
}
yvalue[search]=yvalue[count-3]+(4*h/3)*((2*function[count-2])-(function[count-1])+(2*function[count]));
final=equation(search,yvalue[search]);
yvalue[search]=yvalue[count-1]+(h/3)*((function[count-1])+(4*function[count])+final);
printf("\n\n"
fprintf(fp,"\n\n"
printf("Approximate value is : %.4f ",yvalue[search]);
fprintf(fp,"Approximate value is : %.4f ",yvalue[search]);
fclose(fp);
getch();
}
1>------ Построение начато: проект: ad, Конфигурация: Debug Win32 ------
1> ad.cpp
1>c:\users\eng\documents\visual studio 2010\projects\ad\ad\ad.cpp(19): warning C4996: 'fopen': This function or variable may be unsafe. Consider using fopen_s instead. To disable deprecation, use _CRT_SECURE_NO_WARNINGS. See online help for details.
1> c:\program files\microsoft visual studio 10.0\vc\include\stdio.h(234): см. объявление "fopen"
1>c:\users\eng\documents\visual studio 2010\projects\ad\ad\ad.cpp(20): error C3861: clrscr: идентификатор не найден
1>c:\users\eng\documents\visual studio 2010\projects\ad\ad\ad.cpp(75): error C2108: индекс имеет нецелый тип
1>c:\users\eng\documents\visual studio 2010\projects\ad\ad\ad.cpp(76): error C2108: индекс имеет нецелый тип
1>c:\users\eng\documents\visual studio 2010\projects\ad\ad\ad.cpp(77): error C2108: индекс имеет нецелый тип
1>c:\users\eng\documents\visual studio 2010\projects\ad\ad\ad.cpp(80): error C2108: индекс имеет нецелый тип
1>c:\users\eng\documents\visual studio 2010\projects\ad\ad\ad.cpp(81): error C2108: индекс имеет нецелый тип
========== Построение: успешно: 0, с ошибками: 1, без изменений: 0, пропущено: 0 ==========
vt
Исправлено модератором.
Уважайте читателей сообщества - заключайте код в соответствующий тег!