Разыменование пустого указателя

#include <stdio.h>
#include <locale.h>
#include <math.h>
#include <malloc.h>
 
int main(void)
{
    setlocale(LC_ALL, "ru");
    int n;
    printf("Введите n(1-10):\nn=");
    scanf_s("%i", &n);
    if ((n <= 0) || (n > 10))
    {
        printf("\nНеверно введено n");
    }
    else
    {
        double** a = (double**)malloc(n * sizeof(double*)), * sum = (double*)calloc(n, sizeof(double)), p, q, r;
        int* k = (int*)malloc(n * sizeof(int)), jp = -1, rj, s = 0, t = 0;
        if ((a) && (sum) && (k))
        {
            for (int i = 0; i < n; i++)
            {
                a[i] = (double*)malloc(n * sizeof(double));
            }
            printf("Введите матрицу:\n");
            for (int i = 0; i < n; i++)
            {
                for (int j = 0; j < n; j++)
                {
                    scanf_s("%lf", *(a + i) + j);
                }
            }
            printf("Введите начало диапазона:\np=");
            scanf_s("%lf", &p);
            printf("Введите конец диапазона:\nq=");
            scanf_s("%lf", &q);
            if (p >= q)
            {
                printf("\nВведены неверные p и q");
            }
            else
            {
                printf("Введённый диапазон:\n[%g,%g]\n", p, q);
                for (int j = 0; j < n; j++)
                {
                    for (int i = 0; i < n; i++)
                    {
                        if (jp == -1)
                        {
                            if ((*(*(a + i) + j) >= p) && (*(*(a + i) + j) <= q))
                            {
                                jp = j;
                            }
                        }
                    }
                }
                printf("\nИсходная матрица:\n");
                for (int i = 0; i < n; i++)
                {
                    for (int j = 0; j < n; j++)
                    {
                        printf("%g ", *(*(a + i) + j));
                    }
                    printf("\n");
                }
                if (jp != -1)
                {
                    printf("Номер первого столбца, содержащего хотя бы один элемент из диапозона [%g,%g]: %i", p, q, jp);
                }
                else
                {
                    for (int i = 0; i < n; i++)
                    {
                        for (int j = 0; j < n; j++)
                        {
                            if (**a == *(*(a + i) + j))
                            {
                                s++;
                            }
                        }
                    }
                    if (s == n * n)
                    {
                        printf("\nВ данной матрице нельзя упорядочить строки по возрастанию сумм элементов этих строк");
                    }
                    else
                    {
                        for (int i = 0; i < n; i++)
                        {
                            for (int j = 0; j < n; j++)
                            {
                                *(sum + i) += *(*(a + i) + j);
                            }
                            *(k + i) = i;
                        }
                        for (int i = 0; i < n; i++)
                        {
                            for (int j = 0; j < n - 1; j++)
                            {
                                if (*(sum + j) > * (sum + j + 1))
                                {
                                    r = *(sum + j);
                                    *(sum + j) = *(sum + j + 1);
                                    *(sum + j + 1) = r;
                                    rj = *(k + j);
                                    *(k + j) = *(k + j + 1);
                                    *(k + j + 1) = rj;
                                }
                            }
                        }
                        printf("\nМатрица, в которой упорядочены строки по возрастанию сумм элементов этих строк:\n");
                        for (int i = 0; i < n; i++)
                        {
                            for (int j = 0; j < n; j++)
                            {
                                printf("%g ", *(*(a + *(k + i)) + j));
                            }
                            printf("\n");
                        }
                    }
                }
            }
        }
    }
}