|
.9e \n", i+1, p, fp);
// 오차한계보다 작을 때까지 지속
if ((fabs(p - p1) < tolerance)) {
printf("\nApproximate solution P = %15.9f \n", p);
printf("with Func(P) = %15.9f\n", fp);
printf("Number of iterations = %3d\n", i+1);
printf("Tolerance = %15.9e\n", tolerance);
ok = false;
} else {
i++;
p0 =
|
|
|
포트란으로 이분법 프로그램 소스입니다. 초보자용으로 교수님이 질문했을경우 답변하기 쉽게
코멘트 넣었구요. 누구나 쉽게 이해할수 있도록 do문 function 없습니다.
고급기술없이 쉽게 만들어서 A+맞은 소스입니다.
참고 하실분은 참
|
|
|
t=[n' xn(1:i)' fx' d']; %% Iteration별 결과 값을 나타내는데 xn은 다른 값들과 행렬을 맞추기 위하여 i까지만 돌게 지정해준다.
disp(result)
plot(xn) %% 수렴 값의 그래프
grid;
title('newton'); %% 그래프 제목을 표시
xlabel('계산 횟수'); %% 그래프 x축 제목을 표
|
|
|
기치 벡터 x를 입력
printf("\nx(0)를 입력해 주세요 : \n");
for(i=0;i<n;i++)
scanf("%lf", &x[i]);
Jacobi(n, Max, esp);
}
void Jacobi(int n, int Max, double esp)
{
int i, j, k;
double sum, err, next_x[50];
for(k=1;k<=Max;k++){
printf("\nk = %d \n", k);
for(i=0;i<n;i++){
sum=0.;
for(j=0;j<n
|
|
|
lt is... \n");
double x4=q4;
printf(" x4 = %lf\n",x4);
double x3=p3*x4 + q3;
printf(" x3 = %lf\n",x3);
double x2=p2*x3+q2;
printf(" x2 = %lf\n",x2);
double x1=p1*x2+q1;
printf(" x1 = %lf\n",x1);
}
|
|
|
1.5000
1.4844
0.0179
1.051
8
1.4688
1.4844
1.4766
-0.0034
0.528
9
1.4766
1.4844
1.4805
0.0073
0.263
10
1.4766
1.4805
1.4785
0.0019
0.135
11
1.4766
1.4785
1.4775
-0.0007
0.068
False-position of Bracketing method
function y = han3(x);
y = x^2 - 0.2333*x - 1.8391221;
end
>>falsi('han3', 0, 2, 0.0
|
|
|
- - - - - - - - - - - - - -");
do{
xi2 = xi - f(xi) / fp(xi);
error = xi2 - xi;
printf("\n %3d %6.4f %6.4f %6.4f %6.4f %6.4f ", count, xi, f(xi), fp(xi), f(xi2), error);
xi = xi2;
count++;
}while(fabs(error)>0.0001 || fabs(f(xi))>0.0001);
printf("\n\n\n\n -- - - - - - - - - - - - - - - - - -
|
|
|
|
#include<stdio.h>
#include<math.h>
#define maxDimension 6
double upper(int i,int j),lower(int i,int j);
double Amatrix[maxDimension][maxDimension];
double b[maxDimension],z[maxDimension],x[maxDimension];
double upperMatrix[maxDimension][maxDimension],lowerMatrix[maxDimension][m
|
|
|
);
printf("Root=%f \n Iteration=%d \n",xn ,it);
}
② 결과값
3. 느낀점
손으로 풀이를 할 때는 표도 그려야 하고 일일이 계산하는게 귀찮고 번거로웠는데 이렇게 프로그램을 짜서 해보니 구하는 값만 정확히 알면 쉽게 구해짐을 알 수 있었습니다. 이번 과
|
|
|
x;
hold on
plot(x,y,'r')
X = [5 6 10 14 16 20 22 28 28 36 38 5];
Y = [30 22 28 14 22 16 8 8 14 0 4 5];
for i=1:1:11
plot(X(i),Y(i),'o');
end
grid on
xlabel('x')
ylabel('y')
title('y=30.7396297 - 0.7719100*x & Data');
legend('y=30.7396297 - 0.7719100*x','data',0);
데이터와 회귀분석선 Graph
의 그래프
|
|
메뉴펼치기
