I showed results when I used vpasolve for this problem,but when I used it for another system of the same type I got this error Error using mupadengine/feval (line 157) MuPAD error: Error: Cannot differentiate the equation. [numeric::fsolve] Error

syms y22 y33 y5 y6 y7 y8 y9 t y(t) y1 v h=2*pi/(300) r=1-sqrt(3/7) s=1+sqrt(3/7) xn=0

y0=1.0 fn(t)=-y(t)+0.001*exp(sqrt(-1)*t) %y0=1.0 y11=(9995*sqrt(i))/(10000) y(v)=-y(v)+0.001*exp(sqrt(-1)*v) eq1=-y22-y0+2*y1+h^2/(10*s*(r - s)*(s^2 - 3*s + 2))*fn(s)-h^2/(10*r*(r - s)*(r^2 - 3*r + 2))*fn(r)-(h^2*(25*r + 25*s - 25*r*s - 28))/(30*(r - 1)*(s - 1))*fn(y1)+ (h^2*(5*r*s - 3))/(60*r*s)*fn(y0) -(h^2*(10*r + 10*s - 5*r*s - 17))/(60*(r - 2)*(s - 2))*fn(y22) eq2=-h*y33 -y0+y1-(h^2*(8*r - 19))/(60*s*(r - s)*(s^2 - 3*s + 2))*fn(s)+ (h^2*(8*s - 19))/(60*r*(r - s)*(r^2 - 3*r + 2))*fn(r) -(h^2*(73*r + 73*s - 65*r*s - 92))/(60*(r - 1)*(s - 1))*fn(y1)+ (h^2*(8*r + 8*s + 5*r*s - 19))/(120*r*s)*fn(y0) -(h^2*(82*r + 82*s - 45*r*s - 145))/(120*(r - 2)*(s - 2))*fn(y22) eq3=-h*y5-y0+y1+(h^2*(7*r - 4))/(60*s*(r - s)*(s^2 - 3*s + 2))*fn(s) -(h^2*(7*s - 4))/(60*r*(r - s)*(r^2 - 3*r + 2))*fn(r) -(h^2*(18*r + 18*s - 25*r*s - 14))/(60*(r - 1)*(s - 1))*fn(y1) -(h^2*(7*r + 7*s - 15*r*s - 4))/(120*r*s)*fn(y0)+ (h^2*(3*r + 3*s - 5*r*s - 2))/(120*(r - 2)*(s - 2))*fn(y22) eq4=-y6+(1 - r)*y0+r*y1+(h^2*r*(r^5 - 6*r^4 + 10*r^3 - 8*r + 3))/(60*s*(r - s)*(s^2 - 3*s + 2))*fn(s)+(h^2*(3*r - 8*s - 8*r*s + 12*r^2*s - 3*r^3*s + 3*r^2 - 7*r^3 + 2*r^4 + 3))/(60*(r - s)*(r - 2))*fn(r) -(h^2*r*(3*r + 7*s - 8*r*s - 8*r^2*s + 2*r^3*s + 3*r^2 + 3*r^3 - r^4 - 4))/(60*(s - 1))*fn(y1) -(h^2*(r - 1)*(5*r + 8*s - 27*r*s + 13*r^2*s - 2*r^3*s + 5*r^2 - 5*r^3 + r^4 - 3))/(120*s)*fn(y0)+(h^2*r*(r - 1)*(r + 2*s - 3*r*s - 3*r^2*s + 2*r^3*s + r^2 + r^3 - r^4 - 1))/(120*(r - 2)*(s - 2))*fn(y22) eq5=-h*y7-y0+y1+(h^2*(3*r^5 - 15*r^4 + 20*r^3 - 8*r + 3))/(60*s*(r - s)*(s^2 - 3*s + 2))*fn(s)+(h^2*(8*s - 60*r^2*s + 60*r^3*s - 15*r^4*s + 40*r^3 - 45*r^4 + 12*r^5 - 3))/(60*r*(r - s)*(r^2 - 3*r + 2))*fn(r)+ (h^2*(7*r + 7*s - 15*r*s + 20*r^3*s - 5*r^4*s - 10*r^4 + 3*r^5 - 4))/(60*(r - 1)*(s - 1))*fn(y1)+ (h^2*(8*r + 8*s - 35*r*s + 60*r^2*s - 30*r^3*s + 5*r^4*s - 20*r^3 + 15*r^4 - 3*r^5 - 3))/(120*r*s)*fn(y0)-(h^2*(2*r + 2*s - 5*r*s + 10*r^3*s - 5*r^4*s - 5*r^4 + 3*r^5 - 1))/(120*(r - 2)*(s - 2))*fn(y22) eq6=-y8+(1 - s)*y0+ s*y1 -(h^2*(3*s - 8*r - 8*r*s + 12*r*s^2 - 3*r*s^3 + 3*s^2 - 7*s^3 + 2*s^4 + 3))/(60*(r - s)*(s - 2))*fn(s) -(h^2*s*(s^5 - 6*s^4 + 10*s^3 - 8*s + 3))/(60*r*(r - s)*(r^2 - 3*r + 2))*fn(r) -(h^2*s*(7*r + 3*s - 8*r*s - 8*r*s^2 + 2*r*s^3 + 3*s^2 + 3*s^3 - s^4 - 4))/(60*(r - 1))*fn(y1) -(h^2*(s - 1)*(8*r + 5*s - 27*r*s + 13*r*s^2 - 2*r*s^3 + 5*s^2 - 5*s^3 + s^4 - 3))/(120*r)*fn(y0)+(h^2*s*(s - 1)*(2*r + s - 3*r*s - 3*r*s^2 + 2*r*s^3 + s^2 + s^3 - s^4 - 1))/(120*(r - 2)*(s - 2))*fn(y22) eq7=-h*y9-y0+y1-(h^2*(8*r - 60*r*s^2 + 60*r*s^3 - 15*r*s^4 + 40*s^3 - 45*s^4 + 12*s^5 - 3))/(60*s*(r - s)*(s^2 - 3*s + 2))*fn(s)-(h^2*(3*s^5 - 15*s^4 + 20*s^3 - 8*s + 3))/(60*r*(r - s)*(r^2 - 3*r + 2))*fn(r)+ (h^2*(7*r + 7*s - 15*r*s + 20*r*s^3 - 5*r*s^4 - 10*s^4 + 3*s^5 - 4))/(60*(r - 1)*(s - 1))*fn(y1)+ (h^2*(8*r + 8*s - 35*r*s + 60*r*s^2 - 30*r*s^3 + 5*r*s^4 - 20*s^3 + 15*s^4 - 3*s^5 - 3))/(120*r*s)*fn(y0) -(h^2*(2*r + 2*s - 5*r*s + 10*r*s^3 - 5*r*s^4 - 5*s^4 + 3*s^5 - 1))/(120*(r - 2)*(s - 2))*fn(y22) eq8=-h*y11-y0+y1-(h^2*(8*r - 3))/(60*s*(r - s)*(s^2 - 3*s + 2))*fn(s)+ (h^2*(8*s - 3))/(60*r*(r - s)*(r^2 - 3*r + 2))*fn(r)+(h^2*(7*r + 7*s - 15*r*s - 4))/(60*(r - 1)*(s - 1))*fn(y1)+ (h^2*(8*r + 8*s - 35*r*s - 3))/(120*r*s)*fn(y0) -(h^2*(2*r + 2*s - 5*r*s - 1))/(120*(r - 2)*(s - 2))*fn(y22)

[y1,y9,y8,y7,y6,y5,y33,y22]= vpasolve([eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8])