Ok mi scuso con lo staff 
Ecco qui:
codice:
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
int main(){
/* vengono definite le variabili*/
int n=0,n1=0,n2=0;
double a,b,x,x1,x2;
double fa,fb,fc,fx,fx1,fx2;
double errx,errx1,errx2;
printf("Occorrono due estremi a e b di un intervallo all'interno del\n");
printf("quale cercheremo un valore x soluzione dell'equazione f(x)=0\n");
printf("*** Attenzione! a e b devono essere tali che f(a)*f(b)<0 ***\n");
printf(" Estremo inferiore?\n");
scanf("%lf", &a);
fa = 200*sin(a)+25*cos(a)-160;
printf(" Estremo superiore?\n");
scanf("%lf", &b);
fb = 200*sin(b)+25*cos(b)-160;
printf("\n");
fc=fa*fb;
if(fc < 0.){
do {
x =(a + b) * 0.5;
fx = cos(x)-sin(x);
n=n+1;
if( fx == 0. ) {
printf("Trovo uno zero in %25.20g\n",x);
return false;
}
if( fa*fx < 0. ) {
b = x;
fb = fx;
}
else {
a = x;
fa = fx;
}
errx = fabs( (b-a) * 0.5 );
} while(errx >= 1.e-15);
printf("Trovo uno zero in %25.20g +/- %g \n",x,errx);
printf("Converge in %d passi\n",n-1);
}
/*else (fc>0);
{
printf("***Attenzione! Nell'intervallo non esistono due valori di x\n");
printf("per cui f(x)=0\n");
x =(a + b) * 0.5;
fx = cos(x)-sin(x);
n1=n1+1;
n2=n2+1;
if( fx == 0. ){
printf("Trovo uno zero in %25.20g\n",x);
return false ;
}
if( fa*fx < 0. ) {
b = x;
fb = fx;
do {
x1=(a+b)*0.5;
fx1=cos(x1)-sin(x1);
n1=n1+1;
if( fa*fx1 < 0. ) {
b = x1;
fb = fx1;
}
else {
a = x1;
fa = fx1;
}
errx1 = fabs( (b-a) * 0.5 );
} while(errx1 >= 1.e-15);
}
printf("Trovo uno zero in %25.20g +/- %g \n",x1,errx1);
printf("Converge in %d passi\n",n1);
if( fb*fx < 0. ){
a = x;
fa = fx;
do {
x2=(a+b)*0.5;
fx2=cos(x2)-sin(x2);
n2=n2+1;
if( fa*fx2 < 0. ) {
b = x2;
fb = fx2;
}
else {
a = x2;
fa = fx2;
}
errx2 = fabs( (b-a) * 0.5 );
} while(errx2 >= 1.e-15);
}
printf("Trovo uno zero in %25.20g +/- %g \n",x2,errx2);
printf("Converge in %d passi\n",n2);*/
system ("PAUSE");
}
/*} */
Rispondi quotando
spero di poter imparare stando a contatto con voi...
