Question

matlab程序问题

建立M文件eq1．m
function dy=eq1(x,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=1/5*sqrt(1+y(1)^2)/(1-x);
取x0=0，xf=0．9999，建立主程序ff6．m如下：
x0=0，xf=0．9999
[x,y]=ode15s('eq1',[x0 xf],[0 0]);
plot(x,y(:,1),’b．')
hold on
y=0:0．01:2;
plot(1,y,’b*')

Answer 1

不知道你的程序是要什么样的结果，单纯从代码来看，就是全角和半角的问题。我已经修改如下：

建立M文件eq1．m

     function dy=eq1(x,y)

                  dy=zeros(2,1);

                  dy(1)=y(2);

                  dy(2)=1/5*sqrt(1+y(1)^2)/(1-x);

取x0=0，xf=0．9999，建立主程序如下： 

    x0=0;

       xf=0.9999;

      [x,y]=ode15s('eq1',[x0 xf],[0 0]);

      plot(x,y(:,1),'b.') 

      hold on

      y=0:0.01:2;

       plot(1,y,'b*')

得到的图像是

Answer 2

你的主程序输入中有些小数点等是全角字符
x0=0
xf=0.9999
[x,y]=ode15s('eq1',[x0 xf],[0 0]);
plot(x,y(:,1),'b.')
hold on
y=0:0.01:2;
plot(1,y,'b*')
画出的图很难看

Answer 3