cp(1)=-0.4370;
cp(2)=0.1955;
ca(1)=0.3;
ca(2)=-1.5;
ce(1)=1;
ce(2)=-0.7;
w=0.0000000000000000001;
n=1;
c(1)=1;
d(1)=1;
p(1)=0.5353;
p(2)=0.0983;
p(3)=0.2938;

 
while sqrt(c.^2+d.^2)>=w
    
solve(' cp(n+1)=c(n)*ca(n+1)+d(n)*ce(n+1)','cp(n)=c(n)*ca(n)+d(n)*ce(n)','c(n)','d(n)');

我想请问下这种我前面已经给过了几个值的这种方程应该如何求解呢?(其中ca(n+1),ca(n),ce(n+1),ce(n),cp(n+1),cp(n+1)是通过循环得到的值)