function f=fc(y)
global D;
a=50;
b=1;
c=5;
f1=(1+D)*y(1)+3*y(1)*(2*(y(2)+y(1)^2)-y(1)^2);
f2=c*(15*y(1)*y(2)^2+10*y(2)*y(1)^3+y(1)^5)-2*(1+D)*y(1)^3;
f=[f1; f2];
clear;
clc;
s=[];
global D;
for D=0:0.01:1
s0=fsolve('fc',[rand;rand]);
s=[s;s0'];
end
把所有解都存在s里了
加一个简单的循环就可以实现了。