t=0.1; h=0.01; eps=0.000001; k=100; A=115.9375006; B=118.7349801; C=121.2134001; Da=5; Db=6; Dc=7; i=dnetn(5,eps,t,h,x,k); printf("\n"); printf("i=%d\n",i); printf("\n"); for(i=0;i<=4;i ) printf("x(%d)=%f\n",i,x[i]); printf("\n"); } void dnetnf(x,y,n) int n; double x[],y[]; { y[0]=A*cos(x[0])-x[3]; y[1]=B*cos(x[1])-x[3]; y[2]=C*cos(x[2])-x[3]; y[3]=tan(x[0])-x[0]-tan(x[1]) x[1]-(Da-Db)/(x[3]*cos(x[4])); y[4]=tan(x[1])-x[1]-tan(x[2]) x[2]-(Db-Dc)/(x[3]*cos(x[4])); n=n; return; } |