안녕하세요

매트랩을 이번에 배워서 써먹어보려는데 뭔가 계산이 꼬이네요ㅠ

잘해서 멋진 그래프 그려보고싶어요 도와주세요!

 

식은 이렇고요

 

제가 매트랩 식을 짠거는

 

clc
clear

k=[0.598
1.65
3.1
3.87
2.33
1.29
0.462
0.328
10.17
0.0865
0.114]';

 

ramda= [1.77
0.5774
0.06743
0.006214
0.0004739
0.0000481
0.000005344
5.716E-07
1.036E-07
2.959E-08
75850000000]';

 

p=3800;
t=57600000;
ts=[0:10^7]';

 

for ts=0:10^7
    if 0<=ts<10^3
        K=0.1;
    P=p./200.*sum(k.*exp(-ramda.*ts));
    pfp=(1.+K).*P.*(t.+ts);
    elseif 10^3<=ts<10^7
        K=0.2;
    P=p./200.*sum(k.*exp(-ramda.*ts));
    pfp=(1.+K).*P.*(t.+ts);
           
    else
        K=0;
    P=p./200.*sum(k.*exp(-ramda.*ts));
    pfp=(1.+K).*P.*(t.+ts);
       
    end
    end   
pfp
plot(ts,pfp)

 

일케했어요

plot으로 ts를 x축 pfp를 y축으로 하려는데

 

어디가 잘못된걸까요?

for 문 안에 if를 넣은게 잘못일까요?