Escalonamento de matriz

Gente estou com problemas na hora de depurar esse código e não sei onde está o erro, gostaria da ajuda de vocês para saber onde estou errando

a mensagem de erro é "segmentaion fault"

#include <stdio.h>

int main()
{
int i,j,k,N;

printf("Qual e a ordem do sistema?\n");
scanf("%d", &N);
float B[N][N+1];
for(i=1;i<N+1;i++)
{
for(j=1;j<N+2;j++)
{
printf("Digite o valor do coeficiente A[%d", i);
printf("][%d",j);
printf("] da matriz gaussiana:");
scanf("%f*", &B[i][j]);
}
}
for(i=1;i<=N;i++)
{
for(j=i+1;j<=N;j++)
{
for(k=N+1;k>=i;i--)
{
B[j][k]=B[i][k]*B[j][i]/B[i][i]-B[i][k]; -B[i][k] *\para substituir a linha diretamente/*
}
}
}
for(i=1;i<N+1;i++)
{
for(j=1;j<N+1;j++)
{
printf("O elemento de A[%d", i);
printf("][%d",j);
printf("] escalonado é%f*", B[i][j]);
}
}

return 0;
}

O meu objetivo e escalonar uma matriz mas ta com problema para 2x2 em diante e nãi sei pq, pois ao meu ver está correto

Por favor, edite sua postagem original e coloque a tag [code] imeditamanete antes do seu código, e [/code] imediatamente depois. Com isso, vai facilitar muito na hora de ler -- especialmente porque o fórum normalmente considera que [i] é um sinal para iniciar texto em itálicos, e isso não acontece quando você marca devidamente o início e o final do código fonte.

E lmebre-se que, em C, se você diz que um array tem N elementos, os índices desse elementos variam de 0 a N-1. Veja o programa abaixo.

int fib[10];

int i;



fib[0]=fib[1]=1;



/* ‘i’ conta de 2 a 9 */

for(i=2; i<10; i++)

  fib[i]=fib[i-1]+fib[i-2];



/*

  ‘i’ conta de 0 a 9, mas na hora de se imprimrir o número ordinal,

  usa-se ‘i+1’, até porque não existe “zerésimo”.

*/

for(i=0; i<10; i++)

  printf("O %dº número da sequencia de Fibonacci é %d.\n", i+1, fib[i]);