dear all,
I am solving two equation using newton method and meets problem in title. My whole procedure is below:
solve fi([setaf,fi],[dsetaf,dfi],solver=“UMFPACK64”)=
96 : int3d(th2)(cnewtondsetafnf*((6fi1^5-15fi1^4+10fi1^3)rawl+rawi-(6fi1^5-15fi1^4+10fi1^3)rawi)wf1(dx(setaf)+dy(setaf)+dz(setaf)))
97 : +int3d(th2)(cnewtondsetafnf*((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(cl-ci)setaf(30fi1^4-60fi1^3+30fi1^2)wf1(dx(fi1)+dy(fi1)+dz(fi1)))
98 : -int3d(th2)(cnewtondsetafnf*((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(cl-ci)setaf1(30fi1^4-60fi1^3+30fi1^2)wf1(dx(fi1)+dy(fi1)+dz(fi1)))
99 : +int3d(th2)(cnewtonnf((6fi1^5-15fi1^4+10fi1^3)hl+(1-(6fi1^5-15fi1^4+10fi1^3))hi)(dx(dsetaf)dx(setaf)+dy(dsetaf)dy(setaf)+dz(dsetaf)dz(setaf)))
100 : -int3d(th2)(cnewtondfi(30fi1^4-60fi1^3+30fi1^2)nf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(cl-ci)(log(setaf1)-log(277.15)+2)setaf)
101 : +int3d(th2)(cnewtondfi(30fi1^4-60fi1^3+30fi1^2)nf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(cl-ci)(log(setaf1)-log(277.15)+2)setaf1)
102 : -int3d(th2)(cnewtondfi(30fi1^4-60fi1^3+30fi1^2)nf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(sl-si)seffsetaf/277.15)
103 : +int3d(th2)(cnewtondfi(30fi1^4-60fi1^3+30fi1^2)nf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(sl-si)seffsetaf1/277.15)
104 : +int3d(th2)(cnewtondsetafnf*((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)3gamma(2+12fi1^2-12fi1)wf1(dx(fi1)+dy(fi1)+dz(fi1))fi/(widthrawl))
105 : -int3d(th2)(cnewtondsetafnf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)3gamma(2+12fi1^2-12fi1)wf1(dx(fi1)+dy(fi1)+dz(fi1))fi1/(widthrawl))
106 : +int3d(th2)(cnewtondsetafnf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)3gamma(2fi1+4fi1^3-6fi1^2)wf1(dx(fi)+dy(fi)+dz(fi))/(widthrawl))
107 : +int3d(th2)(cnewtondsetafnf((rawl-rawi)/(rawlrawl))(30fi1^4-60fi1^3+30fi1^2)wf1(dx(fi1)+dy(fi1)+dz(fi1))6gammawidth*(dx(fi)dx(fi1)+dy(fi)dy(fi1)+dz(fi)dz(fi1)))
108 : +int3d(th2)(cnewtondsetafnf((rawl-rawi)/(rawlrawl))(30fi1^4-60fi1^3+30fi1^2)wf1(dx(fi1)+dy(fi1)+dz(fi1))3gamma(2fi1+4fi1^3-6fi1^2)/width)
109 : +int3d(th2)(cnewtondsetafnf((rawl-rawi)/(rawlrawl))(120fi1^3-180fi1^2+60fi1)fiwf1(dx(fi1)+dy(fi1)+dz(fi1))3gammawidth(dx(fi1)dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1)))
110 : -int3d(th2)(cnewtondsetafnf((rawl-rawi)/(rawlrawl))(120fi1^3-180fi1^2+60fi1)fi1wf1(dx(fi1)+dy(fi1)+dz(fi1))3gammawidth(dx(fi1)dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1)))
111 : +int3d(th2)(cnewtondsetafnf((rawl-rawi)/(rawlrawl))(30fi1^4-60fi1^3+30fi1^2)wf1(dx(fi)+dy(fi)+dz(fi))3gammawidth*(dx(fi1)dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1)))
112 : +int3d(th2)(cnewtondsetafnf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(cl-ci)(setaf1-277.15)fi(120fi1^3-180fi1^2+60fi1)wf1(dx(fi1)+dy(fi1)+dz(fi1)))
113 : -int3d(th2)(cnewtondsetafnf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(cl-ci)(setaf1-277.15)fi1(120fi1^3-180fi1^2+60fi1)wf1(dx(fi1)+dy(fi1)+dz(fi1)))
114 : +int3d(th2)(cnewtondsetafnf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(cl-ci)(setaf1-277.15)(30fi1^4-60fi1^3+30fi1^2)wf1(dx(fi)+dy(fi)+dz(fi)))
115 : +int3d(th2)(cnewtondsetaf(120fi1^3-180fi1^2+60fi1)finf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(sl-si)seffwf1(dx(fi1)+dy(fi1)+dz(fi1)))
116 : -int3d(th2)(cnewtondsetaf(120fi1^3-180fi1^2+60fi1)fi1nf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(sl-si)seffwf1(dx(fi1)+dy(fi1)+dz(fi1)))
117 : +int3d(th2)(cnewtondsetaf(30fi1^4-60fi1^3+30fi1^2)nf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(sl-si)seffwf1*(dx(fi)+dy(fi)+dz(fi)))
118 : +int3d(th2)(cnewtondfiwf1*(dx(fi)+dy(fi)+dz(fi))/(vkf/(6rawlwidth*(sl-si))))
119 : -int3d(th2)(cnewtondfi(120fi1^3-180fi1^2+60fi1)finf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(cl-ci)(setaf1log(setaf1/277.15)+setaf1-277.15))
120 : +int3d(th2)(cnewtondfi*(120fi1^3-180fi1^2+60fi1)fi1nf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(cl-ci)(setaf1log(setaf1/277.15)+setaf1-277.15))
121 : -int3d(th2)(cnewtondfi*(120fi1^3-180fi1^2+60fi1)fi((rawl-rawi)/rawl)(3gammawidth*(dx(fi1)dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1))+(pow(fi1,2)pow(1-fi1,2)3gamma/width)))
122 : +int3d(th2)(cnewtondfi(120fi1^3-180fi1^2+60fi1)fi1((rawl-rawi)/rawl)(3gammawidth*(dx(fi1)dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1))+(pow(fi1,2)pow(1-fi1,2)3gamma/width)))
123 : -int3d(th2)(cnewtondfi(30fi1^4-60fi1^3+30fi1^2)((rawl-rawi)/rawl)(9gammawidth(dx(fi)dx(fi1)+dy(fi)dy(fi1)+dz(fi)dz(fi1))+fi((2fi1+4fi1^3-6fi1^2)3gamma/width)))
124 : +int3d(th2)(cnewtondfi*(30fi1^4-60fi1^3+30fi1^2)((rawl-rawi)/rawl)(fi1((2fi1+4fi1^3-6fi1^2)3gamma/width)))
125 : +int3d(th2)(cnewtondfi*(30fi1^4-60fi1^3+30fi1^2)((rawl-rawi)/rawl)(9gammawidth(dx(fi)dx(fi1)+dy(fi)dy(fi1)+dz(fi)dz(fi1))))
126 : -int3d(th2)(cnewtondfi3gammafi(2+12fi1^2-12fi1)/width)
127 : +int3d(th2)(cnewtondfi3gammafi1*(2+12fi1^2-12fi1)/width)
128 : -int3d(th2)(cnewton6gammawidth(dx(dfi)dx(fi)+dy(dfi)dy(fi)+dz(dfi)dz(fi)))
129 : +int3d(th2)(dsetafnf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)wf1(dx(setaf)+dy(setaf)+dz(setaf)))
130 : +int3d(th2)(dsetafnf*((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)3gamma(2fi1+4fi1^3-6fi1^2)wf1(dx(fi1)+dy(fi1)+dz(fi1))/(widthrawl))
131 : +int3d(th2)(dsetafnf*((rawl-rawi)/(rawlrawl))(30fi1^4-60fi1^3+30fi1^2)wf1(dx(fi1)+dy(fi1)+dz(fi1))3gammawidth*(dx(fi1)dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1)))
132 : +int3d(th2)(dsetafnf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(cl-ci)(setaf1-277.15)(30fi1^4-60fi1^3+30fi1^2)wf1(dx(fi1)+dy(fi1)+dz(fi1)))
133 : +int3d(th2)(dfi(30fi1^4-60fi1^3+30fi1^2)nf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(sl-si)seffwf1*(dx(fi1)+dy(fi1)+dz(fi1)))
134 : +int3d(th2)(nf*((6fi1^5-15fi1^4+10fi1^3)hl+(1-(6fi1^5-15fi1^4+10fi1^3))hi)(dx(dsetaf)dx(setaf1)+dy(dsetaf)dy(setaf1)+dz(dsetaf)dz(setaf1)))
135 : +int3d(th2)(dfiwf1(dx(fi1)+dy(fi1)+dz(fi1))/(vkf/(6rawlwidth*(sl-si))))
136 : -int3d(th2)(dfi*(30fi1^4-60fi1^3+30fi1^2)nf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(cl-ci)(setaf1log(setaf1/277.15)+setaf1-277.15))
137 : -int3d(th2)(dfi*(30fi1^4-60fi1^3+30fi1^2)nf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(sl-si)seffsetaf/277.15)
138 : +int3d(th2)(dfi*(30fi1^4-60fi1^3+30fi1^2)nf((6fi1^5-15fi1^4+10fi1^3)rawl+(1-(6fi1^5-15fi1^4+10fi1^3))rawi)(sl-si)seff)
139 : -int3d(th2)(dfi(30fi1^4-60fi1^3+30fi1^2)((rawl-rawi)/rawl)(3gammawidth(dx(fi1)dx(fi1)+dy(fi1)dy(fi1)+dz(fi1)dz(fi1))+(pow(fi1,2)pow(1-fi1,2)3gamma/width)))
140 : +int3d(th2)(dfi3gamma(2fi1+4fi1^3-6fi1^2)/width)
141 : +int3d(th2)(6gammawidth*(dx(dfi)dx(fi1)+dy(dfi)dy(fi1)+dz(dfi)dz(fi1)))
142 : -int2d(th2,upp,down)(dsetafN.x(nf((6fi1^5-15fi1^4+10fi1^3)hl+(1-(6fi1^5-15fi1^4+10fi1^3))hi)(dx(setaf1)))+dsetafN.y*(nf*((6fi1^5-15fi1^4+10fi1^3)hl+(1-(6fi1^5-15fi1^4+10fi1^3))hi)(dy(setaf1)))+dsetafN.z*(nf*((6fi1^5-15fi1^4+10fi1^3)hl+(1-(6fi1^5-15fi1^4+10fi1^3))hi)(dz(setaf1))));
143 : //-int2d(th2,upp)(dsetafN.y(nf*((6fi1^5-15fi1^4+10fi1^3)hl+(1-(6fi1^5-15fi1^4+10fi1^3))hi)(dy(setaf1))))
144 : //-int2d(th2,upp)(dsetafN.z(nf*((6fi1^5-15fi1^4+10fi1^3)hl+(1-(6fi1^5-15fi1^4+10*fi1^3))hi)(dz(setaf1))));
145 :
146 : //+on(upp,down,setaf=263.15);
147 : setaf1[]-=setaf[];
148 :
149 : err=setaf1[].linfty;
150 : if(err< 1e-2) cnewton=1;
151 : plot(setaf);
152 : //plot(fi);
153 :
154 : cout <<Ra << " " << step << " err ="<< err << endl;
155 : if( err < tol) break;
156 : }
157 : sizestack + 1024 =193416 ( 192392 )