• 0
Sign in to follow this  
PauloHM

Algoritmo interpolador de lagrange

Question

O código abaixo realiza interpolação de uma função através do método de lagrange,de fórmula genérica:

http://upload.wikimedia.org/math/a/7/8/a78...875653bbd2a.png

Minha solução compila sem apresentar erros, porém , o resultado não é exibido em tela. Se alguém encontrar alguma irregularidade, fico grato.

#include <stdio.h>
#include <stdlib.h>

int main(int argc, char *argv[])
{
  
  int tamanho,cont,i,j;
     
     printf("Insira o numero de valores conhecidos da funcao que deseja interpolar:\t");
     scanf("%d",&tamanho);
     printf("\n"); 
     
  float vetfx[tamanho],vetx[tamanho],x,pdtparcial,pdtfinal,fx;
     
     printf("Valores de x pre existentes:\n");
     
     for(cont=0;cont<tamanho;cont++){
        printf("Insira o valor de x%d :\t",cont);
        scanf("%f",&vetx[cont]);
     }                                
     printf("\n");   
        
          
     printf("Valores de f<x> pre existentes:\n");    
     
     for(cont=0;cont<tamanho;cont++){
        printf("Insira o valor de y%d :\t",cont);
        scanf("%f",&vetfx[cont]);
     }                                
     printf("\n");   

     printf("Insira o valor de x para o qual deseja obter f<x>: ");
     scanf("%f",&x);
     printf("\n");
     
     fx=0;        
     pdtfinal=1;
          
     for(i=0;i<tamanho;i++){
       for(j=0;j<=tamanho;j++){
          if(j=i){j++;}
          pdtparcial=0;
          
          pdtparcial=(x-vetx[j])/(vetx[i]-vetx[j]);          
          
          if(j!=tamanho){ pdtfinal*=pdtparcial; }
       }
       
       fx=fx+vetfx[i]*pdtfinal;
        
       pdtfinal=0;                                            

     } 

     printf("O valor de f<%f> e : %.4f ",x,fx);
     printf("\n\n");
     

  system("PAUSE");    
  return 0;
}

[]'s

Share this post


Link to post
Share on other sites

2 answers to this question

Recommended Posts

  • 0
O código abaixo realiza interpolação de uma função através do método de lagrange,de fórmula genérica:

http://upload.wikimedia.org/math/a/7/8/a78...875653bbd2a.png

Minha solução compila sem apresentar erros, porém , o resultado não é exibido em tela. Se alguém encontrar alguma irregularidade, fico grato.

#include <stdio.h>
#include <stdlib.h>

int main(int argc, char *argv[])
{
  
  int tamanho,cont,i,j;
     
     printf("Insira o numero de valores conhecidos da funcao que deseja interpolar:\t");
     scanf("%d",&tamanho);
     printf("\n"); 
     
  float vetfx[tamanho],vetx[tamanho],x,pdtparcial,pdtfinal,fx;
     
     printf("Valores de x pre existentes:\n");
     
     for(cont=0;cont<tamanho;cont++){
        printf("Insira o valor de x%d :\t",cont);
        scanf("%f",&vetx[cont]);
     }                                
     printf("\n");   
        
          
     printf("Valores de f<x> pre existentes:\n");    
     
     for(cont=0;cont<tamanho;cont++){
        printf("Insira o valor de y%d :\t",cont);
        scanf("%f",&vetfx[cont]);
     }                                
     printf("\n");   

     printf("Insira o valor de x para o qual deseja obter f<x>: ");
     scanf("%f",&x);
     printf("\n");
     
     fx=0;        
     pdtfinal=1;
          
     for(i=0;i<tamanho;i++){
       for(j=0;j<=tamanho;j++){
          if(j=i){j++;}
          pdtparcial=0;
          
          pdtparcial=(x-vetx[j])/(vetx[i]-vetx[j]);          
          
          if(j!=tamanho){ pdtfinal*=pdtparcial; }
       }
       
       fx=fx+vetfx[i]*pdtfinal;
        
       pdtfinal=0;                                            

     } 

     printf("O valor de f<%f> e : %.4f ",x,fx);
     printf("\n\n");
     

  system("PAUSE");    
  return 0;
}

[]'s

Alguém ? =S

Share this post


Link to post
Share on other sites
  • 0

//Verifique as linhas marcadas em vermelho. Seu código agora está 100%!!!

 

#include <stdio.h>
#include <stdlib.h>

int main(int argc, char *argv[])
{
  
  int tamanho,cont,i,j;
     
     printf("Insira o numero de valores conhecidos da funcao que deseja interpolar:\t");
     scanf("%d",&tamanho);
     printf("\n"); 
     
  float vetfx[tamanho],vetx[tamanho],x,pdtparcial,pdtfinal,fx;
     
     printf("Valores de x pre existentes:\n");
     
     for(cont=0;cont<tamanho;cont++){
        printf("Insira o valor de x%d :\t",cont);
        scanf("%f",&vetx[cont]);
     }                                
     printf("\n");   
        
          
     printf("Valores de f<x> pre existentes:\n");    
     
     for(cont=0;cont<tamanho;cont++){
        printf("Insira o valor de y%d :\t",cont);
        scanf("%f",&vetfx[cont]);
     }                                
     printf("\n");   

     printf("Insira o valor de x para o qual deseja obter f<x>: ");
     scanf("%f",&x);
     printf("\n");
     
     fx=0;        
     pdtfinal=1;
          
     for(i=0;i<tamanho;i++){
       for(j=0;j<=tamanho;j++){
          if(j==i){j++;}
          pdtparcial=0;
          
          pdtparcial=(x-vetx[j])/(vetx-vetx[j]);          
          
          if(j!=tamanho){ pdtfinal*=pdtparcial; }

       }
       
       fx=fx+vetfx*pdtfinal;
        
      // pdtfinal=0;                                            
       pdtfinal=1;                                            
     } 

     printf("O valor de f<%f> e : %.4f ",x,fx);
     printf("\n\n");
     

  system("PAUSE");    
  return 0;
}

1 minuto atrás, rfabriciors disse:

//Verifique as linhas marcadas em vermelho. Seu código agora está 100%!!!

 

#include <stdio.h>
#include <stdlib.h>

int main(int argc, char *argv[])
{
  
  int tamanho,cont,i,j;
     
     printf("Insira o numero de valores conhecidos da funcao que deseja interpolar:\t");
     scanf("%d",&tamanho);
     printf("\n"); 
     
  float vetfx[tamanho],vetx[tamanho],x,pdtparcial,pdtfinal,fx;
     
     printf("Valores de x pre existentes:\n");
     
     for(cont=0;cont<tamanho;cont++){
        printf("Insira o valor de x%d :\t",cont);
        scanf("%f",&vetx[cont]);
     }                                
     printf("\n");   
        
          
     printf("Valores de f<x> pre existentes:\n");    
     
     for(cont=0;cont<tamanho;cont++){
        printf("Insira o valor de y%d :\t",cont);
        scanf("%f",&vetfx[cont]);
     }                                
     printf("\n");   

     printf("Insira o valor de x para o qual deseja obter f<x>: ");
     scanf("%f",&x);
     printf("\n");
     
     fx=0;        
     pdtfinal=1;
          
     for(i=0;i<tamanho;i++){
       for(j=0;j<=tamanho;j++){
          if(j==i){j++;}
          pdtparcial=0;
          
          pdtparcial=(x-vetx[j])/(vetx-vetx[j]);          
          
          if(j!=tamanho){ pdtfinal*=pdtparcial; }

       }
       
       fx=fx+vetfx*pdtfinal;
        
      // pdtfinal=0;                                            
       pdtfinal=1;                                            
     } 

     printf("O valor de f<%f> e : %.4f ",x,fx);
     printf("\n\n");
     

  system("PAUSE");    
  return 0;
}

Ah... obrigado pelo código... salvou minha vida 😉

Share this post


Link to post
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.

Guest
Answer this question...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

Sign in to follow this