c语言多元方程组编程,求大家帮我把这个非线性方程组编程一个C语言。非常感谢！...

a1,a2,a3,b1,b2,b3,k1-k23是已知的，求解N1-N29一共29个未知数，就是29个方程解29个未知数，求大神帮忙，我的C语言忘完了

a1=0.3;a2=0.07;a3=0.015;b1=0.465;b2=0.05;b3=0.1;，

k1=659.1868765;

k2=18189.20742;

k3=6012.066351;

k4=30.06055264;

k5=71.99830518;

k6=139.8574962;

k7=268.7249903;

k8=13030895032409;

k9=5324.185271;

k10=7.569522933;

k11=38.3802878;

k12=12.59153449;

k13=6.684150735;

k14=4.193535028;

k15=4.065965375;

k16=101.9568377;

k17=4.664740682;

k18=1243.435699;

k19=1.042678642;

k20=1907.157919;

k21=89917.24523;

k22=239542.0967;

k23=31859855.35;

f(1)=N(1)+N(2)+N(3)+N(4)+N(5)+N(6)+N(7)+N(8)+N(9)+N(10)+N(11)+N(12)+N(13)+N(14)+N(15)+N(16)+N(17)+N(18)+N(19)+N(20)+N(21)+N(22)+N(23)+N(24)+N(25)+N(26)+N(27)+N(28)+N(29)-1=0;

f(2)=a3*(N(1)+N(7)+N(8)+N(9)+N(16)+N(17)+N(20)+2*N(23)+2*N(24)+N(25)+N(26)+2*N(27))-a1*(N(3)+N(15)+N(19)+N(22))=0;

f(3)=a3*(N(2)+N(10)+2*N(11)+N(12)+6*N(13)+7*N(14)+N(18)+N(21)+3*N(23)+N(24)+N(25))-a2*(N(3)+N(15)+N(19)+N(22))=0;

f(4)=a3*(0.5*N(4)+N(7)+2*N(28)+3*N(9)+N(10)+N(11)+3*N(12)+N(13)+12*N(14)+N(15)+N(24)+2*N(25)+N(26)+N(27)+2*N(28)+3*N(29))-b1*(N(3)+N(15)+N(19)+N(22))=0;

f(5)=a3*(0.5*N(5)+N(16)+2*N(17)+N(18)+N(19)+N(26)+N(27)+N(28)+N(29))-b2*(N(3)+N(15)+N(19)+N(22))=0;

f(6)=a3*(0.5*N(6)+2*N(20)+N(21)+N(22))-b3*(N(3)+N(15)+N(19)+N(22))=0;

f(7)=N(7)-k1*N(1)*N(4)=0;

f(8)=N(8)-k2*N(1)*N(4)^2=0;

f(9)=N(9)-k3*N(1)*N(4)^3=0;

f(10)=N(10)-k4*N(2)*N(4)=0;

f(11)=N(11)-k5*N(2)^2*N(4)=0;

f(12)=N(12)-k6*N(2)*N(4)^3=0;

f(13)=N(13)-k7*N(2)^6*N(4)=0;

f(14)=N(14)-k8*N(2)^7*N(4)^12=0;

f(15)=N(15)-k9*N(3)*N(4)=0;

f(16)=N(16)-k10*N(1)*N(5);

f(17)=N(17)-k11*N(1)*N(5)^2;

f(18)=N(18)-k12*N(2)*N(5)=0;

f(19)=N(19)-k13*N(3)*N(5)=0;

f(20)=N(20)-k14*N(1)*N(6)^2=0;

f(21)=N(21)-k15*N(2)*N(6)=0;

f(22)=N(22)-k16*N(3)*N(6)=0;

f(23)=N(23)-k17*N(2)^3*N(1)^2=0;

f(24)=N(24)-k18*N(1)^2*N(2)*N(4)=0;

f(25)=N(25)-k19*N(1)*N(2)*N(4)^2=0;

f(26)=N(26)-k20*N(1)*N(4)*N(5)=0;

f(27)=N(27)-k21*N(1)^2*N(4)*N(5)=0；

f(28)=N(28)-k22*N(1)^2*N(4)^2*N(5)=0;

f(29)=N(29)-k23*N(1)^2*N(4)^3*N(5)=0;