метод милна

[Bun4ecTep]

код тут

#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 ==========

меняю fopen na fopen_s , пишет что функция не принимает два аргумента

vt

Исправлено модератором.
Уважайте читателей сообщества - заключайте код в соответствующий тег!

@темы: C++

Comments

[Караидель]

msdn.microsoft.com/ru-ru/library/z5hh6ee9.aspx

[alba-longa]

Пользоваться тегами для оформления кода религия не позволяет? Переделайте, так код абсолютно нечитабелен.
Что сверху в инклюдах подключается?
Переменная search какой тип имеет? А каким может быть тип индекса массива?