高分求一个使用MKL解方程的源代码

[wxixi]

A为对称正定稀疏矩阵：
A=[ 2 0 0 1
0 8 3 0
0 3 6 0
1 0 0 9];
b=[1
2
3
4];

方程为：Ax=b

如何使用MKL的稀疏求解器Sparse Solver解上面方程？请教高手，源码给分！

[wxixi]

少说了个要求：要用C或C++语言。
高手们，快来解决这个实际问题啊！

[intel_cyu]

MKL目录下，有个例子程序 \examples\solver\source\pardiso_sym_c.c
把输入的数据改一下就可以。