Hi all, Im working on a code that simulates a stock price over two years and then compares it to a strike price (call options). I want to simulate the stock price 3000 times for each starting value (30 to 150), then take the average payout (for that particular staring price) and plot that vs the starting price. The idea is that if the final price is higher than the strike price (starting price) you get a payout equal to the difference between the final and strike price, if not you get 0.
I keep getting an error (Error: Index in position 2 exceeds array bounds (must not exceed 730)) It gives the error in the line
S(k+1, n) = S(k, n) + S(k, n).*r.*dt + sqrt(dt).*sigma.*S(k, n).*eta(k, n);
Here is my code
r = 0.05;T = 2;dt = 1/365;sigma = sqrt(r);rng('shuffle')eta = randn(1, 365*T);n_iterations = 3000;S = zeros(365*T, n_iterations); for S_0 = 30 : 150 S(1, :) = S_0 ; n = 1 : n_iterations; for k = 1 : 365*T S(k+1, n) = S(k, n) + S(k, n).*r.*dt + sqrt(dt).*sigma.*S(k, n).*eta(k, n); end strike = S(1, :); price = S((365*T + 1), :); eurocallreturn(strike, price); avereturn = mean(eurocallreturn); hold on plot(strike, avereturn) xlabel('Price') ylabel('Return')end plot([50:130], exp(0.05*2)*blsprice([50:130], 90, 0.05, 2, sqrt(0.05))) function R = eurocallreturn(strike, price)if price > strike R = price - strike; else if price <= strike R = 0; endend end %What am I doing wrong? And how do I get matlab to
% store all the outputs for avereturn (for each starting value) in a vector?
%Thanks in advance
Best Answer