restart; eq1 := 531062-S/(70*365)-(.187*(1/365))*(H+C+C1+C2)*S/N = 0; eq2 := (4/365*(T+C))*S/N-(.187*(1/365))*(H+C+C1+C2)*T/N-(1/(70*365)+1/(5*365))*T = 0; eq3 := (.187*(1/365))*(H+C+C1+C2)*S/N-(4/365)(T+C)*H/N-(1/(70*365)+1/(4*365))*H = 0; eq4 := (.187*(1/365))*(H+C+C1+C2)*T/N+(4/365*(T+C))*H/N-(1/(70*365)+3/(8*365)+.2*(1/365)+.1)*C = 0; eq5 := .1*C-(1/(70*365)+1/(4*365)+1/60+1)*C1 = 0; eq6 := (1/60)*C1-(1/(70*365)+1/(4*365)+1/210+1)*C2 = 0; eq7 := C1-(1/(70*365)+1/60+0.1e-2)*CT1 = 0; eq8 := C2-(1/(70*365)+1/210+(1/9)*(0.1e-2*7))*CT2+(1/60)*CT1 = 0; eq9 := N-S-T-H-C-C1-C2-CT1-CT2 = 0; soln := solve({eq1, eq2, eq3, eq4, eq5, eq6, eq7, eq8, eq9}, {C, C1, C2, CT1, CT2, H, N, S, T})

@Axel Vogt 

If I were to give values for eta_1 and eta_2, do I have a shot at a solution?

@abscissa 

@Axel Vogt 

This is actually my MATLAB code, just to show the equations. My Maple code does not use the dot, and I have numbers in place of Lambda and gamma.

So what do you think?

@acer 

That's a moot point now - I have given up on this endeavor because neither Maple nor Mathematica were able to handle it!

@Carl Love 

Wow. Maple can do a lot of things! This makes my work so much easier.

Thanks again

@Carl Love 

Thanks! It works now. 

Is there a way I can evaluate the whole Jacobian at a point?