?? calculateyfunc.m
字號:
function [yy]=calculatey(eta,sita,m1,m2)
eps1=10^-6;
if eta==0
yy=sqrt(abs(-1/(2*cos(sita))));
else
y0=initialyfunc(eta,m1);
[f1,d1]=oadfunc(y0,eta,sita,m1,m2);
k=1;
y(1)=y0;
y(2)=y(1)-f1/d1;
while abs(f1)>eps1
u=1;
k=k+1;
[f1,d1]=oadfunc(y(k),eta,sita,m1,m2);
y(k+1)=y(k)-u*f1/d1;
[f2,d2]=oadfunc(y(k+1),eta,sita,m1,m2);
while abs(f2)>=abs(f1)
u=u/2;
y(k+1)=y(k)-u*f1/d1;
[f2,d2]=oadfunc(y(k+1),eta,sita,m1,m2);
end
end
yy=y(k);
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -