|
#include<stdio.h>
#include<math.h>
#define TRUE 1
double function(double x) // 함수 지정
{
double y;
y=-256.25*x+1650;
return(y);
}
double function2(double x) // 미분함수 미리 지정
{
double y;
y=-256.25;
return(y);
} 없음
|
|
|
|
#include<stdio.h>
#include<math.h>
#define TRUE 1
#define Pi 3.141592 // 파이값 미리 설정
double function(double x) // 함수값 설정
{
double y;
y=-Pi*pow(x,3)+30*Pi*pow(x,2)-3000;
return(y);
}
double function2(double x) // 미분함수 설정
{
double y;
y=-3*Pi*pow(x,2)+6
|
|
|
- - - - - - - - - - - - - - - - - - - - - - - - - - \n");
}
double f(double x)
{
return (10+(3.592/x)/x)*(x-0.04267)-(0.082054*300);
}
double fp(double x)
{
return 10-(3.592/x)/x+((2*3.592*0.04267/x)/x)/x;
} 1. 이상기체 몰랄 볼륨 구하기
2. 실제기체 몰랄 볼륨 구하기 - 뉴틈랩슨 방법
|
|
|
뉴튼 랩슨법
clear all
clc
f='1/3*(sin(x).^3)+1/2*(cos(x).^2)-x'
df='(sin(x).^2)*(cos(x))-(cos(x))*(sin(x))-1'
n=-1; x0=0; err=0.0000001; num=1;
fprintf(' 횟수 해 오차 \n ');
fprintf(' x(n) x(n+1) |X(n+1)-X(n)| \n ');
while num>=err
x=x0;
a=eval(f); b=eval(df);
x1=x0-a/b;
num=abs(x1-x0);
x=x1;
|
|
|
보이기
System.out.println("\n-----------최종값-------------");//최종 출력
System.out.println("\n총 시행횟수 = "+iter+"\nx = "+str1+"\nea = "+ea);
}
}
2. 결과값 도출 1. 문제 및 결과
(a) 황금분할탐색법(xl = -2, xu = 4, εs = 1%)
(b) 뉴튼(Newton) 법(X0 = 3, εs = 1%)
|
|
메뉴펼치기
