|
2 = p1 - f(p1)*(p1-p0)/(f(p1)-f(p0));
p0 = p1;
p1 = p2;
n++;
printf("n= %d Pn= %f f(Pn)= %f \n",n,p2,f(p2));
}
printf("Solution is %f\n",p2);
return 0;
}
float f(float x)
{
return x*x*x - 2*x*x -5;
}
float absol(float x)
{
if(x>=0)return x;
else return (-1) * x;
}
(2)Result
2.3-3.b
(1) Source cod
|
|
|
()
{
double x1=0;
double x2=1;
double x3;
double x3o=0;
double f1=-1;
double f2=-0.12319;
while(1)
{
x3=x2-f2*(x2-x1)/(f2-f1);
printf("%lf\n",x3);
double f3=func(x3);
if (fabs((x3-x3o)/x3)*100 < 0.0005 )
break ;
if (fabs((x3-x3o)/x3)*100 > 0.0005)
{
x3o=x3;
x1=x2;
x2=x3;
f1=func(x1);
f2
|
|
|
일반적으로 기둥부재는 초기변형이 없는 직선이며, 하중은 기둥단면의 도심에 작용하여 편심이 없는 부재로 가정하고 이에 대한 거동을 살펴보았다. 그러나 이러한 이상적인 부재는 실제로는 존재하지 않는다. 일반적으로 부재의 제작이나
|
|
|
nitial guess를 0으로 주고 iteration 했을 때는 근을 찾지 못하고 종료하는 것을 보았다. 이로써 적절한 initial guess를 얻고 대입하는 것이 modified Secant Method를 사용하는데 있어서 또 하나의 과제라고 하겠다. 반대로 이 문제의 경우 friction factor의 방
|
|
|
개구간법이 튜튼랩슨, 시컨트 , 심플원포인트방법 사용
|
|
메뉴펼치기
