수치해석레포트3
본 자료는 미리보기를 지원하지 않습니다.
닫기
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
해당 자료는 2페이지 까지만 미리보기를 제공합니다.
2페이지 이후부터 다운로드 후 확인할 수 있습니다.

소개글

수치해석레포트3에 대한 보고서 자료입니다.

목차

3-1 largrange interplation method

3-2 newton iterpolation method

본문내용

#include
#include
#include
#define fa(x) (exp(x))
int k,m,n1;
double s[10],x[10],f[10],p[10];


void lagrange(void);

void main(void)
{

int i;
double error;

/* 파일 입출력을 위한 파일 지정 부분 */
FILE *in, *out;
in=fopen(\\\"data.in\\\",\\\"r\\\");
out=fopen(\\\"data.out\\\", \\\"w\\\");

fprintf(out, \\\"Number of inputed data : \\\");
fscanf(in, \\\"%d\\\",&n1);
fprintf(out,\\\"%d\\\",n1);
for(i=1;i<=n1;i++){
fscanf(in,\\\"%lf%lf\\\", &x[i],&f[i]);
fprintf(out, \\\"\n%lf %lf\\\", x[i],f[i]);
}
fprintf(out,\\\"\nNumber of outputed data : \\\");
fscanf(in, \\\"%d\\\", &m);
fprintf(out,\\\"%d\\\",m);

for(i=1;i<=m;i++) fscanf(in,\\\"%lf\\\",&s[i]);
fprintf(out,\\\"\n****** Solution ******\n\\\");
fprintf(out,\\\" NO. Point Value Error\\\");
lagrange();

for(i=1;i<=m;i++){
error=p[i]-fa(s[i]);
fprintf(out,\\\"\n %2d %7.3lf %7.5lf %15.6le\\\", i,s[i],p[i],error);
}
/* 파일닫기! */
fclose(in);
fclose(out);

}

void lagrange(void)
{

int i,j;
double sc, sp;
for (k=1;k<=m+3;k++){
sc=0;
for(i=1;i<=n1;i++){
sp=1.;
for(j=1;j<=n1;j++){
if(i!=j)
sp=sp*(s[k]-x[j])/(x[i]-x[j]);
}
sc+=(sp*f[i]);
}
p[k]=sc;
}
}

키워드

  • 가격1,000
  • 페이지수8페이지
  • 등록일2004.11.29
  • 저작시기2004.11
  • 파일형식압축파일(zip)
  • 자료번호#275539
본 자료는 최근 2주간 다운받은 회원이 없습니다.
청소해
다운로드 장바구니