I make some changes:
% split the solution vector
pas=size(A11);
v1=v(1:pas);
v2=v(pas(1)+1:2*pas(1));
v3=v(2*pas(1)+1:3*pas(1));
v4=v(3*pas(1)+1:4*pas(1));
% boundary conditions
for i=1:pas
vv1(i,1)=1;
vv2(i,1)=0;
vv3(i,1)=0;
vv4(i,1)=0;
end
vv1(:,2)=v1;
vv2(:,2)=v2;
vv3(:,2)=v3...
Hi Greg,
Thanks for your reply.
I think my issues start when I try to split the big vector v, into 4 small vectors (one for each equation:v1,v2,v3,v4).
After that, the boundary conditions(vv1,vv2,vv3,vv4) is not correctly implemented (it's need to be created dynamically).
In conclusion, there...